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Chapter I. Introduction 

Polymers have been widely used as wiring electrical insulation materials in space/air-craft. The 
dielectric properties of insulation polymers can change over time, however, due to various aging 
processes such as exposure to heat, humidity and mechanical stress. Therefore, the study of polymers 
used in electrical insulation of wiring is important to the aerospace industry due to potential loss of life 
and aircraft in the event of an electrical fire caused by breakdown of wiring insulation. 

Part of this research is focused on studying the mechanisms of various environmental aging process of 
the polymers used in electrical wiring insulation and the ways in which their dielectric properties change 
as the material is subject to the aging processes. The other part of the project is to determine the feasibility 
of a new capacitive nondestructive testing method to indicate degradation in the wiring insulation, by 
measuring its permittivity. 

1 . Motivation 

Dielectric wiring insulation is used to separate electrical conductors by preventing the flow of charge 
between wires. Insulation materials function to maintain a continuous and specified value of permittivity 
over a specified range of electromagnetic field frequency and strength. Another essential property of 
wiring insulation is the dielectric strength, a field at which the material fails to resist the flow of current 
and arcing occurs. The dielectric properties of potential wiring insulation materials are always carefully 
considered to guarantee that the selected materials satisfy requirements of the operating environment. 

Both the dielectric permittivity and dielectric strength of wire insulation may change over time, 
however, due to various degradation processes such as thermal aging, moisture exposure and mechanical 
degradation. For example, wiring may be improperly installed and maintained, increasing the risk of 
damage due to heat, moisture and chafing [1]. Such damage mechanisms may act acutely, or act to ‘age’ 
the insulation material over many cycles of aircraft operation. These mechanisms by which wire systems 
insulation may be degraded produce what are known as a ‘soft’ faults, which act to modify the impedance 
of the affected region of the coated wire structure, when viewed as a transmission line, rather than a ‘hard’ 
fault such as an open or short in the conductor itself. It has been reported [2] that aircraft suffer from 
undiagnosed wiring degradation which may cause short-circuiting, fire and loss of control function. 
According to Captain Jim Shaw, manager of the in-flight fire project for the United States Air Line Pilots 
Association (ALP A), there are on average three fire and smoke events in jet transport aircraft each day in 
USA and Canada alone, and the vast majority are electrical. It was presented in Air Safety Week, 2001, 
that aircraft were making emergency landings, suffering fire damage to the point of being written off etc, 
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at the rate of more than one a month based on the experience of the previous few months. These issues 
remain a concern for new aircraft. 

Motivated by these concerns, the contribution of this work is to explore and record change in dielectric 
properties of wire insulation due to various degradation processes. 

2. History of wiring insulation 

Table 1 shows wiring insulation materials applied in commercial aircraft since the 1960s [1]. PVC 
(polyvinyl chloride) and Nylon were the main insulation materials from the 1960s to the 1980s. However, 
in the next decades, PI (polyimide) was almost the only wiring insulation polymer used in the listed 
airplanes. After the 1900s, another two materials, TKT 1 (Teflon -Kapton -Teflon) and Tefzel® ETFE 
(ethylene-tetrafluoroethylene), have been widely used. This work focuses on three polymers: PI, PTFE 
and ETFE. More detailed information about these polymers will be introduced in Chapters II, section 3. 

3. Technical approach 
Permittivity 

The permittivity is a parameter that indicates the relative charge storage capability of dielectrics in the 
presence of an electric field. In general, permittivity is complex, denoted s* = s' — je ” . Complex 
permittivity measurements have been made on the sample materials investigated here, before and after 
degradation, to explore the changes in permittivity and dielectric relaxations in response to degradation. 
Two instruments were employed to measure complex permittivity of the polymers. The first one is a 
Novocontrol Spectrometer, which is capable of measurement over frequency range from 1 pHz to 3 GHz. 
A temperature-controlled sample cell also permits measurements at temperatures from -200 °C to 400 °C. 
The other one is an Agilent E4980A LCR meter coupled with a 16451 dielectric test fixture, which is 
available from 20 Hz to 2 MHz at room temperature. 

Breakdown voltage 

Another essential property of dielectric insulators is the dielectric breakdown voltage, the point at 
which the applied voltage causes current flow in a device (transistor, capacitor etc) to increase 
uncontrollably. Breakdown in a capacitor results in the replacement of a reactive insulating component by 
either a low-resistance short circuit or open circuit, usually with disastrous consequences as far as the 
overall circuit function is concerned. The probability of its occurrence must therefore be kept to an 


1 Teflon is a trade name for PTFE (polytetrafluoroethylene) and Kapton is a trade name for polyimide. 
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absolute minimum. Dielectric breakdown of insulation polymers before and after degradation processes 
have been measured by a DIELECTRIC RIGIDITY 6135 which can supply voltage up to 60 kV.. 

Thermal exposure 

Thermal exposure can significantly influence properties of polymers by changing microstructure, phase 
morphology, chemical composition, etc. The effect of thermal exposure in air on the permittivity of PI, 
PTFE and PI has been explored, which will be discussed in Chapters II, III, and IV, respectively. 


Table 1 History of wiring insulation used in commercial aircraft. 


Wiring insulation 
material 

Applied 

years 

Applied aircraft 

PVC/ Nylon 

1960-1980 

707, 727, 737, and DC-8 

Slash 6 

1965-1985 

DC-9 

Poly-X and Stilan 

1970-1980 

747, DC-10/MD-1 1 

PI 

1970 later 

727, 737, 757, 767, MD-80/-90, DC-10/MD-1 1, Lockheed L- 
1011, Airbus 

Tefzel 

1975 later 

727, 737, 757, 767, MD-80/-90, DC-10/MD-1 1 

Tefzel/PI 

1980-1990c 

747 

TKT 

1990 later 

737, 757, MD-80/-90, DC-10/MD-1 1 


Water and saline exposure 

There are several physical consequences of water absorption to wire insulation material including 
plasticization, swelling, and changes in dielectric properties. Even though polyimide has very good 
electrical and physical properties, it is very susceptible to humidity, which can give rise to cracks in the 
insulation and cause electrical malfunctions. In response to the concern that aircraft which serve in navy 
are exposed to sea water, the effect of water and saline exposure on dielectric properties of polyimide has 
been studied. Effect of saline exposure on insulating properties of PI will be presented in Chapter II. 

Mechanical stress 

During cycles of aircraft operation and due to improper installation, wiring insulation materials may be 
exposed to mechanical stress, which can result in structural changes and consequently influence the 
dielectric properties of the insulation. Given this concern, the influence of mechanical strain on the 
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permittivity of PTFE is investigated in Chapter III. A system capable of measuring dielectric permittivity 
while a polymer sample is simultaneously under tensile strain is designed and applied. 

Material characterization 

Thermal analysis instruments, such as a TGA (thermogravimetric analyzer) Q50 instrument, a DMA 
(dynamic mechanical analyzer) Q800 instrument and a DSC (differential scanning calorimeter) Q20 
instrument, are used to investigate thermal properties of the polymers. TGA uses heat to induce chemical 
and physical changes in materials and performs a corresponding measurement of mass change as a 
function of temperature or time. In some advanced instruments, residual gases released from materials can 
be analyzed using TGA-tandem instruments, such as TGA-FTIR or TGA-Mass Spectrometry, to 
determine the identity of the released gas and give insight into the weight loss mechanism. DMA 
measures the mechanical properties of polymer material as function of temperature and frequency, which 
reveals molecular relaxations in polymers. DSC is used to measure temperatures and heat flow during 
thermal transitions (glass transition, crystallization and melting) in polymeric materials. The degree of 
crystallinity of semi-crystalline polymers can also be obtained from the crystallization exotherm. Those 
methods have been applied to investigate thermal properties of the three polymers, which will be 
presented in Chapters II, III and IV. 

X-ray diffraction (XRD) and Infrared (IR) spectroscopy are also utilized. Both of these analysis 
methods are widely used to determine properties of polymers. XRD turns out to be a convenient and 
reliable method to investigate crystalline structure. The degree of crystallinity of polymers, which plays 
an important part in determining their dielectric properties, has been measured by XRD. IR spectroscopy 
is one of the most common spectroscopic methods applied to analyze organic compositions. It utilizes a 
Michelson interferometer and is based on IR absorption by dipolar molecules as they undergo vibrational 
and rotational transitions. Each peak in an IR spectrum indicates characteristic absorption regions for 
some commonly observed bond strength and bending deformations. It has been used to detect signs of 
oxidation due to thermal exposure of PI and ETFE. 

4. NDE of wire insulation 

The theme of this research is focused on evaluating wiring insulation status through capacitive methods. 
Insulation status can be characterized by its dielectric properties. Model-based capacitive methods 
developed in this research relate quantitatively the measurable capacitance to the dielectric properties of 
wires under test, and therefore allow for effective determination of wire insulation status. Experimental 
studies on realistic aircraft wires showed that dielectric property changes in wiring insulation due to 
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thermal and hydrolytic exposures can be successfully detected using the capacitive methods developed in 
this research, for wire type MIL-W-81381/12. 

Motivation 

This work is motivated by the effective evaluation of degradation status of air- and space-craft wiring 
insulation. Degradation in electrical wiring insulation has the potential to cause aviation catastrophe due 
to consequent short-circuiting or loss of control function [88]. Different wire inspection techniques have 
been developed over the past decade, for the purpose of replacing the traditional visual inspection method. 

Causes of failure and aging in aircraft wiring 

In [89], causes and modes of failure in legacy aircraft wiring have been categorized. These causes 
include chemical degradation such as corrosion of current carriers and hydrolytic scission of polymer 
chains, electrical degradation of wiring insulation that may be due to concentrated electric fields at sites 
of electrical stress and different kinds of arcing, and mechanical degradation due to vibration, over 
bending and other kinds of mechanical stress. 

Inspection techniques 

Visual inspection is probably the most widely used method for aircraft wire inspection. It is highly 
laborious while giving little quantitative information about the condition of the wires. Different physics- 
based wire inspection techniques have been developed over the past decade to replace this traditional 
inspection method, of which a summary is given here. 

Methods that can be applied for wiring conductor inspection 
These wire inspections methods can qualitatively determine if the wiring is faulty but are not suitable 

for inspection of aging aircraft wiring. Resistance measurement methods differentiate broken wires from 
good ones by measuring the end-to-end cable resistance. High resistance indicates broken wires (open 
circuit) while low resistance means the wiring is healthy (short circuit). The low-voltage resistance tests 
and dielectric-withstand-voltage tests can detect faults but are not suitable for miniaturization or 
pinpointing the fault [90], 

One of the most commonly used physics-based techniques for aircraft wiring testing is reflectometry, in 
which a high frequency electrical signal is sent down the wire and any impedance discontinuities in the 
testing wire results in reflected signals. The location of the fault can be determined from the time or phase 
delay between the incident and reflected signals whereas the impedance of the discontinuity is obtained 
from the magnitude of the reflection coefficient. An excellent review paper that compares different 
reflectometry methods is [91]. Reflectometry, however, is not capable of inspecting the insulation 
conditions. Reflectometry methods are distinguished by the types of incident voltage used. Time domain 
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reflectometry (TDR) uses a short rise time voltage step as the incident voltage. This method is susceptible 
to noises and is not optimal for live wire testing [92] [93] [94], Frequency domain reflectometry (FDR) 
uses a set of stepped-frequency sine waves as the incident voltage. A conceptual design of a "smart wiring 
system" based on FDR methods that can be used for on-board testing of aging aircraft wiring has been 
described in [90]. Phase-detection frequency-domain reflectometer (PD-FDR) has also been applied for 
locating open and short circuits in a Navy F-18 flight control harness [95], Sequence time domain 
reflectometry (STDR) and spread spectrum time domain reflectometry (SSTDR) use pseudo noise 
sequence and sine wave modulated pseudo noise code as the incident voltage, respectively. Testing 
systems based on these two techniques are capable of testing live wires and therefore have the potential to 
be used on energized aircraft to locate intermittent faults. The parameters that control the accuracy, 
latency, and signal to noise ratio for SSTDR in locating defects on live cables has been examined in [96], 
and the feasibility of spread spectrum sensors for locating arcs on realistic aircraft cables and live wire 
networks has been demonstrated in [97] and [98], 

In [99], linear relationships between the capacitance/inductance of open-/short- circuited wires (parallel 
insulated round wires, twisted-pair wires, and coaxial cables) and their length have been demonstrated 
and enables the determination of cable length from measured capacitance/inductance values. 

Methods that can be applied for wiring insulation inspection 

Infrared thermography systems and pulsed X-radiography systems have been developed as 

nondestructive testing methods of aircraft wiring [100]. Infrared thermography has the benefits of rapidly 
examining large areas of wiring and can serve as a global testing method, whereas a portable pulsed X-ray 
system can be used to obtain a radiographic image of the portion of the wire or cable. 

Ultrasonic methods have also been developed to obtain quantitative information about aircraft wire 
insulation [101]. These methods, by modeling insulated wires as cylindrical waveguides, have been able 
to relate extensional wave phase velocity to heat damage or exposure in wire insulation and thus provide 
quantitative information about the insulation condition. 

Acoustic and impedance testing methods aiming at locating intermittent faults in aircraft wires and the 
widely used Mil-Std-1553 data bus system have been reported in [102]. Micro-fabricated current sensors 
that could be located in key areas of the electrical wiring and interconnects systems have been reported in 
[103], Partial discharge (PD) analysis methods for diagnosing aircraft wiring faults are explored in [104], 
where a simulation of PD signal based on high-voltage insulation testing standard [105] has been detailed, 
followed by wavelet based analysis to de-noise the PD signals. 
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Capacitance methods developed in this research 

Deficiencies of the above methods suitable for wiring insulation inspection include the need of complex 
instruments in the measurement and not being able to provide quantitative information about the 
insulation condition at specific locations. A favorable solution to these deficiencies is capacitive methods, 
from which quantitative information about the permittivity of wiring insulation at specific locations can 
be obtained using not so complicated equipment. 

A curved patch capacitive sensor, with electrodes that conform to cylindrical test-piece surfaces, has 
been developed for wiring insulation evaluation. Numerical models have been developed and verified for 
both the homogeneous dielectric cylinder structure and the cylindrical structure of dielectric-coated 
conductors. Experimental studies on realistic aircraft wires showed that dielectric property changes in 
wiring insulation due to thermal and hydrolytic exposures can be successfully detected using the curved 
patch capacitive sensors, for wire type M1L-W-81381/12. 
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Chapter II. Polyimide 

Polyimide (PI) is widely used as an insulation material for machines and wiring, and is effective at 
temperatures up to 400 °C. Given the fact that polyimide may be exposed to extreme temperatures during 
unusual events in service, its thermal degradation kinetics and the effect of thermal degradation on its 
permittivity and electrical breakdown behavior have been studied. The lifetime of polyimide under 
electrothermal multi-stress is predicted by using a short term technique. As polyimide is commonly 
immersed in salt water while serving in navy aircraft, effect of water/saline exposure on its permittivity 
and electrical breakdown behavior is also investigated in this chapter. 

1 . Introduction 

Kapton UN is a polyimide film developed by DuPont which has been successfully used as electrical 
insulation in a wide range of temperatures, from -269 °C to +400 °C (4 K - 673 K) [3], The chemical 
name for Kapton HN is poly (4,4'-oxydiphenylene-pyromellitimide), and its chemical structure is shown 
in Figure 1. Kapton® Polyimide is produced from the condensation of pyromellitic dianhydride and 4,4'- 
oxydiphenylamine. In addition to its very light weight and advanced mechanical properties compared to 
other insulator types, Kapton HN polyimide has good dielectric properties, such as high breakdown field, 
low dielectric constant and low loss factor. Selected properties of 125 pm thick Kapton HN film are listed 
in Table 2 [3]. However, polyimide is very susceptible to hydrolytic degradation, which can give rise to 
cracks in the insulation and cause electrical malfunctions [4] . 



Melcher et al [5] explored the effect of moisture on the complex permittivity of polyimide film in a 
temperature range from 80 to 325 K. It is presented that the imaginary part e" measured for a film which 
was dried for two days shows only one maximum in the temperature range (the lowest curve). Absoiption 
of water alters this behavior for different water contents at fixed frequency 10 kHz. The peak height 
increases with water content and an additional smaller loss peak appears at its lower temperature shoulder. 

The shape of the larger peak, which is notified as the high-temperature peak, is the same for all film. 
The influence of the high-temperature peak can be subtracted because its shape is essentially independent 
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of the water content. But the height of the peak increases with higher water content. The second peak, the 
low-temperature peak, is considered to be strongly overlapped by the high-temperature peak. According 
to a statement in [6], since the high-temperature peak is present even at low humidity levels, it is proposed 
to be associated with water absorbed at the carbonyl groups. And the low-temperature peak is only visible 
at higher humidity, it is likely caused by water absorbed at the ether linkage. As the two loss peaks can be 
removed by drying the film, it is concluded that the water dipole causes this relaxation process and not an 
intrinsic dipole of the polyimide chain. 

It is also presented in reference [5] that the increase of the real part s' near room temperature correlates 
to the peaks in s". 

Thermal exposure of polyimide has also been explored [7] [8] [9] [10]. No significant changes in the 
dielectric properties of polyimide were observed after thermal exposure in air or N 2 from 200 to 350 °C 
for up to 5000 hours. 


Table 2 Selected properties of 125 pm thick Kapton®HN film. 


Property 

Value 

Dielectric strength 

154 kV/mm at 60 Hz, 23 °C and 50% RH 

Dielectric constant 

3.5 at 1kHz, 23 °C and 50% RH 

Dissipation factor 

0.0026 at 1kHz, 23 °C and 50% RH 

Melting point 

none 

Glass transition 

360 to 4 1 0 °C 

p-transition 

60 to 127 °C 

y-transition 

-118 to -28 °C 

Ultimate tensile strength 

231 MPa at 23 °C 

Impact strength 

78 N cm at 23 °C 

Yield point at 3% 

69 MPa at 23 °C 


2. Sample material 

All the PI samples under investigation in this chapter were cut from large sheets of 125-pm-thick 
Kapton® HN PI film obtained from Dupont. 

The storage modulus E', loss modulus E" and tan 8 of PI film were measured from -150 to 180 °C at 1 
Hz by DMA and analyzed by the software ‘TA Universal Analysis’ which can estimate values of peaks 
and shoulders in the curve to a tenth of a degree, as shown in Figure 2. The tan 8 curve reveals two 
molecular relaxations at approximately 60°C and 350 °C. The relaxation at 350 °C is attributed to the 
glass transition that occurs in the amorphous phase [3], while the relaxation at approximately 60 °C is 
associated with the P-transition, which is a sub-T g relaxation that takes place at temperatures between 60 
and 127 °C [11]. It is considered to be a result of torsional oscillations of the phenylene ring, involving 
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imide groups of PI [12]. In addition, PI exhibits y-transition in temperature range between -118 and 
-28 °C [11] due to increase in the vibration of aromatic groups as intra- and intermolecular interactions 
decrease in the presence of absorbed moisture [12]. 
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Figure 2 Results of dynamic mechanical analysis on dried PI. 


Since the permittivity and electrical breakdown strength of Kapton PI film changes significantly with 
moisture content [13], dry samples are needed in order to obtain baseline (control) values of dielectric 
strength for comparing with those obtained following thermal exposure and immersion in water. In order 
to determine a heating temperature to effectively remove water from PI, weight loss of a PI sample was 
monitored by thermogravimetric analysis (TGA) while it was heated from 30 to 900 °C in air at 
30 °C/min. As shown in Figure 3, an initial weight loss of 1% was observed at approximately 200 °C, 
which is attributed to loss of water from the sample during the heating process. Therefore, PI samples for 
baseline breakdown measurement were dried by heating at 200 °C for 1 hr, which, moreover, cannot give 
rise to degradation of PI. 

One dry PI control sample was coated with gold paint immediately upon removal from the isothermal 
furnace. Its complex permittivity was measured at frequencies from 1 Hz to 1 MHz over temperatures 
increasing from -140 °C to 180 °C in increments of 10 °C, by using a Novocontrol Dielectric 
Spectrometer with temperature-controlled standard sample cell. The real permittivity, s ' , and loss factor, 
tan 8, of the dry PI control sample over the whole temperature and frequency ranges are plotted as surface 
plots in Figures 4a and 4b, respectively. Below approximately 0 °C, e' increases with temperature; while 
it decreases with increasing temperature at higher temperatures. Given that polyimide is a polar polymer, 
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the dependence of its real permittivity on temperature is determined by how much its intra- and 
intermolecular interactions change with temperature [13]. If intra- and intermolecular interactions is 
independent of temperature, the permittivity would decrease with temperature. However, if intra- and 
intermolecular interactions change significantly with temperature, the dependence of permittivity on 
temperature would be governed by the change in intra- and intermolecular interactions and s' would 
increases with temperature. Therefore, it is speculated here that the intra- and intermolecular interactions 
of polyimide changes significantly with temperature below 0 °C, while at higher temperatures their 
change with temperature is reduced. It can also be seen that 8 ' decreases with frequency at each 
temperature; while tan 8 as a function of temperature reveals two peaks. The higher-temperature peak 
corresponds to the [S-rclaxation of PI. On the other hand, the lower-temperature peak can be ascribed to 
the y-relaxation of PI. The highest values of tan 8 observed at low frequencies and high temperatures are 
considered as a result of electrode polarization [14]. 



Figure 3 Percentage weight loss of PI as a function of temperature measured at 30 °C/min heating rate 
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Figure 4 The real permittivity (a) and loss factor (b) of dry PI over frequency range 1 Hz to 1 MHz and temperature range —140 

to 180 °C. 


3. Thermal degradation 


Thermal degradation kinetics 


Method and experiment 

The most common tool for analyzing polymer degradation is thermogravimetry (TG). TG measures the 
degree of degradation (as measured by mass loss) with respect to time (t) and temperature ( T ) [15]. The 
degree of degradation (a) for the case of total decomposition with zero rest mass can be defined as: 


a = f(t,T ) 


wt% 

100 


( 1 ) 


where wt% is the relative mass obtained directly from the TG experiment. 

TG experiments capturing the polymer degradation at different heating rates provide data that can be 
used to obtain degradation kinetic parameters, such as activation energy, for various reaction models. In 
this work, Kapton is analyzed by TG in an air environment to investigate the degradation in oxidative 
environments. Through the use of isoconversional kinetics, the advanced model mechanisms are 
identified. A mathematical model representing degradation is developed with an excellent statistical fit to 
the experimental TG data and is used to compare isothermal data. Finally, Fourier Transformed Infrared 
(FTIR) analysis and Mass Spectroscopy (MS) analysis of the exit gases identifies the breakdown 
components of Kapton to verify the complex degradation of Kapton. 
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A thermogravimetric (TG) analyzer, model Q50 from TA Instruments (New Castle, DE), was used for 
all of the TG experiments. The experiments were conducted from room temperature to 900 °C at five 
separate ramp rates: 2, 5, 10, 20, and 30 Kmin Under the controlled environment of the TG instrument, 
the samples were degraded in an air atmosphere using a balanced purge gas flow rate of 40 mL/min and a 
sample purge gas flow rate of 60 mL/min. Samples were placed on a platinum pan during the degradation 
process. Kinetic analysis was performed with the Netzsch Thermokinetics 2 program (version 2004.05) 
and standard statistical and plotting programs. Further study was conducted through evolved gas analysis, 
a technique utilizing MS and FT1R on exit gases from the TG experiments, to verify the degradation 
breakdown components and paths. Each test sample was punched out of the film using a circular punch, 
5 mm in diameter, ensuring reproducible sample weight and shape. The sample masses were 3.6 ± 0.5 mg. 


Kinetic modeling 

In degradation kinetics, the degree of degradation (Eqn. 1) varies from 0 (no mass loss) to 1 (complete 
mass loss). When modeling, two separate functions are assumed; K(T) and /(a), such that the governing 
differential equation has the following form: 

Y t = K(X)f{ a) (2) 

where da/dt is the rate of degradation, K(T) is the temperature-dependent rate constant, and /(a) 
corresponds to the reaction model [16]. The temperature-dependent rate constant is commonly described 
by the Arrhenius equation: 

K(T)=Ae~ E/RT (3) 

where R is the universal gas constant, E is the activation energy, and A is a pre-exponential factor [17]. 

When heating at a constant rate, Eqn. 2 can be redefined to eliminate the time-dependence by dividing 
through by the heating rate: 


da 

dT 


= p {a)e -E,ET 


( 4 ) 


where /? = dT/dt is the heating rate. 

Through linear transformation, the kinetic parameters ( A and E ) can be obtained by the time- 
independent rate equation: 




E 

RT 


( 5 ) 
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Eqn. 5 follows the linear formy = a 0 + a 1 x (withx = 1/7) and optimal fit of the kinetic parameters is 
determined using linear regression. By calculating these parameters through linear regression at several 
different mass losses, the variation in the kinetic parameters as a function mass loss is determined. 

In one approach for kinetic degradation modeling, constant activation energy and pre-exponential 
factors are assumed [17]. The model-free isoconversional method allows for varying kinetic parameters 
by assuming both the activation energy and pre -exponential factor are a function of the degree of 
degradation [18]. Freidman’s method, a well -known technique, obtains the activation energy by plotting 
the logarithmic fonn of the rate equation for each heating rate: 

< 6 > 

where the subscripts a and i represent the value at a particular degree of degradation and the data from 
the given heating rate experiment, respectively [17]. The activation energy at each degree of degradation 
is calculated with linear regression from a plot of Zn[/?j (da/dT) a i \ versus 1 /T ai across all of the 
heating rates tested. The Friedman plot not only provides confirmation of the multi-step processes during 
the reaction but also provides insight into the type of reaction steps. The type of reaction can be 
determined by comparing the slope of the constant fractional mass loss trend line to the slope of the 
constant heating rate data at each peak. The peak slope specifically refers to the slope of the linear portion 
to the right side of each peak. Comparing the relative magnitude of each negative slope, three types of 
reactions are defined: normal, accelerated, and retarded. A normal reaction corresponds to slopes of equal 
magnitude in both the fractional mass loss trend line and the peak slope of the constant heating rate 
data — Figure 5 A. An accelerated reaction corresponds to a steeper peak slope in the constant heating rate 
data compared to the fractional mass loss trend line — Figure 5B. A retarded reaction corresponds to a 
steeper fractional mass loss trend line compared to the peak slope in the constant heating rate data — 
Figure 5C. Similar to the Friedman method, kinetic parameters can also be calculated by the Ozawa and 
Flynn-W all integral isoconversional method [19] [20], 

Expanding the kinetic analysis from a single-step reaction to a multistep reaction, the differential 
equations are separated based on each step of the reaction. The overall degree of degradation is 
constructed as follows: 


a = 1 - I; Wj-aj (7) 

where a is the total fractional mass loss, a / is fractional mass loss of a specific reaction step, Wj is the 
contribution of a specific reaction step into total mass loss, and j represents the given reaction step [21]. 
The sum of the contributions of all steps is equal to 1 : 
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( 8 ) 


Each fractional mass loss of a specific reaction step can be written as an individual differential equation 
modeling the degradation of the reaction step such as [22]: 


d ( a ;^ a j+ 1) 
dt 



( 9 ) 


The rate of reaction for a degradation from A B (step 1) is given by d( a 1 -* a 2 )/dt. The rate of 
reaction for the degradation from B C (step 2) is given by d( a 2 -* a 3 )/dt. The rate of reaction for the 
degradation from C D (step 3) is given by d(a 3 -» a 4 )/dt. In this format of differential equations the 
values a l5 a 2 , a 3 and a 4 are the formal concentrations of the formal substances A, B, C, and D. A is the 
educt, B is the product of the first step and educt for the second step, C is the product of the second step 
and the educt for the third step, and D is the product of the third step which is the final product of the 
whole process. Each value of a, changes from 0 to 1. The initial state corresponds to a x = 1 , a 2 =0, a 3 =0 
and a 4 =0, and final state D corresponds to a 1 =a 2 =a 3 =0, and a 4 =l. If the reaction steps are completely 
separated, then the intermediate state after the first step corresponds to a^O, a 2 =l and a 3 =a 4 =0 and the 
intermediate state right after the second step corresponds to a 1 =a 2 =0, a 3 =l and a 4 =0. The degradation 
continues to follow the analogy of chemical kinetics, where step 2 follows step 1, step 3 follows step 2, 
but may occur before complete conversion of A to B. 

Results and discussion 

The TG scans for five different heating rates began at room temperature and the data can be seen in 
Figure 6. Like most polyimides, Kapton is extremely stable at intermediate temperatures [23], The onset 
of degradation increases with increasing heating rate and involves a rapid and complete degradation. The 
derivative of the weight with respect to temperature provides better insight into the mechanism of 
degradation. For a specific heating rate, the number of peaks in the derivative thermograms (DTG) 
represents the minimum number of reaction steps involved. By varying the heating rates, the degradation 
steps can be separated and isolated. At higher heating rates, for Kapton, the reaction mechanisms can be 
separated for better kinetic model understanding. Figure 7 shows the DTG curves. The peaks of the DTG 
help to determine the reaction steps. Based on the DTG curves below, heating at 30 K-min-1, a minimum 
of three reaction steps, or three peaks, can be seen. 

The Friedman plot for Kapton can be found in Figure 8. A multi-step reaction is again evident from the 
curvature of the plot. For each heating rate, there are separate reaction peaks. This indicates the 
probability of a multiple step reaction. Model-free analysis shows a complex process with three peaks for 
curves 30 Kmin 1 and 20 Kmin 1 and only two peaks for 2, 5, and 10 Kmin' 1 . The fluctuation in the 
number of peaks indicates that the mechanism of the decomposition changes with heating rate. 
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Furthermore, the type of reaction can be determined by comparing the fractional mass loss trend lines 
discussed previously with Figure 8. The fractional mass loss trend lines are the solid linear curves in 
Figure 8, and are found from linear regression at specific values of a ranging from 0.2 to 0.8. In all cases 
for Kapton, the peak slope is steeper than the fractional mass loss trend line indicating an accelerated 
reaction, probably autocatalysis. For autocatalysis, the generic governing differential equation, presented 
in Eqn. 2, defines the reaction model, f(aj), such that: 

f( a j) = (1 - «;) n ( 1 + Kcat^j) (10) 

where n represents the reaction order and K cat represents the autocatalysis constant. 

The Friedman analysis is used to calculate the activation energy (E a ) and the pre-exponential factor (A a ) 
from the slope and the y-intercept of the fractional mass loss trend lines, respectively [18] [19] [20], The 
activation energy and pre-exponential factor are shown in Figure 9 and presents activation energies from 
20 kJ/mol to 190kJ/mol. The plot of the activation energy with respect to the amount of degradation 
again confirms the multistep reaction by presenting non-constant activation energy throughout the entire 
degradation process. The fluctuating activation energy indicates an overlap of multiple reactions. As the 
reaction begins, the activation energy is about 190 kJ/mol and then shifts to 40 kJ/mol for a fractional 
mass loss of about 0.35. The activation energy increases to 60 kJ/mol for a fractional mass loss of 0.45, 
and then decreases to 20 kJ/mol for mass loss 0.8, and finishes by trending upward in the final moments 
of decomposition. The error bars show that the activation energy for the beginning of the reaction can be 
well-defined. For the last steps at the fractional mass loss 0.7 the error bar of activation energy is much 
higher and the lower value can reach almost zero kJ/mol. The error bars are calculated using standard 
error from the linear regressions defined by the Friedman Analysis. 

A physical meaning for the mass loss dependent activation energy from the Friedman Analysis is 
difficult to identify with confidence because of the independence of overlapping degradation mechanisms. 
Rather, the Friedman Analysis is useful in identifying multistep reactions. Given the complexity of 
backbone structure in polyimide, the chemical structure can rearrange in tandem with the degradation 
through aroyl migration or hydrolysis of the imido group. Dine -Hart et al. have proposed possible 
degradation pathways in their studies of polyimide film [24] [25]. 

An integral isoconversional method called Ozawa-Flynn-Wall Analysis was also used to calculate the 
activation energy as a function of fractional mass loss [19] [20]. Similar to the differential method used 
in the Friedman Analysis, the activation energy can be extracted using isoconversional trend lines. The 
benefit of comparing these two methods for activation energy provides insight into the type of reaction 
step to best model the degradation. Since the integral method for calculating activation energy cannot 
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utilize separation of variables, degradation kinetics involving competitive reactions show variations 
between the activation energies between the Friedman and Ozawa-Flynn-Wall Analysis. In conjunction 
with DTG peaks, the experimental data suggests a minimum of three steps with a combination of 
competitive and consecutive steps. 



Figure 5 Friedman plot for a single step (A) normal reaction, (B) accelerated reaction, and (C) retarded reaction. 
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Temperature (°C) 


Figure 6 TG curves broaden as the rate increases from 2 to 30 Kmin" 1 . 



Temperature (°C) 


Figure 7 DTG curves for the data shown in Figure 6. 
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Figure 8 Friedman plot from the data shown in Figure 6. 


For the simulation, a model of three parts was used. The schematic representation of the mechanisms 
can be seen in Figure 10. The first part is the process from reactant A to reactant B, which proceeds along 
two different paths 1 and 2. The second part is the one elementary reaction from reactant B to reactant C, 
and the third part is the process from C to D, which also follows two parallel paths. Two different paths 
for the third part of the model are necessary because the experimental data, Figures 5 and 6 show that the 
decomposition mechanism for the last 60% mass loss depends on the heating rate. A multivariate version 
of the Borchardt and Daniels method was used to determine to optimal fit of the kinetic parameters by 
multiple linear regression [26]. The results of the model fit can be seen in Figure 11, with parameters 
given in Table 3. 

These parameters come from the combination of Eqn. 2 with the autocatalytic reaction model found in 
Eqn. 10 for each step of the reaction diagramed in Figure 10. The Arrhenius parameters, E and A, are 
related to the temperature sensitivity of the reaction [27]. The reaction order and autocatalytic constant 
provide additional description of the chemical and physical reactions. The autocatalytic constant describes 
the extent in which the degradation reaction itself acts as a catalyst for that reaction. In thermodynamics 
of gases and liquids, the reaction order is an integer of stoichiometric equivalence. Flowever, if the 
reaction takes place in the solid-solid, solid-liquid, and solid-gaseous phases, physical processes influence 
the reaction rate such as diffusion, phase-boundary reactions, or nucleation. Therefore, the direct 
evaluation of experimental data with u nkn own reaction order gives non-integer values. With respect to 
the first 3 steps, the effects are minimal and an approximate reaction order of n = 1 can be used without 
effecting the model drastically. Yet, in the final two steps, there are significant variations to the reaction 
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order that cannot be approximated away. Therefore, the physical processes influencing the degradation of 
the final competitive reaction steps differ from the stoichiometric coefficients and play a significant role 
in the degradation. 

From a statistical perspective, the model follows the data with an r 2 value of 0.99991. The large 
activation energies of the initial steps and very small values of activation energies for the last steps are in 
agreement with results of Friedman analysis, if the error bars are taken into account. 

The kinetic parameters obtained for the multi-step model were then used to develop an isothermal 
model that would represent degradation across a 4-hour isothermal exposure for an air atmosphere, also 
obtained by TG analysis. We modeled a temperature spread of 425 to 475°C in 25°C increments. Figure 
12 shows a comparison of the isothermal experimental data and the mathematical models. Isothermal tests 
intrinsically have significant uncertainty because of the variability to achieve the set temperature. While 
the isothermal data and the model prediction differ by up to — 16%, the model captures the general trend 
and magnitude of weight loss shown in the experimental data. This error is within the bounds of a typical 
isothermal test, but could be brought down with improvements to isothermal experimental environment 
such as tighter temperature control and gradients. However, the isothermal model was able to accurately 
capture the general trend of the mass loss. 

To further explore the degradation process, evolved gas analysis was conducted during the TG 
experiments. This coupled FTIR and MS analysis involved monitoring the exit gases as a function of time 
during the TG experiments at different heating ramp rates. The 3-dimensional FTIR data can be seen in 
Figure 13 and the 2-dimensional FTIR data constructed for analysis is shown in Figure 14. The MS data 
is shown in Figure 15. 

There are four peak groups of interest as seen in the four sections in the FTIR figure. The first major 
group of peaks revolves around the 3500 to 3800 cm 1 wavenumber range. These peaks indicate the bond 
stretching of functional groups with removable hydrogen (bonds such as NH 2 , NH, COOH, and OH) [28], 
The second major group of peaks occurs in the wavenumber range 2100 to 2400 cm" 1 . These peaks carry 
the largest intensity and indicate bonds such as CO and C0 2 [28], This is to be expected since carbon 
dioxide and carbon monoxide are dominant products of degradation. The third set of peaks appears at 
wavenumbers around 1000 to 2000 cm" 1 . Since these fall close to the “fingerprint region” it is difficult to 
claim specifically what these peaks indicate. However, this peak area signifies bonds such as N=0, N=C, 
N0 2 , and fragmented aromatic rings [28]. The final set of peaks are found at around 720 cm" 1 and 
indicates nitro and nitroso compounds (NO bonds) [28], These conclusions correspond well to typical 
breakdown mechanisms proposed for Kapton. Dine-Hart showed that under oxidative degradation, 
Kapton initially evolves C0 2 through hydrolysis of the imido group followed by the decarboxylation of 
the resulting acid group, evolves C0 2 through an aroyl migration, and CO through extrusion of the imido 
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group [24] [29]. Blumenfeld also shows the evolution of C0 2 , CO, and H 2 0 through the initial oxidative 
degradation of the aryl ether groups while leaving behind the diimide group [30], Further degradation at 
elevated temperatures cleaves the diimide residual into its nitrogen compounds [30]. 

The mass spectroscopy data corresponds closely to the FTIR data as well. The first major peak 
identified is mass number 44. This peak, which spreads across the whole TG degradation curve, 
corresponds to carbon dioxide, which was shown in the FTIR data to be a major exit gas contributor [28], 
Peaks 17 and 18 correspond to water and more generally to the removable hydrogen of the functional 
groups also found in the FTIR data above [28]. Peak 22 is the double ionized peak for carbon dioxide and 
therefore is overshadowed by the first carbon dioxide peak [28]. Peaks 30, 45, and 46 are indicative of 
the nitro and nitroso compounds correlating to the FTIR peak in the 720 cm 1 range [28], One peak not 
shown, mass number 28, is the mass number for carbon monoxide [28]. While the FTIR data indicates 
evolution of carbon monoxide, the mass number also corresponds to nitrogen in its diatomic form and 
consequently the nitrogen in the air environment overshadows the carbon dioxide emission. 

Both the MS and 3-dimensional FTIR data validate the multistep reaction by presenting multiple peaks 
in the exit gas analysis. The 3-dimensional FTIR data, shown in Figure 13, clearly shows three peaks 
along the time axis, at around 2400 cnf 1 wavenumber. In addition, the MS data indicates overlapping 
peaks, shown in Figure 15, as the reaction progresses providing support for the multistep reaction chosen 
in Figure 10. 
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Figure 9 Activation energy plot for air atmosphere from (a) Friedman Analysis and (b) Ozawa-Flynn-Wall Analysis. 
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Figure 10 Schematic representation of the multistep reaction. 
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Figure 11 Best fit model of the TG data for the four-step reaction models in Fig. 6, with parameters given in Table 3. The curves 

represent the model and the symbols represent the experimental data. 



Figure 12 TG isothermal curves of experimental data and model prediction. 
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Table 3 Parameters used in the kinetic model. 



Activation Energy 
(kJ/mol) 

Log (Aj/s' 1 ) 

Reaction Order 

Log (K cat ) 

Step 1 

189.62 

8.40 

0.77 

0.94 

Step 2 

189.48 

8.38 

0.75 

0.93 

Step 3 

168.73 

7.73 

0.89 

0.09 

Step 4 

0.10 

- 3.43 

3.17 

1.65 

Step 5 

9.79 

- 3.17 

0.40 

1.38 
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Figure 13 3-dimensional FTIR data for exit gases of a 30 K-min-1 ramp rate TG on degradation onset. 
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Wave Number (cm-1) 

Figure 14 FTIR data for exit gases of a 30 Kmin-1 ramp rate TG on degradation onset for four different spectral ranges. Peak 
intensities are only proportional within each range and should not be compared from range to range. 
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Figure 15 MS data for exit gases of a 30 Kmin-1 ramp rate TG. 


Effect of thermal degradation on permittivity 

It has been presented previously in the literature that dielectric properties of PI are not changed 
significantly by heating at temperatures up to 350 °C for as long as 5000 hours either in air or N 2 [7] [8] 
[9] [10]. However, given the fact that electrical wiring insulation still has to work under extreme 
temperatures, such as may occur in the vicinity of extreme events during service, it is also valuable to 
explore the effect of thermal degradation at higher temperatures in air on the dielectric permittivity of PI. 
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The Kapton PI film was sectioned into 3 cm by 10 cm rectangular samples, which were heated at 400, 
425, 450 and 475 °C for 1, 2, 3, 4 and 5 hours in an isothermal muffle furnace. In order to eliminate 
possible distortion of the samples during the heating process, the two shorter edges of each PI sample 
were stabilized by a pair of mirror-finish stainless steel plates. Consequently, the majority of each sample 
was exposed to oxygen and eligible for permittivity measurement. After heating, the PI samples were 
removed from the furnace and cooled in air. With higher temperatures or longer heating times, the 
samples became more brittle, darker and thicker, which can be associated with the formation of 
oxidized layers on the sample surfaces during the initial degradation of PI [24] [29]. Thermal degradation 
at temperatures higher than 475 °C or for exposure times longer than 5 hr was not conducted because the 
PI became so badly deformed that its permittivity could not be measured successfully. 

The thic kn ess of all treated samples was measured by using a micrometer with uncertainty of 1 pm after 
cooling. Directly following thic kn ess measurement, the samples were coated with silver paint, serving 
both as film electrodes for permittivity measurement and to prevent significant moisture exchange 
between the samples and air. The silver paint dries quickly at room temperature, but requires curing for 
16 hr at room temperature to achieve a volume resistivity as low as 5x1 O' 6 to 1.25 x 10' 5 Q m. When the 
silver paint was totally cured, each sample was cut into three smaller samples for permittivity 
measurement. The permittivity of all treated samples was measured from 1 kHz to 2 MHz at room 
temperature, using the Agilent E4980A LCR meter coupled with a 16451 test fixture. The three 
identically degraded samples under all degrading conditions revealed highly similar values of s' and tan 5, 
the standard deviation measured for both is less than 1%. 

Figure 16 shows the loss factor of PI degraded at 475 °C for 3 hr as a function of temperature and 
frequency. Compared with Figure 4b, the height of the p-relaxation peak is substantially increased, 
indicating tremendous increase in the intensity of the p-relaxation [31]. The P-relaxation cannot be 
observed at frequencies higher than 500 Hz after thermal degradation. To explore the effect of thermal 
degradation on the y-relaxation, the frequency dependence of tan 5 at room temperature for the degraded 
PI is plotted and compared with that of the dry PI in Figure 17. Even though the intensity of the y- 
relaxation is not obviously changed, the peak for y-relaxation moves to higher frequencies after the 
thermal degradation. 

As shown in Figure 18, the real permittivity of PI measured at 1 kHz increases with time and 
temperature of thermal degradation in air. The real permittivity of the PI sample heated at 475 °C for 3 
hours is 4.15, versus 3.77 for the sample heated at 450 °C for 5 hours and 3.24 for the control sample. PI 
samples degraded at 475 °C for 4 and 5 hours were highly brittle and wrinkled, which meant that it was 
not possible to measure their permittivity successfully. On the other hand, tan 5 of PI reveals little change 
after thermal degradation at 400, 425 or 450 °C and is not shown here. After thermal degradation at 
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475 °C for 3 hr, however, tan 8 of PI as function of frequency is changed significantly, as shown in Figure 
17. 



Figure 16 The loss factor of PI degraded at 475 °C for 3 hr, measured over frequency from 1 Hz to 1 MHz and temperature from 

-140 to 180 °C. 



Frequency (Hz) 


Figure 17 The loss factor of dry PI and PI degraded at 475 °C for 3 hr as a function of frequency at room temperature. 
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Thermal exposure time (hr) 


Figure 18 Effect of thermal degradation time and temperature on s' of PI at 1 kHz. Error bars indicate the standard deviation in 

measurements on three nominally-identical samples. 

The activation energy for the p- and y-relaxations can be calculated from the Arrhenius equation: 


E — k In 


m 



(ii) 


where E is the activation energy of molecular relaxation; k is the Boltzmann constant; 7) and T 2 are the 
absolute temperatures of transition at frequencies /) and f 2 , respectively. The equation can be rewritten as 
a linear function relating 1/T and In co: 

i r E 1 

ln f = — - - ( 12 ) 


where C is a constant. By plotting the dependence of In f on 1/T, the activation energy can be obtained 
from slope of the line. Based on the measured complex dielectric spectra as shown in Figures 4b and 16, 
the Arrhenius plots of both p- and y-relaxations of the dry and degraded PI (at 475 °C for 3 hr) are given 
in Figure 19. For the dry PI, the activation energy is determined to be 0.48 ± 0.02 eV for the y-relaxation 
and 0.85 ± 0.05 eV for the p-relaxation. The uncertainties arise from estimated uncertainty in the slopes 
obtained via linear regression. Those values of activation energy determined here agree, to within 
experimental uncertainty, with values presented in the literature: 0.47 eV for the y-relaxation [32], and 
0.87 eV for the p-relaxation [33] [11]. However, p-relaxation of the degraded PI is only observed below 
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500 Hz and reveals no variation as function of temperature. Therefore, the Arrhenius plot for p-relaxation 
is a vertical line as shown in Figure 19, and its activation energy could not be determined. On the other 
hand, the activation energy of the y-relaxation is determined to be 0.46 ± 0.02 eV for degraded PI, which 
is not influenced by thermal degradation at 475 °C for 3 hr, despite the shift of the y-relaxation to higher 
frequency. 

Dine-Hart et al have studied the process of oxidative thermal degradation of PI for a short period in air 
by thermal analysis methods, saponification testing and infrared spectroscopy [24] [29]. It was observed 
that the chemical structure of PI film was significantly changed by thermal degradation [24]. Figure 20 
demonstrates pyrolysis process of imide groups during oxidative thermal degradation of PI, which gives 
rise to C0 2 and CO, breaks chemical bonding in imide groups and forms free radicals [29]. Thus, we 
hypothesize that the significant increases in s' of PI following oxidative thermal degradation at 400 to 
475 °C could arise from the formation of free radicals, which enhance electrical polarization in PI. 
Moreover, due to scission in chemical bonding in imide groups, the torsional oscillations of the phenylene 
ring become much easier, increasing the intensity of the p-relaxation. As discussion in section 3.1, 
release of gas was detected while a Kapton PI sample was heated from 30 °C to 900 °C with heating rate 
30 °C/min. Extracting data from Figure 13, Figure 21 shows the FTIR spectra of the PI sample at 30, 400, 
450 and 480 °C. Compared with the spectrum at 30 °C, the FTIR spectra at 400, 450 and 480 °C reveal 
new absorbance peaks at 2170, 2320 and 2370 cm' 1 , corresponding to chemical bonds of CO and C0 2 
[24], The mass also confirmed the release of C0 2 and CO at temperatures above 400 °C; and 
consequently endorse the degradation path of PI shown in Figure 20. Thus, we can conclude that the 
significant increase in the real permittivity and intensity of the p-relaxation are results of pyrolysis of 
imide groups of PI during oxidative thermal degradation. 
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Figure 19 Arrhenius plot for (3- and y- relaxations of dry PI and PI degraded at 475 °C for 3 hr. 
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Figure 20 Pyrolysis process of imide groups of PI during heating [29], 



Figure 21 FTIR spectra of Kapton polyimide at 30, 400, 450 and 480 °C. 


31 


Effect of thermal degradation on electrical breakdown behavior 

The effect of heating at temperatures up to 400 °C on electrical breakdown behavior of PI has been 
studied previously. In reference [34], decreases in both AC and DC dielectric strength of Kapton PI as 
the testing temperature increased from 20 to 250 °C were presented. Diaham et al [35], who tested at 
temperatures up to 400 °C, also observed a decrease in DC dielectric strength of PI fdm with increasing 
temperature. In contrast, it has been stated in references [7] and [9] that the dielectric strength of PI was 
enhanced after heat treatment below 350 °C for up to 300 hours. 

However, considering that PI as electrical insulation often has to function at even higher temperatures, 
such as in spacecraft [36] or in the vicinity of extreme events during service [37], it is also valuable to 
study changes in its breakdown behavior due to thermal exposure at higher temperatures. In this paper, 
the influence of short-period thermal exposure at 450 to 480 °C in air on AC dielectric strength of PI is 
presented. Since none of the papers mentioned above have explored the reason why dielectric strength 
was changed by the heating processes, the relationship between thermal exposure mechanisms and 
electrical breakdown of PI is explored here. 


Weibull distribution 

Given the random distribution of electrical weak points in insulation material, a different dielectric 
strength can be measured on nominally identical samples [38]. As a result of this behavior, statistical 
distribution analysis has been adopted to analyze the breakdown behavior of PI. Several statistical 
distribution methods are available for the description of electrical breakdown behavior, including Weibull, 
Gumbel, log-normal and recently proposed empirical failure probability methods [39], among which the 
Weibull method is the most commonly used for solid insulation materials [38]. The cumulative 
distribution function (CDF) [40] for two-parameter Weibull distribution is shown in the following 
equation [38]: 


F(E;a,j3 ) = l-exp 


f-Y 

V a) 


(13) 


where F is the probability of failure at an electrical field up to E, E is the measured breakdown strength of 
insulation, a is the scale parameter, and [> is the shape parameter, which is a measure of the dispersion of 
E and decreases with increasing scatter. Both a and p are always positive and obtained through least- 
squares regression in this work. In the least-squares regression method [38], 7, and X, values are assigned 
by the following equations: 


Y i =HE i ) 


( 14 ) 
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X, = ln[-ln[l-F(/,iV)]} 


(15) 


Then, 7, is a linear function of X„ 

Y i =mX 1 +c (16) 

in which the slope, m, and intercept, c are given by the following equations: 


Z(x,-x)(y,-y) 


m = 




c-Y - mX . 


(17) 


Parameters f> and a are estimated from m and c, respectively: 

j3 — —— , and (X — exp(c) . 
m 


(18) 


The CDF function as described in equation (34) is supposed to create a reasonably straight line in a 
typical Weibull plot, where logarithmic dielectric strength, E, is plotted on the horizontal axis; and the 
failure probability, F, is plotted on the vertical axis. (For example, see Figure 22 later.) 

Two-parameter Weibull distribution is a special case of three-parameter Weibull distribution, as shown 
in the following equation, where the local parameter y is zero. 


F(E\a,(3) = 1-exp 


E_ 

\aj 


(19) 


Weibull firstly employed mean ra nk s to estimate the probability of failure in 1951 [41], but later he 
recognized Benard’s approximation for median rank, which assigns the CDF of failure F(i, N), 




/ -0.3 
W + 0.4 


( 20 ) 


where N is the total number of measured dielectric strength values and i is the assigned rank when the 
dielectric strengths are listed in order from smallest to largest. Benard’s approximation is considered to be 
of good accuracy for plotting the Weibull distribution and estimating a and P; its uncertainty is only 0.1% 
for N = 50 and 1% for N = 5 [38]. The goodness of fit can be represented by a simple correlation 
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coefficient, P, which is a number between -1 and 1. As dielectric strength always gives a linear fit with 
positive slope, the P in this study locates between 0 and 1. The closer the P to 1 , the better the fit is. 

p _ ll^CXj-XXyYi-Y) 

In this paper, the two-parameter Weibull distribution is applied to statistically analyze the influence of 
thermal exposure and immersion in water on the electrical breakdown behavior of PI. As stated in 
reference [38], 20 measured values of breakdown strength for each sample under test is considered 
sufficient to provide meaningful results. Here, breakdown tests were conducted at 20 different points on 
each sample. 

Experiment 

The samples were cut into 6 cm by 10 cm sheets, large enough to measure dielectric strength 20 times 
at different points. The mass of each sample was measured using a balance with uncertainty 0.01 mg. The 
thickness of the samples was measured using a micrometer with uncertainty 1 pm. The PI samples for 
thermal exposure were heated at 475 °C for 1, 2, 3 and 4 hours in an isothermal muffle furnace to 
investigate the effect of thermal exposure time; and at 450, 460, 465 and 480 °C for 4 hours to study the 
effect of temperature. Thermal exposure for longer times or at higher temperatures was not conducted as 
the PI was so brittle and wrinkled that its dielectric strength could not be measured successfully. Directly 
after thermal exposure, the samples were removed from the furnace and cooled to room temperature. It 
was observed that the PI samples became more brittle, darker and thicker with increasing thermal 
exposure times or temperatures. Directly after cooling, the mass of the samples was measured again and 
compared with the value obtained before thermal exposure. 

Breakdown voltage was measured using a Dielectric Rigidity 6135.054 instrument designed according 
to standards ASTM D-149 [42] and D-876. The Dielectric Rigidity instrument consists of an oil- 

insulated transformer, capable of supplying AC voltage up to 60 kV, and a shut-proof resistant test 
chamber, which houses a bath for insulating liquid and a mobile support for electrodes with diameter of 
6.34 mm. In this work, Envirotemp® FR3® fluid, a vegetable oil with antioxidant additive, cold flow 
additive and colorant, was applied as an insulating and cooling medium. During the measurement process, 
the voltage across the PI samples was increased from zero at a rate of 0.5 kV/s and stopped at a value 
where the sample failed, which was recorded as its breakdown voltage, E BD . Twenty measurements of 
dielectric strength were made on each sample. 
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Results and discussion 

The two-parameter Weibull cumulative distribution function of dry PI and the samples heated at 475 °C 
for 1 to 4 hours is plotted in Figure 22. The symbols represent experimentally measured dielectric 
strengths and the solid lines represent a linear fit to the 20 data points, for each exposure condition. Since 
the symbols in Figure 22 reveal “S” shape curves around their linear fit, it is hypothesized that three- 
parameter Weibull distribution could give a better fit [40], The correlated coefficients of two-parameter 
and three-parameter Weibull distributions for the experimentally measured dielectric strengths have been 
calculated and are compared in Table 4. It is observed that for each thermal exposure condition, the 
highest correlation coefficients of the three-parameter Weibull distribution are very similar or even lower 
than those of two-parameter Weibull distribution. For 0 hr and 3 hr, the local parameters are negative, 
which is physically unreasonable as dielectric strength is always a positive number. As stated in reference 
[16], the correlated coefficient of the three-parameter Weibull distribution should be much higher than 
that of the two-parameter Weibull distribution to be an equally good fit. Therefore, it is concluded that 
three-parameter Weibull distribution cannot give a better fit than two-parameter Weibull distribution, and 
statistical analysis of electrical breakdown behavior of all PI samples investigated in this study will be 
conducted by using two-parameter Weibull distribution. 

The CDF of two-parameter Weibull distribution of the samples heated at 450 to 480 °C for 4 hours is 
compared with that of dried PI in Figure 23. The dielectric strength values measured for dry PI (0 hr) are 
in the range from 130 to 160 kV/mm, agreeing well with values for PI films with equivalent thic kn ess 
presented in references [35] [43] [44]: 100 to 200 kV/mm. 

With longer heating time or higher temperature, the linear fits to the measured data shift towards lower 
dielectric strength and reveal decreasing slope, indicating greater dispersity in the measured dielectric 
strength values. The scale parameter (a) and the shape parameter (P), obtained through least-squares 
regression, are applied here to quantitatively analyze the effect of thermal exposure on breakdown 
behavior of PI. The dependence of a and p on heating time and temperature is plotted in Figures 24 and 
25, respectively. Both a and p decrease with increasing heating time or temperature. Statistically, the 
short-period thermal exposure at 450 to 480 °C significantly decreases the dielectric strength of the PI 
samples and gives rise to much more randomly dispersed weak points in PI film. 

It has been stated previously [45] [46] that low-level degradation in polymers due to thermal exposure, 
such as chemical aging and physical aging, increases probability of electrical breakdown of the material. 
The chemical aging of PI during thermal exposure, which gives rise to pyrolysis process of imide groups 
and release of C0 2 and CO gases, has been discussed in section 2. On the other hand, as listed in Tables 5 
and 6, weight loss in PI increases with heating time or temperature. It is 2.90% for the sample heated at 
450 °C for 4 hr, and 21.78% for the one heated at 480 °C for 4 hr. The large weight loss following 
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thermal exposure of PI further supports the occurance of chemical aging of PI during thermal exposure. 
To explore the possibility of physical aging during thermal exposure, enthalpic relaxation measurement 
on a PI sample heated at 480 °C for 4 hr was conducted by using DSC. However, no endothermic peak 
was observed on the DSC curve, revealing no sign of physical degradation in the degraded PI. Therefore, 
the increased susceptibility of PI film to electrical breakdown following the heat treatment is considered 
to be a result of consequential chemical aging of PI, which gives rise to free radicals, as shown in Figure 
20 . 



Figure 22 The cumulative distribution function of the measured dielectric strength of PI samples heated at 475 °C for up to 4 
hours. Symbols represent experimental data and lines are obtained by least-squares fitting to the data. 


Table 4 Correlated coefficient of 2-parameter and 3-parameter Weibull distribution of dielectric strength of PI heated at 475 °C 


Time of heating at 475 °C 

Correlation coefficient, P 

Two-parameter Weibull 

Three-parameter Weibull 

Ohr 

0.9876 

0.9876, y is negative 

1 hr 

0.9903 

0.9902 

2 hr 

0.9938 

0.9933 

3 hr 

0.9451 

0.9602, but y is negative 

4 hr 

0.9770 

0.9862 


36 



Figure 23 As for Figure 41 but for PI samples heated for 4 hours at various temperatures from 450 to 480 °C. 



Heating time (hr) 


Figure 24 The Weibull-statistical scale parameter (a) and the shape parameter ((3) as functions of heating time. 
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Heating temperature (°C) 

Figure 25 The Weibull-statistical scale parameter (a) and the shape parameter ((3) as functions of heating temperature for 4 hr 

heating time. 

Table 5 Weight loss of PI samples heated at 475 °C for up to 4 hours. 


Heating time (hr) 

Weight loss (%) 

1 

1.39 

2 

3.04 

3 

7.08 

4 

14.51 


Table 6 Weight loss of PI samples heated for 4 hours at various temperatures from 450 to 480 °C. 


Heating temperature (°C) 

Weight loss (%) 

450 

2.90 

460 

4.18 

470 

8.21 

480 

21.78 


Lifetime prediction under electrothermal multi-stress 

Theory > and method 

The concept of using thermogravimetric analysis in conjunction with electrical stress testing to predict 
lifetime of insulation material was first proposed in 1971 by Toop [47]. The lifetime model relates 
electrical life theory to kinetic decomposition theory through TG analysis. Toop’s model is not applicable 
for all wire insulation materials, in particular these with extremely complex decomposition mechanisms. 
However, Toop has shown that the lifetime of a material such as Formvar, wire enamel with a simple 
thermogram, can be described as accurately from his proposed electrothermal model as it can by 
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following the failure criteria and test procedures given in ASTM D2307. This ASTM standard is 
developed for testing relative thermal endurance of film-insulated round magnet wire. 

Toop’s work is built on the chemical rate theory devised by Dakin, which provides a relationship 
between the material’s degree of degradation and physical properties [48], More recently, Montanari and 
Simoni have further explored the aging, degradation, and modeling with regards to electrical insulating 
materials including models such as tree-growth [49] , multi and space charge [49] which again build on 
the chemical rate theory of Dakin [50] [51]. 

In the calculation of electrothermal lifetime, the degree of degradation is linked to electrical breakdown 
voltage. By an empirical correlation, mathematical modeling of the kinetics of degradation can be 
combined with knowledge of the electrical breakdown voltage to obtain a single lifetime prediction. 

Now replacing a in equation (2) by electrical breakdown voltage V means that f(V) is the reaction 
model dependent on the breakdown voltage mechanism of the material. The function f(V) embodies the 
chemical and physical effects of applied voltage V on the test piece. The rate of change in voltage 
breakdown is constrained by the general differential equation for a changing material property, 

f = Kcnnv) (22) 

such that dV/dt is the change in electrical voltage breakdown and f(V) is the reaction model dependent 
on the voltage breakdown mechanism of the material. 

Rearranging and integrating the governing equation with r replaced by V gives: 

Qm =K(n X' dt (23) 

where limits V 0 and Vf are the breakdown voltages corresponding to the pristine and the failed samples 
respectively, and t f is the time for the breakdown voltage to reach the failure criterion. Taking equation 
(23) and applying these limits, another function of the breakdown voltage can be realized: 

[l^ = \amf v i = K(T)t r (24) 

and substituting the Arrhenius function for K(T~) gives: 

g(V f )-g(V 0 ) = At f e- E / RT (25) 

from which: 

t f = ^r[9(yf)-9(y 0 )\. ( 26 ) 

Similar to the development given in equations (22) through (26) for electrical breakdown voltage, the 
generalized differential equation, equation (2), is also a starting point from which to develop g(a), a 
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function of degradation analogous to g(V ). By taking the rearrangement in (4), and defining the 
parameter x = E/RT, the separation of variables procedure isolates the temperature and mass loss in the 
function to give: 

ftf(a) = !e~ x x~ 2 dx + C. (27) 

The initial and final conditions are T 0 , x 0 , and a 0 and T', x', and a' respectively. Ultimately, at failure, 
T 1 = 77. Therefore, equation (27) can be written as: 

~ A [g(a')~ g(a 0 )\ = ^ [p(x') - p(x 0 )] (28) 

where: 


p(x) = — f'- — dx. 

J x J x x 


(29) 


As seen in [47], when taking the logarithm of (26) the y-intercept is equivalent to the y-intercept of the 
logarithm of (28). Consequently, these equations can be assimilated as follows: 


+ < 30 > 

With the correlation found through p(x), the temperature, T , in equation (30) is defined as the failure 
temperature of the thermogram, Tr [47]. Since p(x) correlates back to the mass loss in equation (4), the 
failure temperature is defined as the temperature at which a corresponding mass loss intersects the slowest 
heating thermogram. See Figure 6 for 7V. 

A thermogravimetric analyzer, model Q50 from TA Instruments (New Castle, DE), was used for the 
TG experiments. The TG experiments were conducted in an oxidative environment with a nitrogen 
balance purge gas flow rate of 40 mL/min and an air sample purge gas flow rate of 60 mL/min. The tests 
were performed over a temperature range of 35 to 900°C at five different heating ramp rates: 2, 5, 10, 20, 
and 30 Kmin '. Kinetic analysis was performed with the Netzsch Thermokinetics 2 program (version 
2004.05). 

Samples for electrical breakdown testing were isothermally heated at 425, 450, 460, 465, 470, 475, and 
480 °C for up to 5 hr in an isothermal muffle furnace. These temperatures were chosen to provide uniform 
sample degradation because Kapton is thermally stable below 400 °C. Prior to degradation, the samples 
were all dried in a standard convection oven for 1 hour at 150 °C. A simple isothermal TG analysis of 
Kapton HN at 150 °C for 1 hour shows 2-3% mass loss and then stabilizes with no additional mass loss; 
indicating the removal of water from the sample. The samples were then placed in a sealed, pre-weighed, 
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plastic bag. The dried Kapton was weighed to obtain initial weight and then placed in the muffle furnace 
for degradation. On completion of isothermal heating the samples were again placed in the pre-weighed 
bag and weighed again to obtain their mass loss. After that, the breakdown voltage of the samples was 
measured. 

A 60 kV dielectric rigidity device (Ceast part number 6135.053 /6135.054) was used for breakdown 
voltage testing. The experiments were conducted at room temperature in a flame-resistant bio-oil liquid 
(Envirotemp® FR3TM Fluid) insulation bath. The tests were conducted in accordance with ASTM 
Standard D149-09 “Standard Test Method for Dielectric Breakdown Voltage and Dielectric Strength of 
Solid Electrical Insulating Materials at Commercial Power Frequencies.” Breakdown voltage was 
measured at twenty separate locations, evenly distributed across each 60 x 90 mnr specimen, and the 
values were averaged for each specimen. 

Experiment, results and discussion 

The kinetic pathway for the thermal degradation of Kapton was shown in our earlier work to fit a five 
step model with both competitive and consecutive reaction steps [21]. The activation energy shown in this 
work is the overall activation energy for the complete degradation of Kapton. When the activation 
energies for the individual steps are separated, five different constant activation energies are obtained. 
The initial reaction steps, a competitive model, contain activation energy of 189kJ/mol [21]. Since the 
degradation for the electrical breakdown analysis is from the initial degradation steps, this activation 
energy is used for calculating the lifetime prediction for Kapton. 

Measured breakdown voltage of degraded Kapton is plotted in Figure 26. A linear fit of the breakdown 
voltage as a function of degradation is used to represent the trend of the data. This trend line is useful in 
calculating the lifetime to electrical failure. Each data point in the figure represents the average of twenty 
test measurements on a degraded sample sheet. The error bars represent the standard error with a 95% 
confidence interval between the twenty test values. 
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Figure 26 Measured voltage breakdown of degraded Kapton Film (symbols) with best linear fit (solid line). 

In the electrothermal lifetime prediction model of (30), time to failure, tf, is calculated from test 
temperature and voltage, TG analysis parameters E and Tf, and the breakdown voltage value. First, a 
breakdown voltage value is selected (e.g. 12kV) and the relevant fractional mass loss is read from 
Figure 5 (20%). 

This fractional mass loss defines failure temperature, Tf , from Figure 6 — the temperature of the slowest 
thermogram corresponding to the relevant fractional mass loss (541 °C or 814 K). The failure 
temperature along with the activation energy of 189 kJ/mol are then utilized in (29) to calculate the 
correlation fitting parameter p(x). 

Using equation (30), the lifelines shown in Figure 27 are calculated based on two particular sets of 
environmental conditions; voltages of 12 and 14.7 kV (10 and 20% mass loss respectively) and a 
temperature range from 250 to 400 °C. Lifelines for other sets of conditions can also be obtained but these 
two are chosen for demonstration here. The values of Tf, a , Vf, and E used are found in Table 7. 
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Figure 27 Predicted time to failure at 12 and 14.7 kV and for isothermal temperatures ranging from 250 to 400 °C. 
Table 7 Parameters values used to obtain the calculated lifelines shown in Figure 27. 


v f 

(kV) 

a 

(%) 

(K) 

E 

(kJ/mol) 

P 

K/hr 

t f at 350°C (hr) 

12 

20 

814 

189 

120 

1200 

14.7 

10 

801 

189 

120 

729 


As stated previously [21], the model presented in this work is not applicable to all materials; 
specifically those with extremely complex degradation mechanisms. The Kapton film studied here may, 
however, be treated in this way by considering only the initial part of the weight loss curves (within the 
first 25%) where a simple mechanism is responsible for degradation. Further, the thermogram for Kapton 
is no more complex than the thermogram for Formvar, treated in [21], in which a polyimide film was also 
treated. 

Within the first 25% of weight loss, Kapton remains functional as wire insulation. Beyond 25% weight 
loss, the thin film failed to maintain the mechanical properties required for standard wire insulation and it 
was not possible to collect breakdown voltage data. At the elevated temperatures of this degradation study, 
the the surface of the Kapton film oxidizes. The oxidative layer can be seen on the surface of the material 
as darkening of the material color and alligator cracking of the surface. The oxidation causes the material 
to become brittle and deformed as the size of the oxidative layer expands — causing the material to fold 
over on itself. 

Beyond 25% weight loss, the breakdown voltage is expected to drop rapidly, showing a large decrease 
in breakdown voltage for minor weight losses. This expectation comes from the loss of material 
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properties due to the extensive oxidation of Kapton. In fact, beyond 25% weight loss, there are indications 
that the Kapton polymer backbone is degraded enough that most of the linear chains are broken. The 
change in activation energy seen in Figure 9 for the first 25% weight loss indicates the rupture of the 
Kapton backbone — the sharp drop in activation energy is related to the breaking of the stable covalent 
bonds. When the linear chain is ruptured repeatedly, commonly through aroyl migration or hydrolysis of 
the imido group, the free radicals, mobile molecules, and oxygen molecules provide pathways for 
continued degradation at lower activation energies as observed between 30 and 100% weight loss [21]. 
These types of new pathways have been observed in the mass spectroscopy data from TG experiments in 
which the peak indicating carbon dioxide emission clearly indicates two unique mechanisms. The first 
mechanism (backbone rupture) endures for the first 25% of degradation and the second unique 
mechanism (new pathway) appears at around 25% weight loss [21]. Additionally, the nitrogen 
compounds have distinguishable mechanisms that again indicate a transition at 25% weight loss from the 
rupture of the major chain to the breaking of subcomponents. 

To summarize, Kapton film may be treated by electrothermal analysis based on the criteria used in [47], 
provided that the lifetime prediction stays within the bounds of Kapton’s functional mechanical properties, 
i.e. within 25% weight loss. In this 25% weight loss, electrical breakdown characteristics and physical 
changes are related through free volume and free radicals. Beyond the 25% threshold, the mechanical 
properties are destroyed through the rupturing of the stable covalent bonds in the polymer backbone. This 
is seen through observing the physical appearance of the degraded Kapton, the reduction in activation 
energy around 25%, and through the evolved gas analysis reported in [21]. 


4. Water/saline exposure 

Effect of water/saline exposure on permittivity 

All PI samples for saline exposure experiments were cut into strips 3 cm wide by 9 cm long and 
immersed in distilled water or saline solutions for various times, including 0.5, 1, 2, 4, 6, 12 and 24 days. 
Saline solutions with three different concentrations were made by mixing distilled water with sodium 
chloride: 0.45 g/1, 6 g/1 and 80 g/1, with conductivities 1 mS/cm, 10 mS/cm and 100 mS/cm [52], 
respectively. Immediately after removal from water or saline solutions, the samples were wiped to remove 
water on their surfaces, their thickness was measured by using the micrometer, and the samples were 
coated with silver paint introduced in section 2. When the silver paint was totally cured, each sample was 
cut into three smaller samples for permittivity measurement. The permittivity of all treated samples was 
measured from 1 kHz to 2 MHz at room temperature, using the Agilent E4980A LCR meter coupled with 


44 


a 16451 test fixture. The three identically degraded samples under all degrading conditions revealed 
highly similar values of s' and tan 5, the standard deviation measured for both is less than 1%. 

Figure 28 shows the real permittivity and loss factor measured at 1 kHz as a function of sample 
immersion time in water and saline solutions. In the case of both water and saline, the real permittivity 
and loss factor of PI exhibit rapid increases in the first 4 days, but neither of them changes obviously after 
6 days. Figures 29 and 30 show how s' and tan 5 of PI behave as a function of frequency from 1 kHz to 2 
MHz and immersion time, in water and 80 g/1 saline, over 4 days, s' and tan 8 of PI immersed in water are 
increased by approximately 11% and 40%, respectively; and by 12% and 40% for the sample immersed in 
80 g/1 saline. As shown in Figures 29b and 30b, the intensity of the y-relaxation increases with time of 
immersion, because the PI samples absorb more moisture with longer time of immersion. This 
observation agrees well with the findings by mechanical testing methods [32], The frequency of the y- 
relaxation, however, is not changed as a function of immersion time. 
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Figure 28 The real permittivity (a) and loss factor (b) of PI immersed in water and saline solutions, measured at 1 kHz. Error 
bars indicate the standard deviation in measurements on three nominally-identical samples. 
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Figure 29 The real permittivity (a) and loss factor (b) of PI following immersion in distilled water. 
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Figure 30 The real permittivity (a) and loss factor (b) of PI following immersion in 80 g/1 saline. 


In order to further explore the effect of dissolved sodium chloride on the complex permittivity of PI, the 
differences in s' and tan 8 between the samples immersed in water and those immersed in saline are 
plotted as function of immersion time in Figure 31. It can be concluded that the variation in salinity does 
not change complex permittivity of PI significantly and that the increases in s' and tan 5 are as a result of 
moisture absorption. 
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The reason for the increase in real permittivity and loss factor of PI by moisture absorption can be 
attributed to the formation of polar groups and chain scission during hydrolytic degradation. Deiasi and 
Russell [53] observed a decrease in number-average molecular weight of PI caused by chain scission due 
to a chemical reaction between water and imperfections in the polymer chain. The chain scission 
mechanism is shown schematically in Figure 32 [54], which agrees with the statement in reference [6] 
that water molecules can bond to the carbonyl groups of PI. It has also been stated that the chemical 
reaction is rapid in the early stages of aging, and then slows as the number of susceptible sites in the 
polymer chain decrease [55], This explains why the real permittivity and loss factor of PI are increased 
significantly during the first 4 days of immersion, after which there is little further change. 




Figure 31 Effect of dissolved sodium chloride on the real pemiittivity (a) and loss factor (b) of PI, measured at 1 kHz. 
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Figure 32 Chain scission mechanism of PI hydrolysis through interaction of H 2 0 with the carbonyl groups [54], 


Effect of water exposure on electrical breakdown behavior 

PI is susceptible to moisture absorption. It has been reported that, while a MIL-W-81381 wire was 
undergoing accelerated aging in a humidity chamber, cracks propagated from the PI insulation tape to the 
conductor, acting as a current leakage path, and finally caused destructive failure in the circuit [4]. 
Increases in both the dielectric constant and the loss factor of PI with increasing water uptake have been 
observed [5]. However, even though Dupont has reported linear decrease in AC dielectric strength of PI 
as a function of humidity [3], few have statistically studied the effect of moisture absorption on 
breakdown behavior of PI. This paper statistically analyzes how AC dielectric strength of PI changes as 
the material is immersed in water for up to 96 hours by using Weibull distribution analysis. 

As concluded in section 3.2 that the effect of dissolved sodium chloride on permittivity of PI immersed 
in salt water is minor, it is speculated that electrical breakdown behavior could not be influenced by 
variation in salinity significantly either. Thus, study of electrical breakdown behavior of PI after saline 
exposure is not investigated here. 

Experiment 

The samples were cut into 6 cm by 10 cm sheets, large enough to measure dielectric strength 20 times 
at different points. The mass of each sample was measured using a balance with uncertainty 0.01 mg. The 
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thickness of the samples was measured using a micrometer with uncertainty 1 pm. The PI samples were 
immersed in distilled water at room temperature for 4, 8, 16, 24, 48, 72 and 96 hours. The surfaces of the 
samples were wiped dry after they were removed from the water. Directly after wiping, the mass of the 
samples was measured again and compared with the value obtained before immersion in water. Twenty 
measurements of dielectric strength were made on each sample by using the Dielectric Rigidity 6135.054 
instrument. 

Results and discussion 

The Weibull cumulative distribution function of the PI immersed in water for 0, 4, 8, 16 and 24 hours is 
shown in Figure 33. The straight line fit to each data set moves towards lower values of dielectric strength 
with increasing immersion time. Moreover, slopes of the lines decrease as immersion time increases. As 
described in section 4.2, a and (3 parameters are obtained via least-squares regression. Figure 34 shows 
how a and (3 decrease with longer PI immersion in water. 

As shown in Table 7 and Figure 35, no significant further change in either weight gain or dielectric 
strength of PI was observed after immersion in water for more than 24 hours. This observation is in line 
with the statement in reference [5] that PI is saturated with moisture after immersion in water for 25 
hours. 

Formation of side ionic groups due to dissociation of polymer chains is considered to be the reason why 
the insulation capability of PI film is deteriorated due to water uptake. Figure 32 shows the hydrolysis 
scheme for PI, in which polymer chains are broken and chemical composition is changed. OH- from 
water molecules bonds to either of the four carbonyl groups and H+ bonds to the amine groups, 
generating ionic groups at sides of the new polymer chain [4] [6] [53], The side ionic groups can greatly 
promote charge transport in PI insulation, which is termed as “intrinsic ionic conduction” in reference 
[45]. Consequently, electrical failure of the PI film is significantly accelerated. 

It has been stated that chemical changes in PI due to moisture absorption are fully reversible by 
annealing for 1 hour [56]. Thus, in order to investigate whether physical degradation occurred in PI 
during immersion in water as well as chemical degradation, a sample immersed in distilled water for 24 
hours was dried again by heating at 200 °C for 1 hour in order to eliminate the influence of ionic side 
group formation on dielectric strength of PI. The Weibull plot of breakdown for the annealed sample (not 
shown here) revealed no obvious change when compared with that of the control sample. As listed in 
Table 8, the scale parameter, a, of the sample is not significantly changed. Therefore, it is concluded that 
immersion in water did not give rise to significant physical degradation of the PI film. The change in the 
shape parameter, (3, could be a result of variation in dispersity of electrical weak points in the two samples. 
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Figure 33 The cumulative distribution function of the measured dielectric strength of PI samples immersed in water for 0, 4, 8, 

16 and 24 hours. 




Figure 34 The Weibull-statistical scale parameter (a) and the shape parameter ((3) as functions of time of PI immersion in 

distilled water. 
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Table 8 Weight gain of PI samples immersed in distilled water for up to 96 hours. 


Immersion time (hr) 

Weight gain (%) 

4 

0.87 

8 

1.32 

16 

1.75 

24 

1.80 

48 

1.80 

72 

1.81 

96 

1.80 


Table 9 The Weibull-statistical scale parameter (a) and the shape parameter ((3) of dry PI. 


sample 

a (kV/mm) 

P 

The original dry PI 

152.2 ±0.3 

32 ±2 

2 nd dry PI 

154.4 ±0.2 

40 ±2 



Figure 35 The cumulative distribution function of the measured dielectric strength of PI samples immersed in water for 24, 48, 

72 and 96 hours. 
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5. Summary 

This chapter investigates the effects thermal degradation in air and immersion in water/saline on the 
insulating performances of Kapton PI film. FT1R and MS exit gas analysis verified the degradation 
components and the the multistep degradation process of Kapton. After either the heating process or 
immersion in water, the real permittivity and loss factor of PI are increased substantially. Intensity of the 
p-relaxation is substantially increased due to scission of chemical bonding in imide groups during thermal 
degradation at 475 °C for 3 hr. As PI samples absorbed more moisture with longer time of immersion in 
water or saline, the intensity of y-relaxation following saline exposure is increased. On the other hand, 
using statistical analysis by two-parameter Weibull distribution, we observed that the dielectric strength 
of PI is greatly decreased and individual measurements of breakdown voltage are much more randomly 
distributed, with increasing time or temperature of thermal exposure and increasing time of immersion in 
water. The deterioration of insulating properties of PI after the degrading processes is explained in terms 
of chemical degradation due to pyrolysis of imide groups during thermal exposure and formation of ionic 
side groups during the hydrolysis process of PI while immersed in water. However, dissolved sodium 
chloride shows minor influence on dielectric properties of PI immersed in saline solutions. A short term 
technique based on thermal analysis and electrical breakdown measurement is verified to be a good 
solution to predict lifetime of PI under electrothermal multi-stress. 
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Chapter III. Polytetralluoroethylene (PTFE) 

Due to its excellent dielectric properties and hydrophobicity, PTFE (polytetrafluoroethylene) is widely 
used as an outer coating material for wiring insulation [1], Therefore, it is always exposed to heat and 
mechanical stress during service. In this chapter, the effect of thermal exposure and tensile strain on 
permittivity of PTFE is investigated. 

1 . Introduction 

PTFE, fCFj-h,, is a fluorocarbon polymer, typically with a very high molecular weight. The 
substitution of fluorine for hydrogen causes the material to exhibit extreme properties. Due to the C-F 
bonds, it exhibits special properties surpassing those of most polymers due to the substitution of fluorine 
for hydrogen, such as very high melting temperature and good chemical resistance. The main physical 
and chemical properties of PTFE are compared with those of polyethylene (PE) in Table 9 [57]. 

In PTFE, closed connected amorphous and crystalline phases coexist due to its extremely high 
molecular weight (~ 1 x 1 0 6 ) [58], The temperature-pressure phase behavior of crystalline PTFE is shown 
in Figure 36 [59, 60], At atmospheric pressure, the room temperature crystalline structure of PTFE 
(phase IV) transfers to phase II below 19 °C [61] and to phase I above 31 °C [62]. The first-order 
transition at 19 °C from phase II triclinic to phase IV hexagonal reflects an untwisting in the helical 
conformation from 13 atom/180 degree turn [63] [64] to 15 atoms/tum [63] [65] [66] and an increase in 
the hexagonal lattice spacing. Above 3 1 °C the individual polymer chains lose their well defined helical 
repeat unit [63] [67]. Further rotational disordering and untwisting of the helices produces a pseudo- 
hexagonal structure. At room temperature, phase II transfers to phase III orthorhombic above -0.65 GPa 
[68]. Amorphous PTFE has the same repeat atomic structure as the crystalline domains but without 
significant order [60], 

PTFE has excellent electrical insulating properties due to its low relative permittivity (2. 0-2. 2 [69]), 
low dielectric loss (<0.0002-0.0005 [70]), good frequency stability in a wide spectral range (up to 
10 GHz), and high breakdown strength (19.2 kV/mm [70]). Reference [71] shows the complex 
permittivity of PTFE as a function of temperature at microwave frequencies f ~ 11.5 GHz. The real 
permittivity of PTFE decreases from 2.12 to 2.01 as temperature increases from -151 to 102 °C and the 
dissipation factor is less than 10' 3 even at high temperatures. As presented in reference [72], the real 
permittivity of PTFE decreases from 2.35 to 2.09 as temperature increases from -150 to 250 °C at 10 kHz. 
Additionally, because of its good mechanical and thermal stability over temperatures up to 260 °C [73], 
PTFE has been extensively applied as wire insulation material in air- and spacecraft. 


P (GPa) 
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Table 10 Properties of PTFE and PE. 


Property 

PTFE 

Polyethelyne 

Density (g/cm3) 

2.2-2. 3 

0.92-1 

Melting Temperature (°C) 

342±10 (1st) 
327±10 (2nd) 

105-140 

Tg(a-PTFE) (°C) 

126 

- 

Tg (p-PTFE) (°C) 

30 

- 

Tg (y-PTFE) (°C) 

80 

- 

Surface Energy (dynes/g) 

18 

33 

Resistance to Solvents and 

Excellent, no known 

Susceptible to hot 

Chemicals 

solvent. 

hydrocarbons. 
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Figure 36 Temperature-pressure phase diagram of crystalline PTFE with the inter- and intra-polymer chain crystalline structures. 


2. Sample material 


In this work, all PTFE samples were cut from the same large sheet supplied by Fluoro-Plastics, which is 
nominally 610x610x1 mm 3 . The melting point of as-received PTFE was determined to be 327 °C by a 
DSC (differential scanning calorimetry) test. The PTFE sample appears opaque white and relatively hard 
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below its melting temperature, as shown in Figure 37a; while turned to be translucent and very easily 
pierced with a needle above the melting temperature, Figure 37b. 



Figure 37 PTFE (a) below the melting temperature 327 °C; and (b) above the melting temperature. 



Figure 38 Results of dynamic mechanical analysis on as-received PTFE. 


Figure 38 reveals the storage modulus E', loss modulus E" and tan 8 = E"/E' of an as-received PTFE 
sample as a function of temperature, which was measured by DMA (dynamic mechanical analysis) at 1 
FIz. The shape and intensity of peaks in the tan 8 curve, as well as the magnitude of the storage modulus, 
are determined by the degree of crystallinity in the PTFE. Thermal transitions at 36.1 °C (P transition) 
and 129.3 °C (primary a transition) can be observed according to peaks of tan 8 analyzed by the software 
‘TA Universal Analysis’ which can estimate values of peaks and shoulders in the curve to a tenth of a 
degree. These peaks correspond to 29.2 °C and 120.0 °C in the loss modulus curve. The above values are 
close to the data listed in Table 9 and are discussed more fully in reference [19]. Due to limitation in the 
temperature range of the DMA measurement, the y-transition at around -80 °C is not shown. 
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The real permittivity of as-received PTFE was measured over frequency range from 1 Hz to 1 MHz at 
room temperature by Novocontrol spectrometer, as shown in Figure 39. The error bar was given by the 
Novocontrol spectrometer, which includes the accuracy of the used impedance analyzer and the 
influences of cables and sample holders. The error bar also strongly depends on frequency of 
measurement and the actual sample impedance. As expected, the permittivity of PTFE is extremely stable 
over the whole frequency range. The mean value and standard deviation of the data points in Figure 39 
are 2.214 ± 0.001. Since the dissipation factor of PTFE is so small that it is beyond the measurement 
accuracy of this spectrometer over most of the frequency range, the measured values of loss factor are not 
shown here. 

To explore dielectric performance of PTFE at extreme temperatures, the permittivity was also measured 
over temperature range -150 to 300 °C at frequency from 1 Hz to 1 MHz. During measurements, the 
whole sample cell was housed in a temperature-controlled container, in which the temperature was 
gradually increased in 10 °C increments. Figure 40 shows the real permittivity of as-received PTFE as 
functions of both frequency and temperature. Drawn with data points from Figure 40, Figure 41a shows 
the real permittivity as a function of temperature at 1.15 kHz. At this frequency, s' decreases from 2.287 
to 1.968 as temperature increases from -150 °C to 300 °C. A probable explanation of this behavior is that, 
as a nonpolar polymer, the polarization of PTFE is mainly electronic. In that case, the Clausius-Mossotti 
equation can be applied, as stated in [13], 


where N is the number of polarizable entities per unit volume, a e is electronic polarizability, s 0 is the 
permittivity of free space, s is the relative permittivity of PTFE, M is molar weight and p is mass density. 
Rearranging the above equation gives 


Na e _ (£-1 )M 
3£ 0 0+2)p 


(43) 


£ = 


3 Me 0 


(44) 


M£ 0 -Na E p 


Since density decreases as the temperature increases, it can be postulated from the above equation that s' 
should decrease slightly with temperature, as the measurements show in Figure 41a). 
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Focusing on the temperature range -10 °C to 50 °C in Figure 41b, two step-like decreases in s' are 
observed near 19 °C and 31 °C, which can be associated with structural first-order phase transitions at 
these two temperatures [72]. The room-temperature crystalline structure of PTFE (phase IV hexagonal) 
transfers to phase II triclinic below 1 9 °C by angular displacement [6 1 ] and to phase I pseudo-hexagonal 
above 31 °C by crystal disorder [62], which are both accompanied by step-like changes of the specific 
volume. Due to the thermal expansion of the sample capacitor, s' of PTFE shows two step-like decreases 
in the vicinity of the phase transition temperatures. 



Frequency [Hz] 

Figure 39 Real permittivity of as-received PTFE as function of frequency at room temperature. 



Figure 40 Real permittivity of as-received PTFE as a function of frequency and temperature. 
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Figure 41 Real permittivity of as-received PTFE as a function of temperature at 1.15 kHz. a): -150 to 300 °C; b): -10 to 50 °C. 


3. Thermal degradation 
Thermal degradation kinetics 

Method and experiments 

A thermogravimetric analyzer, model Q50 from TA Instruments (New Castle, DE), was used for all of 
the TG experiments. Experiments involving PTFE were conducted from room temperature to 700 °C at 
five separate ramp rates: 2, 5, 10, 20, and 30Kmin''. Under the controlled environments of the TG 
instrument, the samples were degraded in an air atmosphere using a balanced purge gas flow rate of 40 
mL/min and a sample purge gas flow rate of 60 mL/min. Samples were placed on a platinum pan during 
the degradation process. Kinetic analysis was performed with the Netzsch Thermokinetics 2 program 
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(version 2004.05) and standard statistical and plotting programs. The same kinetic modeling was used as 
studying thermal degradation kinetics of polyimide, as described in Chapter II. 

Results and discussion 

The TG scans for PTFE began at room temperature, and degradation was not observed until elevated 
temperatures were reached. The high temperature TG data is given in Figure 42. Like most 
fluoropolymers, PTFE is extremely stable at intermediate temperatures. The onset of degradation 
increases with increasing heating rate and involves a rapid and complete degradation. From the TG curves, 
the variation between the temperatures of the degradation onset points provides insight into the reaction 
model. The uniform temperature spread of the onset point for the PTFE degradation curves indicates a 
single step for the degradation. 

The derivative of the weight with respect to temperature provides better insight into the mechanism of 
degradation. For a specific heating rate, the number of peaks in the derivative thermogram (DTG) 
represents the minimum number of reaction steps involved. By varying the heating rates, the degradation 
steps can be separated and isolated. Figure 43 shows the DTG curves of PTFE. Notably, only one peak is 
observed for the five heating rates. Based on the single DTG peak, a single step reaction can be seen. 

Once initial observations of the thermograms have been made, a Friedman analysis based on equation 
(6) provides information about multi-step processes during the reaction and also offers insight into the 
type of reaction step. Similar to the DTG plots, the number of peaks indicates the minimum number of 
reaction steps in the degradation. For PTFE there is one single peak as seen in Figure 44. Therefore a 
single step reaction is again evident from the curvature of the plot for PTFE. 

The type of the reaction can also be determined by comparing the slope of the constant fractional mass 
loss trend line to the slope of each peak at the beginning of the reaction step [21]. For PTFE, the first 
peaks have in all cases the acceleration, probably autocatalysis. For autocatalysis, the generic governing 
differential equation, presented in equation 2, defines the reaction model, /(a), such that: 

f{a) = (1 - a) n (l + K cat ) 0 (45) 

where a represents the degree of degradation, n represents the reaction order, K cat represents the 
autocatalysis constant, and X represents the step of the reaction. 

From the Friedman analysis as well as the Ozawa-Flynn-Wall analysis, activation energy as a 
function of fractional mass loss can be obtained [18] [19] [27]. For PTFE, the plot of the activation 
energy with respect to the amount of degradation again confirms the single step autocatalytic reaction by 
presenting linearly declining activation energy seen in Figure 45. Activation energy developed from 
Friedman analysis presents activation energies between 130 kJ/mol and 70 kJ/mol. 
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From the TG curves, DTG curves, Freidman Analysis, and activation energies, a single step 
autocatalytic kinetic model was deemed the best model to represent the PTFE degradation. The model 
fitted data can be seen in Figure 46 and the kinetic parameters can be seen in Table 11. Several other 
reactions were tested in order to validate the single step autocatalytic reaction. These reactions included 
alternative single step reactions as well as a few multiple step reactions. The closest fitting alternative 
models can be seen in Table 12. 

Comparing the single step autocatalytic reaction with the other reactions, an F-test shows that the 
chosen reaction fits the data best. See Table 12 for an F-test comparison and a schematic representation of 
the different models. The F-test is a statistical test used to compare models that have been fit to 
experimental data sets. The F-test identifies the best model fit of the population for which the data was 
sampled by testing the residual variances of the individual models against one another. The equation is as 
follows: 


or: 


Pexpi.fl.fr) = 


Z% 1 Zk=i{ Y jk~Yjk(rnodel 1 )') j 

£f=i£fe=i( i ’ jk-? jk(.model r i) / 

If, 


(46) 


F — test = 


lack-of-fit sum of squares/degrees of freedom 
pure-error sum of squares/ degrees of freedom 


(47) 


where is the degree of freedom for model 1 and f r is the degree of freedom for the reference model. 
Typically, the model with the lowest deviation, the best fit, is chosen as the reference model and 
therefore F exp = 1 . Based on the standard 5% significance level, F critical is determined from the 
F-distribution, the null distribution for the test. 
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Figure 44 Friedman Analysis for PTFE. 



Fract. Mass Loss 


o 

is 

> 

sr 



Fract. Mass Loss 


Figure 45 PTFE activation energy from (a) Friedman Analysis and (b) Ozawa-Flynn-Wall Analysis. 


64 



400 450 500 550 600 650 700 

Temperature (°C) 


Figure 46 Model of the best fit TG data for a single-step reaction model in air for PTFE. In the plot, the curves represent the 

model and the shapes represent the modeled experimental data. 


Table 11 Kinetic Parameters for the single step autocatalytic model for PTFE. 


Activation Energy (kJ/mol) 

Log (Al/s-1) 

Reaction Order 

Log (Kcat) 

169.0 

7.25 

0.74 

1.16 


Table 12 Comparison of kinetic models for PTFE*. 


Reaction Model 

Model Designation 

/(«) 

F 

1 exp 

First Order 

FI 

(1 - a) 

7.21 

Second Order 

F2 

(1 - a) 2 

19.34 

nth Order 

Fn 

(1 - a) n 

1.17 

nth Order with Autocatalysis 

Cn-X 

(1 - a) n (l + K cat X) 

1.0 


*F — 110 

r critical ^ ^ 
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Effect of thermal exposure on permittivity 

Experiment 

For the thermal exposure experiment, rectangular PTFE samples approximately 100 x 80 x 1 mm 3 were 
cut from the large PTFE sheet. Each sample was large enough to provide three smaller samples for 
subsequent permittivity measurements. The PTFE was heated isothermally at 340 ± 3 °C, a temperature 
above the melting point, in an isotemp muffle furnace for 2, 6, 12, 24, 36, 48, 60, 72, 84 and 96 hr. 

The PTFE samples for all exposure times were allowed to cool naturally in situ after the furnaces were 
turned off. Doors of the furnaces were not opened during the cooling process. While in the furnace the 
samples were sandwiched between two stainless-steel plates with mirror finish, to inhibit distortion. 
Once the samples were completely cool, after 48 hr, their permittivity was measured at room temperature 
over the frequency range 1 kHz to 2 MHz. Permittivity of as-received PTFE was also measured as 
baseline data. 

Result and discussion 

Figure 47 shows the way in which s' of PTFE increases as thermal exposure time increases, at four 
frequencies. Thermal exposure at 340 °C for 96 hr increases s' of PTFE by approximately 2%, which is 
significantly larger than the standard deviation between measurements on three nominally identical 
samples, which is less than 0.1%. 

The increase in real permittivity of PTFE after thermal exposure at 340 °C is attributed to enhancement 
of crystallinity following the heating and cooling processes. In semicrystalline polymers, motion of 
charge carriers, such as degenerated monomer molecules, can be blocked at the crystalline/amorphous 
phase boundaries and generate interfacial (Maxwell- Wagner-Sillars) polarization [74], If the crystalline 
phase grows, the interfacial polarization may be enhanced, causing an increase in s'. To investigate this 
hypothesis, X-ray diffraction (XRD) measurements were conducted on as-received PTFE and PTFE 
thermally exposed at 340 °C for 48 hr and 96 hr. After thermal exposure, three samples were allowed to 
cool slowly in situ as described in before. Three other samples were quenched in liquid N 2 to inhibit 
crystallization. While the XRD patterns of the PTFE samples all suggest a hexagonal crystal structure is 
present regardless of thermal history, the degree of crystallinity (as measured by peak area) in quickly- 
quenched PTFE was found to be 0.24 ± 0.01, versus 0.29 ± 0.01 in as-received PTFE, 0.32 ± 0.01 in 
PTFE slowly cooled after isothermal heating at 340 °C for 48 hr and 0.35 ± 0.01 in PTFE slowly cooled 
after isothermal heating at 340 °C for 96 hr. Hence the increase in s' as a function of thermal exposure 
time, shown in Figure 47, can be attributed to increasing crystallinity. The uncertainties derive from the 
standard deviation of measurements on three nominally identical samples. 

The increase in crystallinity observed in slow-cooled thermally-exposed PTFE compared with as- 
received samples supports the observation that s' increases following thermal exposure at 340 °C (which 
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is higher than the melting point 327 °C) with slow cooling (Figure 47). On the other hand, the reduction 
in crystallinity of quickly-quenched, thermally-exposed PTFE compared with as-received samples 
accompanied a reduction in s' of approximately 2%, much larger than the standard deviation between 
measurements of three nominally identical samples (0.2%). This observation further supports the 
hypothesis that increase in crystallinity is accompanied by an increase in s' for the range of conditions 
studied here. 

An alternative explanation for the increase in s' with increasing crystallinity of PTFE can be framed in 
terms of a charge carrier diffusion model. In reference [75] it is suggested that development of a 
heterogeneous structure such as exists in a semicrystalline polymer leads to an increase in permittivity of 
the polymer approaching the formation of a percolating network of the conducting amorphous phase. 

The two proposed models offer similar features since the time to explore the percolation structure by 
the charge carriers is equivalent to the time to transfer charges to the micro-capacitors formed by 
semicrystalline structure in the interfacial polarization model. 



Figure 47 Real permittivity of PTFE as a function of thermal exposure time at 340 °C in air. 


67 


4. Effect of tensile strain on permittivity 

Experiment 

The dielectric permittivity of untreated PTFE was firstly measured over the frequency range from 1 Hz 
to 1 MHz at room temperature using a Novocontrol Dielectric Spectrometer. As shown in Figure 39, the 
real permittivity of PTFE is approximately constant over this frequency range, with average value sr' = 
2.076 ± 0.002. For this reason, subsequent permittivity measurement of all samples in this study was 
conducted only at 1 00 kHz. The imaginary permittivity of PTFE is too small to be accurately measured by 
the dielectric spectrometer, so is not included in this discussion. 

In order to effectively study the influence of tensile strain on the permittivity of PTFE, a system capable 
of in-situ measurement of the permittivity of a polymer sample while under tensile strain was employed. 
This avoids the difficulty that once the polymer sample is removed from a tensile loading frame, the 
polymer chains may recover to some extent, changing the permittivity. As shown in Figure 48, the test 
fixture of an Agilent E4980A LCR meter fits in between the two clamps of a tensile loading frame 
(TestResources Inc., Model 150Q250) in order to measure the capacitance of the sample while it is 
stretched. For this reason, the accessible length of the sample was initially 75 cm or longer. Note that this 
configuration applies longitudinal stress while permittivity is measured orthogonally to the direction of 
the applied stress. The thickness of the samples was measured to uncertainty 1 pm by using a micrometer 
tipped with flat, parallel discs with diameter of approximately 5 mm that contact the sample. The diameter 
of the discs was selected to be similar to that of the capacitor electrodes with diameter of 5 mm, so that 
the thickness was measured over an area at least as small as the sample area. 

All the PTFE samples investigated in this study were cut from a single sheet supplied by Fluoro- 
Plastics Inc.. The dimension of each sample was approximately 90 x 40 x 0.39 mm3. An engineering 
stress-strain curve of one sample, shown in Figure 49, was measured in order to guide the selection of 
strains for simultaneous permittivity measurement. Engineering stress is defined as the ratio of the 
internal force with the original area of the sample surface on which the internal force acts on, in this case 
the cross-sectional area of the sample strip. The sample reveals elastic mechanical behavior below 18 
MPa load, and fails at approximately 31 MPa load, at which the engineering strain is approximately 190%. 
In order to ensure successful measurement without failure of the PTFE sample, mechanical strains of 25, 
50, 75, 100, 125 and 150% were selected for simultaneous permittivity measurement.. 
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Figure 48 Experimental arrangement for permittivity measurement while the sample is under tensile strain, using an Agilent 
E4980A LCR meter and Test Resources, Inc. tensile load frame Model 150Q250. 

To explore the potential effect of strain rate, strains of 50, 100, and 150% were achieved by strain rates 
of 12.7 and 25.4 min/min, respectively. Other strains studied were achieved by strain rate 25.4 mm/min. 
Once the assigned strains were reached, dielectric permittivity of all stretched samples was measured at 
their centers, at which they were most reliably stretched to the desired strain with assigned strain rate. The 
initial dielectric permittivity of each sample was measured at 100 kHz prior to application of tensile strain. 

Additionally, to study the effect of stress relief, samples with strains of 50, 100, and 150% were held in 
the tensile loading frame for up to 24 hours, and their in-situ permittivity was measured at times t = 1, 2, 4, 
6, and 24 hours, respectively. Finally, the load was released and permittivity of all the samples was 
immediately measured again, to explore the effect of recovery of the polymer chains. 

X-ray Diffraction (XRD) was conducted to examine possible structural changes in PTFE as a 
consequence of tensile strain, in order to explore the mechanism by which mechanical strain influences 
the permittivity of PTFE. 
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Figure 49 Engineering stress-strain curve of PTFE. 


Results and discussion 

Figure 50 shows the difference between the real relative permittivity of strained and unstrained PTFE, 
Asr, as a function of mechanical strain. The real permittivity of PTFE decreases with increasing tensile 
strain, e.g. by approximately 19% for 150% strain. As stated in Chapter II, a dielectric material exhibits 
some or all of the four polarization mechanisms: atomic, ionic, dipolar and interfacial polarization [74], 
Given that PTFE is a non-polar semicrystalline polymer, its real permittivity is associated with only two 
polarization mechanisms: atomic polarization due to a shift of the electron clouds in each atom under the 
influence of the applied electric field and interfacial polarization due to accumulation of free charges at 
interfaces between the amorphous and crystalline polymer phases [74], The observed decreases in the 
real permittivity of PTFE therefore may in principal be as a result of decrease in either atomic 
polarization or interfacial polarization, or both. Atomic polarization as defined above, however, cannot be 
influenced by mechanical strain. On the other hand, according to the Fringed Micelle Model [76], the 
chains of a polymer under mechanical load increase their alignment along the direction in which the load 
is applied, hence increasing the degree of alignment in the amorphous phase and producing a greater 
degree of crystallinity. Thus, it is hypothesized here that the interfacial polarization of PTFE declines 
under mechanical strain due to decreased mobility of free charges, as a result of increased ordering of the 
polymer chains. 
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To investigate this hypothesis, the degree of crystallinity of samples with tensile strains of 50, 100 and 
150% was determined by XRD, immediately upon removal of the polymer samples from the load frame. 
It was observed that the degree of crystallinity increases, overall, with mechanical strain, as listed in 
Table 13. The XRD spectra also revealed that the crystal structure of strained PTFE at room temperature 
is hexagonal, the same as that of un-strained PTFE [64] . 

After the samples were released from the tensile loading frame, their permittivity was immediately 
measured and observed to recover somewhat, increasing by up to 8% as shown in Figure 51, suggesting 
partial recovery of the polymer chains. It is reasonable to conclude that the decrease in the real 
permittivity of PTFE is associated with increased ordering of the polymer chains in the presence of 
mechanical strain. 

The potential effect of stress relief was investigated by measuring permittivity on samples strained to 50, 
100 and 150% for 24 hours. Permittivity was measured at times t = 1, 2, 4, 6, and 24 hours. No 
significant change in measured permittivity was observed as a function of time up to 24 hours. Finally, 
the permittivity of PTFE stretched with strain rate of 25.4 mm/min was observed to be quantitatively 
similar to that of PTFE stretched to the same strain with half the strain rate, 12.7 mm/min. Thus, within 
the range of stresses and strains studied here, it is concluded that there is no obvious influence of stress 
relief and strain rate on the permittivity of strained PTFE. 

Table 13 Crystallinity of PTFE for various values of tensile strain, measured by X-ray diffraction immediately upon removing 

the sample from the tensile tester. 


Tensile strain (%) 

Degree of crystallinity (%) 

0 

36 ±3 

50 

42 ±3 

100 

42 ±3 

150 

47 ±3 
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Figure 50 The change in the real relative permittivity of PTFE as a function of tensile strain. Note, g r = 2.076 at 0% strain. 



50 60 70 80 90 100110120130140150 


Strain (%) 

Figure 51 The difference between the real relative permittivity of PTFE under strain (solid symbol) and of released PTFE (open 

symbol), compared with untreated PTFE. 
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5. Summary 

Thermal degradation of PTFE is considered as a single step autocatalytic reaction. Thermal exposure 
studies on PTFE show increase in real permittivity by approximately 2% following isothermal heating at 
340 °C for 96 hr. This observation is associated with increase in crystallinity (from 29% in an as-received 
sample to 35% in a sample isothermally heated at 340 °C for 96 hr) that enhances the interfacial 
polarization between the amorphous phase and crystal phase of the polymer. On the other hand, in the 
presence of mechanical loading, polymer chains of PTFE are oriented along the direction of loading, 
which limits motion of free charges and thereby weakens the interfacial polarization. Consequently, the 
real permittivity of PTFE measured transverse to the direction of the applied load was observed to 
decrease with applied load; by approximately 19% for engineering strain of 150%. An associated increase 
in the degree of crystallinity with increasing tensile strain was measured (from 36% in an as-received 
sample to 47% in a sample with engineering strain 150%), supporting the hypothesis that limited motion 
of free charge carriers is responsible for the observed reduction in permittivity transverse to the applied 
load direction. Upon removal from the load frame, the oriented polymer chains of PTFE recovered to 
some extent, as indicated by an associated increase in the permittivity by up to 8%. On the other hand, the 
strain rate by which the PTFE samples were stretched to various strains and the stress relief process 
investigated while the PTFE samples were held at constant strain for up to 24 hr did not change the real 
permittivity of PTFE significantly. 

It is interesting that, in the case of thermal exposure, an increase in crystallinity of approximately 6% is 
observed along with a modest increase in permittivity (up to 2%) in thermally exposed samples for the 
conditions studied whereas, for strained samples, an increase in crystallinity of approximately 1 1 % is 
observed along with a larger (up to 19%) reduction in permittivity. It should be noted that value of 
crystallinity is not the primary factor in determining interfacial polarization relaxation intensity (and 
therefore permittivity). For example, a study on poly( ether ether ketone) [77] concluded that the nature 
and degree of perfection in the crystallized regions play a more dominant role in determining the 
dielectric properties of the polymer. The degradation mechanisms studied here (thermal exposure and 
mechanical stress) likely give rise to crystallizations that are significantly different from each other in 
nature and degree of perfection. Thus, the trends in the measured values of crystallinity and in the 
permittivity for each of the degradation mechanisms described above are likely not in conflict but rather 
explained by significant differences in the crystal nature of those samples. 
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Chapter IV. Ethylene-tetrafluoroethylene (ETFE) 

ETFE (ethylene-tetrafluoroethylene) has been applied as the only wiring insulation material in Boeing 
747, 767 and 777 aircraft since the 1990s. This chapter describes an investigation of the way in which the 
real permittivity and loss factor of ETFE change after thermal exposure in air for up to 96 hours. 

1 . Introduction 

Ethylene-tetrafluoroethylene (1:1 ratio) alternating copolymer (-CH 2 CH 2 CF 2 CF 2 — )„ is a fluoroplastic 
material. At room temperature, it has a molecular conformation in which extended zigzag chains are 
packed in orthorhombic cells with cell dimensions a = 0.96 nm, b = 0 nm, c = 0.50 nm and y = 96° [57]. 
Each molecule is adjacent to four other molecules in which the — CH 2 - groups of a chain are positioned 
next to the -CF 2 - groups of the next chain [78], Essentially, the bulky -CF 2 - groups nestle into the space 
above the smaller -CH 2 - groups of an adjacent chain [79], interlocking the chains as a result. Due to this 
special chain conformation, ETFE has desirable properties such as a stable dielectric constant, low 
dissipation factor, high melting temperature and high elastic modulus. Selected physical and electrical 
properties of ETFE taken from the literature are listed in Table 14. 

Reference [57] has presented the dielectric constant of ETFE over the frequency range from 1 kFIz to 
10 GFIz at 22 °C and 150 °C. At room temperature, the dielectric constant of ETFE is between 2.5 and 2.6 
below 10 MFIz and decreases to approximately 2.3 at 10 GHz. Change in dielectric constant due to 
temperature can also be observed, especially at high frequencies. The dissipation factor of ETFE is less 
than 0.001 below 10 kHz, increasing to 0.023 at 100 MHz and decreasing at higher frequencies [57]. 


Table 14 Properties of ETFE. 


Property 

Value 

Property 

Value 

Density (g/cm3) 

1.75-1.79 

Dielectric constant 

2. 5-2. 6 below 
10 MHz 

Melting Temperature (°C) 

255-280 

Dissipation factor 

Less than 0.007 below 1 MHz 

Molecular relaxation 
temperatures (°C) 

-110(a) 
0-100 (a') 
~-87(y) 

Dielectric strength (kV/mm) 

64 (254 pm film) 
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2. Sample material 

The ETFE samples studied here were discs nominally 1 mm thick and 63 mm in diameter, which were 
extruded from Tefzel 750® ETFE pellets supplied by Dupont. The heat flow and storage modulus as a 
function of temperature on samples cut from different spots of the extruded discs were measured by DSC 
(Differential Scanning Calorimeter) and DMA (Dynamic Mechanical Analysis), respectively. Good 
homogeneity in thermal and mechanical properties of the extruded samples was verified, as shown in 
Figures 52a and 52b. 

The storage modulus E', loss modulus E" and tan 8 of extruded ETFE from -150 °C to 180 °C at 1 Hz 
was measured by DMA, as shown in Figure 53. The storage modulus is 868 MPa at 25 °C, and 57.5 MPa 
at 150 °C. The tan 8 curve reveals three molecular relaxations at 86.0, approximately 50, and -98.1 °C, 
analyzed by the software ‘TA Universal Analysis’. The relaxation at 86.0 °C is the a-relaxation, which 
has been attributed to the glass transition that occurs in the amorphous phase [17]. The relaxation at 
approximately 50 °C is attributed to the a'-transition that is related to a reversible crystal transformation 
from the orthorhombic to the hexagonal lattices. This transformation occurs gradually over a broad 
temperature range (0 to 100 °C). In some ETFE, this transition is ill-defined and can be better described 
as a transition from pseudohexagonal to a hexagonal crystal phase [12]. The peak at -98. 1 °C is associated 
with the y-transition. 

The real permittivity and loss factor of extruded ETFE was measured over the frequency range 1 Hz to 
1 MHz and the temperature range -150 °C to 210 °C by the Novocontrol Dielectric Spectrometer. Three 
dimensional dependences of s' and tan 8 are shown in Figures 54a and 54b, respectively. The highest 
values of s' and tan 8, observed at low frequency (below ~10 Hz) and high temperature (above ~200 °C) 
may be due to one or more extrinsic or intrinsic electrical phenomena, such as electrode polarization, 
interfacial polarization and conductivity effects [80]. 

Extracting data from Figure 54, the dependence of s' and tan 8 on temperature at 1.15 kHz are shown in 
Figure 55. Two peaks in the loss factor of ETFE at -80 °C and 92 °C are observed. These are associated 
with y- and a-transitions respectively [81]. The two peak temperatures are close to the transition 
temperatures listed in Table 14, -87 °C and 110 °C. The differences in these values are likely due to 
measurement of different material properties by different equipment on slightly different samples. 
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Figure 52 Homogeneity of extruded ETFE by (a) DSC; (b) DMA. 
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Figure 53 Results of dynamic mechanical analysis on extruded ETFE. 
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(b) 



Figure 54 Real permittivity (a) and loss factor (b) of extruded ETFE as a function of frequency and temperature. 
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Figure 55 Real permittivity and dissipation factor of extruded ETFE as a function of temperature at 1.15 kHz. 


3. Thermal degradation 

Thermal degradation kinetics 

Method and experiment 

The kinetics of thermal degradation of ETFE was studied by using the method for studying PTFE. The 
extruded ETFE was cleaned prior to experimental testing with a Kim Wipe and ethanol to remove any 
oils or dusts on the sample that might affect the weight and degradation of the material. Each test sample 
was removed from the disk through a circular punch ensuring reproducible sample weight and shape. The 
sample masses were 8.7 ± 0.2 mg. To reduce sample variability, the samples were chosen from the lower 
third of the extruded disk. 

Expounding the kinetic analysis from a single-step reaction (described in Chapter III) to a consecutive 
multistep reaction, the differential equations are as follows: 

^ = -f(a,b)A t e~ E t /RT (48) 

T t = -T t -mOA 2 e~^ (49) 


c = 1 — a — b = a 


(50) 
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where the reaction follows an A B C model [82] [83]. The rate of reaction for the degradation from 
A -> B (step 1) is given by da/dt. The rate of reaction for the degradation from B C (step 2) is given 
by db/dt. In this format of differential equations a, b, and c are analogous to concentrations in chemical 
kinetics and A, B and C represent degradation states; where state A corresponds to a = 1, b = 0 and 
c = 0, state B corresponds to a = 0, b = 1 and c = 0, and state C corresponds to a = 0, b = 0 and c = 1. 
The degradation continues to follow the analogy of chemical kinetics, where step 2 follows step 1, but 
may occur before complete conversion of a to b. 

Results and discussion 

The TG scans for ETFE began at room temperature, and degradation was not observed in either system 
until elevated temperatures were reached. The high temperature TG data is given in Figure 56. Like most 
fluoropolymers, ETFE is extremely stable at intermediate temperatures. The onset of degradation 
increases with increasing heating rate and involves a rapid and complete degradation. From the TG curves, 
the variation between the temperatures of the degradation onset points provides insight into the reaction 
model. For ETFE, the separation between the temperatures of the degradation onset points increases 
between curves as the heating rate increases. This indicates a consecutive step reaction model. Whereas, a 
decreasing separation between the degradation onset points as the heating rate increases indicates a 
competitive model. 

The derivative of the weight with respect to temperature provides better insight into the mechanism of 
degradation. For a specific heating rate, the number of peaks in the derivative thermogram (DTG) 
represents the minimum number of reaction steps involved. For ETFE, it is clearly seen at 2 Kmin 1 that 
there are two distinct peaks. At slow heating rates, the reaction is dominated by a specific mechanism; 
however, as the heating rate increases, this dominating mechanism is replaced with an alternative 
mechanism at higher heating rates. Figure 57 shows the two mechanisms in the DTG curves. Based on the 
DTG curves below, heating at 2 Kmin' 1 , a minimum of two reaction steps, or two peaks, can be seen. 

Once initial observations of the thermograms have been made, a Friedman analysis based on equation 6 
provides information about multi-step processes during the reaction and also offers insight into the type of 
reaction step. Similar to the DTG plots, the number of peaks indicates the minimum number of reaction 
steps in the degradation. For ETFE there are two distinct overlapping peaks seen in Figure 58. Therefore a 
multi-step reaction is again evident from the curvature of the plot for ETFE. 

The type of the reaction can also be determined by comparing the slope of the constant fractional mass 
loss trend line to the slope of each peak at the beginning of the reaction step [21]. For ETFE, the first 
peaks have in all cases the acceleration, probably autocatalysis. For autocatalysis, the generic governing 
differential equation, presented in equation 2, defines the reaction model, f(a), such that: 
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f(a) = (1 - a) n (l + K cat X) (51) 

where a represents the degree of degradation, n represents the reaction order, K cat represents the 
autocatalysis constant, and X represents the step of the reaction. 

From the Friedman analysis as well as the Ozawa-Flynn-Wall analysis, activation energy as a function 
of fractional mass loss can be obtained [18] [19] [27]. The plot of the activation energy for ETFE with 
respect to the amount of degradation again confirms the multistep reaction by presenting fluctuating 
activation energy seen in Figure 59. This indicates an overlap of multiple reactions. Activation energy 
developed from Friedman analysis presents activation energies between 140 kJ/mol and 210 kJ/mol. 

To verify the two-step model the TG curves, DTG curves, Freidman Analysis, and activation energies 
suggested, a two-step consecutive model was introduced to (a) the two slowest heating rates and (b) the 
two fastest heating rates — seen in Figure 60. Table 15 presents the corresponding model parameters. The 
difference between the parameters, primarily the activation energies, indicates the two-step model is an 
insufficient model for predicting thermal degradation and indicates another step must be involved. 
Additionally, the model fit of the slow heating rate DTG seen in Figure 61 identifies a poor fit at the end 
of the reaction and this difference between model and experimental data indicates an additional step 
required. 

Consequently, a three-step consecutive model was used. The model fit is presented in Figure 62. Table 
16 presents the corresponding model fitting parameters. The activation energy for the first step in the final 
model is the same as the activation for the two-step fit of the slow thermograms. This means that the first 
step of the final model is the deciding step for the initial part of the degradation reaction. The second step 
is a secondary process for that deciding step. This secondary process is responsible for large temperature 
separation of the degradation onset points at high heating rates in Figure 56. This second step, a coupling 
step, should not be considered separately from the first step. The third step in the three-step consecutive 
model is responsible for the second DTG peak, or the final stages of the degradation process for both the 
slow and fast heating rates. 

However, comparing the three step consecutive reaction to three-step competitive reactions, an F-test 
shows that the consecutive reaction fits the data better. See Table 17 for F-test comparison and a 
schematic representation of the different models. 
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Figure 56 TG curves for ETFE. 
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Figure 57 DTG curves for ETFE. 
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Figure 58 Friedman Analysis for ETFE. 




Figure 59 Activation energy from (a) Friedman Analysis and (b) Ozawa-Flynn-Wall Analysis. 
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Figure 60 Two-step consecutive model fits to (a) the two slowest heating rates and (b) the two fastest heating rates. 


Table 15 Reaction parameters for two-step consecutive model fits for ETFE. 



Activation Energy (kJ/mol) 

Log (Al/s-1) 

Reaction Order 

Log (Kcat) 

OX 

£ .S 
5 +* 

Step 1 

246.26 

16.47 

1.97 

1.26 

.= C5 

czi « 

3 

Step 2 

181.84 

11.02 

0.96 

- 0.56 

Fast 

Heating 

Step 1 

79.30 

1.34 

0.71 

3.15 

Step 2 

103.52 

5.48 

1.02 

- 4.00 
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Figure 61 DTG model fit of the two slowest heating rates for ETFE. 
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Figure 62 Model of the best fit TG data for the three-step reaction models air atmospheres. In the plots, the curves represent 
the model, the data point shapes represent the modeled experimental data, and the + represents the un-modeled experimental data. 
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Table 16 Reaction parameters for the three-step consecutive model fit. 



Activation Energy (kJ/mol) 

Log (Al/s-1) 

Reaction Order 

Log (K cat ) 

Step 1 

247.25 

16.55 

1.48 

-0.59 

Step 2 

0.03 

-3.50 

0.74 

2.34 

Step 3 

123.03 

6.38 

1.37 

0.63 


Table 17 F-Test statistical analysis of the model fits for three-step reactions for ETFE*. 


Consecutive 

Reaction 

Competitive 
Reaction 1 

Competitive 
Reaction 2 

Competitive 
Reaction 3 

A — ► B — ► C — ► D 

A — > B — C 


A — r> B 

U- C > D 

1.0 

1.78 

11.69 

57.24 


* ^critical — 1-02 

Effect of thermal exposure on permittivity 

Experiment 


Prior to thermal exposure, extruded ETFE discs were annealed isothermally at 60 °C in another isotemp 
muffle furnace for 2 hr prior to thermal exposure. The annealing was done in order to relieve residual 
stress that resulted from the extrusion process and flatten the samples for permittivity measurement. 
When cooled completely, annealed discs were subjected to thermal exposure at 160 ± 2 °C for 2, 6, 12, 24, 
36, 48, 60, 72, 84 and 96 hr. The temperature 160 °C, much lower than the melting temperature of 245 °C, 
was chosen due to significant sample deformation after thermal exposure at higher temperatures that 
prevented successful permittivity measurement. Three discs were exposed for each exposure time. 

During the exposure process, all the samples were sandwiched between two stainless-steel plates with 
mirror finish, to inhibit distortion. After the furnace was turned off, doors of the furnaces were not 
opened until all exposed ETFE samples were naturally cooled down in situ. Once the samples were 
completely cool, after 48 hr, their permittivity was measured at room temperature over the frequency 
range 1 kHz to 2 MHz. Permittivity of as-extruded ETFE and annealed ETFE samples were also 
measured to provide comparative data. 

Result and discussion 

Values for s' and tan 5 of both extruded and annealed ETFE samples are compared in Figure 63. The s' 
of extruded ETFE decreases from 2.599 to 2.539 as frequency increases from 1 kHz to 2 MHz. Loss 
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factor, tan 8, increases over this frequency range but is small (below 0.007) even at 2 MHz. These results 
agree with the values listed in Table 13. It can be observed that s' of extruded ETFE and annealed ETFE 
are similar over the frequency range measured and tan 8 overlaps closely until above 1 MHz. (The 
separation below 20 kHz is due to large uncertainty in the measurements for low frequencies.) It can be 
concluded that the annealing process at 60 °C for 2 hr does not influence the dielectric properties of ETFE 
significantly and it is reasonable to anneal the samples prior to thermal exposure at 1 60 °C for up to 96 hr, 
for the purpose of relieving residual stress. 

Figures 64 and 65 show how thermal exposure duration influences s' and tan 8, respectively, of ETFE at 
four frequencies. Since the measurement uncertainty in tan 8 is quite large below 20 kHz, tan 8 at 1.00 
and 10.6 kHz are not shown in Figure 65. Figure 64 shows how thermal exposure time influences s' of 
ETFE at four frequencies. As exposure time increases up to 96 hr, s' increases by 0.05 (approximately 
2%, which is significantly larger than the standard deviation between measurements on three nominally 
identical samples, 0.5%). From Figure 65, tan 8 is not changed significantly by thermal exposure at 1 12 
kHz. At 1.18 MHz, however, tan 8 decreases by 1 1 % after thermal exposure for 96 hr. 

XRD tests conducted on annealed ETFE and ETFE thermally exposed at 160 °C for both 48 hr and 96 
hr indicated that the degree of crystallinity, 0.65 ± 0.01, was not changed by thermal exposure. This is in 
accordance with results presented in reference [84], 

On the other hand, mid-infrared (IR) analysis on annealed ETFE and ETFE thermally exposed at 
160 °C for 96 hr has revealed chemical changes in ETFE following thermal exposure. The absorbance 
spectra are shown in Figure 66. Similar peaks corresponding to the presence of CF 2 at 1460 cm' 1 , and C- 
H at 2850 cm 1 and 2890 cm 1 are observed in both spectra. However, the spectrum of thermally-exposed 
ETFE contains an extra peak at 1680 cm' 1 , corresponding to the C=0 group, and a much higher and 
broader -OH peak at 3450 cm' 1 . These indicate that chemical changes occurred during thermal exposure 
of ETFE, including oxidation and dehydrofluorination [85], Oxidant thermal exposure accompanied by 
chain breakage and polar group formation can increase dielectric permittivity of polymers by enhancing 
their charge carrier mobility. It has been reported [86] that chain breakage can decrease polymer 
viscosity by reducing average molecular weight, which increases the charge carrier mobility. The 
formation of polar groups can also promote the charge carrier mobility in ETFE. The decrease in tan 8 at 
1.18 MHz, shown in Figure 65, is also attributed to the additional polar groups formed during thermal 
exposure [87]. 
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Figure 63 Real permittivity and dissipation factor of extruded and annealed ETFE as a function of frequency at room 

temperature. 



Figure 64 Real permittivity of ETFE as a function of thermal exposure time at 160 °C. 
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Figure 65 Dissipation factor of ETFE as a function of thermal exposure time at 160 °C. 



Wave number (cm' 1 ) 


Figure 66 Mid-IR spectra of annealed ETFE and ETFE thermally exposed at 160 °C for 96 hr. 


4. Summary 

The thermal degradation of ETFE is considered as a three-step consecutive reaction model with each 
step governed by an autocatalysis reaction model. After thermal exposure at 160 °C for 96 hr in air, 
approximately 2% increase in the real permittivity and 11% decrease in the loss factor of ETFE are 
observed. These observations are considered as a result of oxidation and dehydrofluorination during the 
heating process, which give rise to formation of polar groups. 
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Chapter V. Development of a Capacitive Sensor for Nondestructive 

Evaluation of Wiring Insulation 

Curved patch capacitive sensors have been developed to evaluate the material dielectric properties of 
cylindrical structures. These sensors consist of two curved patch electrodes that are coaxial with and 
exterior to the cylindrical test-piece. Description on the development of capacitive NDE sensor in this 
chapter is divided into three parts. In part 1, a numerical model that relates the permittivity of a 
homogeneous dielectric rod to the capacitance of a curved patch sensor is described. Experimental results 
that verify the validity of the numerical model are provided. In part 2, a numerical model that relates the 
permittivity of a cylindrical dielectric-coated conductor to the capacitance of a curved patch sensor is 
described, together with benchmark experiment results that verify the theory. In the last part of this 
chapter, a prototype capacitive sensor has been fabricated and applied for dielectric property 
measurements of aircraft wires (wire type MIL-W-81381/12). Groups of wires have been exposed to 
different thermal and hydrolytic environments. Experimental studies using the prototype capacitive sensor 
showed that, thermal and hydrolytic exposure induced dielectric property changes in this type of wire 
insulation have been successfully detected. 

1 . Analysis of a capacitive sensor for the evaluation of circular dielectric cylinders 

An arc-electrode capacitive sensor has been developed for the quantitative characterization of 
permittivity of cylindrical dielectric rods. The material property of the cylindrical test-piece can be 
inversely determined from the sensor output capacitance based on a theoretical model. For the modeling 
process, the electrostatic Green's function due to a point source exterior to a dielectric rod is derived. The 
sensor output capacitance is calculated numerically using the method of moments (MoM), in which the 
integral equation is set up based on the electrostatic Green's function. Numerical calculations on sensor 
configuration optimization are performed. Calculations also demonstrate the quantitative relationship 
between the sensor output capacitance and the test-piece dielectric and structural properties. Capacitance 
measurements on different dielectric rods with different sensor configurations have been performed to 
verify the validity of the numerical model. Very good agreement (to within 3%) between theoretical 
calculations and measurement results is observed. 

Introduction 

Increasing demands for dielectric measurements have been observed over the past decade, with new 
applications of advanced composites in modem aircraft, automobiles, and shipbuilding. Specifically, 
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dielectric measurements are important for the characterization of thin films, substrates, circuit boards, 
printed-wiring boards, bulk materials, powders, liquids and semisolids. 

Capacitance methods, because of their simplicity, relatively low cost, and high accuracy, have been 
applied to characterize the dielectric properties of many different materials. Over the past 100 years, 
closed-form solutions for capacitances due to various canonically shaped electrodes have been found, by 
mapping out the electrostatic field in the vicinity of the conductors. Canonical electrode shapes are those 
formed from surfaces easy to describe in standard coordinates, including strips, circular discs, annular 
rings, cylindrical arcs, spherical caps, etc. 

It is convenient to solve capacitance problems associated with electrified strips using Cartesian 
coordinates. The capacitances for two parallel and coplanar infinite strips [106], as well as charged thin- 
strip quadruples [107], have been solved using the triple integral equations. In addition, the potential 
associated with a physically more realistic strip of finite length, the potential due to polygonal plates, and 
the potential due to a charged elliptical plate have been derived using dual integral equations in Cartesian 
coordinates [107]. 

Using cylindrical coordinates, the solution to dual integral equations has been applied to obtain the 
surface charge distribution of a charged disc in free space, and also to obtain the potential due to a circular 
disc placed between two parallel earthed planes [106]. The solution of the Fredholm equation has been 
applied to solve for the capacitance of an electrified disc situated inside an earthed coaxial infinitely long 
hollow cylinder. Also, the field due to two equal coplanar electrified discs has been solved by the method 
of Kobayashi potentials, while the capacitance between two identical, parallel and coaxial discs has been 
obtained by solving Love's integral equation. These solutions are available in the classic book written by 
Sneddon [106]. Furthermore, an axisymmetric problem of several charged coaxial discs has been 
considered by the dual integral equation method, and the solution has been obtained for the case when the 
distance between neighboring discs is large compared to their radii [108]. The potential of a system of N 
charged, arbitrarily located, circular discs has also been considered in [108]. Aside from classic problems 
associated with discs, Cooke's solution of a set of triple integral equations has been applied to solve for 
the field due to a charged annular ring of finite width in free space [106], while the solution to the integro- 
series equations has been applied to obtain the total charge for a capacitor that consists of a coupled disc 
and spherical cap [107]. 

Another set of canonical capacitance problems discussed in the literature are infinitely long axially 
slotted open cylinders. In [107], the capacitance generated by a pair of charged symmetrically or 
asymmetrically placed circular arcs has been calculated in circular cylindrical coordinates, whereas the 
capacitance due to a pair of charged symmetrically placed elliptic arcs has been solved in elliptic cylinder 
coordinates [107]. 
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The canonical capacitance problems, mentioned above, are all discussed in free space situations, and 
need further modifications to be applicable for characterization of material dielectric properties. Other 
semi-analytical and numerical capacitive solutions have been developed over the past decade to keep pace 
with new applications of dielectric materials. For example, interdigital dielectrometry has been applied 
for material dielectric property characterization as one of the most commonly used capacitance sensors. 
An excellent review paper on interdigital sensors is [109], in which sensor modeling, fabrication, 
measurement setup, and applications are discussed in detail. In addition to a widely-used effective semi- 
analytical approach, called the continuum model [110], conformal mapping techniques have also been 
applied to obtain closed-form solutions for the capacitance of interdigital sensors in surface contact with 
multi-layered dielectric structures [111]. Examples of practical applications of interdigital sensors include 
estimation of meat fat content [112] and insulation damage detection in power system cables [113]. 
Multichannel fringing electric field sensors, designed by finite -element (FE) method calculations for 
sensor modeling, optimization and performance evaluation, have been used for material property 
measurements [114], Cylindrical geometry electroquasistatic dielectrometry sensors have been developed 
using semi-analytical models to quantitatively relate the dielectric properties of multi-layered test-pieces 
to sensor output transcapacitance [115]. Concentric coplanar capacitive sensors for nondestructive 
evaluation of multi-layered dielectric structures have been developed in [116], and can be applied to 
detect water ingression and inhomogeneities in aircraft radome structures. In addition, rectangular 
coplanar capacitance sensors have been developed to detect water intrusion in composite materials [117] 
and for damage detection in laminated composite plates [118]. Rectangular capacitive array sensors that 
detect surface and subsurface features of dielectric materials have been developed in [119]. In [99], 
approaches of determining the length of open-circuited aircraft wires through capacitance measurements 
have been presented. 

In this research, a model-based capacitive method is developed for the quantitative dielectric property 
characterization of circular cylindrical dielectric rods. The work is motivated by testing of cylindrical 
components such as wiring insulation or polymeric tubing, and will be developed to deal with those cases 
in future. The capacitance sensor consists of two arc-shaped patch electrodes that are located exterior to 
and coaxially with the cylindrical test-piece. These two sensor electrodes exhibit a measurable 
capacitance whose value depends on both the dielectric and geometrical properties of the dielectric rod. 
The arc-electrode configuration offers a nondestructive and convenient way of determining the dielectric 
constant of cylindrical test-pieces, compared to cutting a slice from the test-piece for a conventional 
parallel-plate capacitance measurement. A numerical method, the method of moments (MoM), is 
employed in the numerical calculations. First, the Green's function for a point source over the surface of a 
dielectric rod is derived in cylindrical coordinates, in the form of modified Bessel functions of the first 
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and second kinds of integer order n : / n (z) and K n (z) . This Green's function then serves as the 

integration kernel in MoM calculations, from which the sensor surface charge distribution is obtained. 
Once the sensor surface charge distribution is known, one can easily calculate the sensor output 
capacitance C through C = Q / V , where Q is the total charge on one electrode and V is the potential 
difference between the two sensor electrodes. Validation of numerical results by benchmark experiments 
has been performed, and very good agreement (to within 3%) between theoretical calculations and 
measurements is observed. The advantage of the arc-electrode capacitive sensor reported in this study, 
compared to existing planar capacitive sensors, is that the arc electrodes conform to the surface of a 
cylindrical test-piece and concentrate the electric field in the material under test. In addition, the physics- 
based model developed in this report allows inverse determination of test-piece permittivity from 
measured arc-electrode capacitance. 

Modeling 

Figure 67 shows the configuration of the arc -electrode capacitive sensor. The capacitive sensor consists 
of two identical arc-electrodes coaxial with a cylindrical dielectric rod, and exhibits a measurable 
capacitance C that is quantitatively related to the permittivity and diameter of the material under test. In 
the theoretical modeling, the cylindrical dielectric rod is assumed to be infinitely long while the arc- 
electrodes are infinitesimally thin. The more general case in which the electrodes and the test-piece have 
different radii, as shown in Figure 67, is considered theoretically. Flowever, in order to achieve maximum 
output capacitance, it is more desirable to have the sensor electrodes in tight contact with the cylindrical 
test-piece in measurements. One practical approach to achieving tight surface contact between the 
electrode and the test-piece is to deposit the sensor electrodes on a compressible dielectric material used 
as the sensor substrate, and press the substrate against the test-piece to conform the arc electrodes to the 
test-piece surface. This approach will be attempted in a future version of the sensor. 


93 


Dielectric rod 




Figure 67 Arc-electrode capacitive sensor. The radii of the sensor electrodes and the cylindrical dielectric rod are denoted /7g 
and a, respectively. The arc-angle of each sensor electrode is <f > 0 (rad). The length of each electrode in the vertical direction is / 

and the width in the horizontal direction is W — (f)§ X p Q . 



Figure 68 Point source outside of a dielectric rod, assumed infinitely long. 


Derivation of Green's function in cylindrical coordinates 

The electrostatic Green's function due to a point charge outside of an infinitely long dielectric rod is 
derived in cylindrical coordinates, to form the integral equations later used in MoM calculations. Figure 
68 shows a point charge placed at (p' ,(f)' ,z'^ exterior to a cylindrical dielectric rod of radius a and 
dielectric constant s 2 . Without loss of generality, the dielectric constant for the medium exterior to the 
dielectric rod is assumed to be Si. The resulting potential v F(/7,^,z) due to such a point charge satisfies 
the Laplace equation in each homogeneous region: 
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and is subject to the interface conditions at the surface defined by p = a \ 
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where superscripts (52) and (53) correspond to the regions defined by p > a and 0 < p <a , respectively. 
To find a suitable solution, one starts with the fundamental solution of the Laplace equation in free space: 

1 


G 0 ( r l r ') : 


(54) 


An\r-r'\’ 

where r corresponds to the observation point at {p,(j),z) and r' corresponds to the source point at 
(p' ,(/)' ,z'). Using the integral [120] 


1 2 


r tc 


poo 

K 0 (/cp)cos(/cz)d/c, 

J 0 


(55) 


where r = \j p + z" and K (] (kp) is the modified Bessel function of the second kind of order zero, one 

can rewrite the fundamental solution in cylindrical coordinates, which is convenient considering the 
cylindrical boundary surfaces of the problem. By using the addition theorem [120] 


*„(«?) = I o( k: P<) K o (^>) + 2X cos W^-^)]A(^<)^, (^>), (56) 

/=1 

where g = p 1 + p' 2 - 2 pp' cos (</> - (/>') , p < is the lesser of p and p , and p y is the greater, (55) is 
transformed to an arbitrary coordinate system. The potential at an observation point ( p,(f),z ) due to a 
source point at (p', </>', z') in free space is expressed as 

\ ( poo ~ 

G °( r ' = ~1k 2 X I Jo °^ P ' P ’’ K J cos ^ z “ z '^ dK 


0° ^ 1 

+2^ cos[/(^ - 0')] J G, ( p , p' , k) cos[/e(z - z')\dic 

t = i 0 J 


(57) 


where 


G, (p,p',K) = I t (icp < )£,(*?>) / = 0,1,2,---, 


(58) 


95 


and I^iep^ is the modified Bessel function of the first kind of order / . It is noted that the difference 

between the Green's function due to a point charge in free space and that due to a point charge exterior to 
an infinitely long dielectric rod lies in the interface conditions at a surface p = constant. To find the 
Green's function due to a point charge outside of a dielectric rod, one needs only to modify the integral 
kernel G t (p, p' ,re) in (57) so that the interface conditions in (53) are satisfied: 


Gf\p,p',K) = Gf\p,p',K), (59) 

and 

Si — G? ) {p,p',K) = £ 2 — G ( t 2) {p,p\K) (60) 

op dp 

To create a kernel that satisfies (59) and (60), define 

Gf\p, p\ re) = I, (Kp K )K t (rep y ) + A(re)K t 0 Kp)K t (rep'), (6 1 ) 

and 

G, (2) (p, p',k) = B(ie)I t ( rep)K t (rep'). (62) 

The first term in (61) represents the primary field due to the point source, while the second term 
represents the reflected field with a reflection coefficient A(re) . In the region 0 < p <a there exists only 
the transmitted field and B(k) in (62) is the transmission coefficient. Substitute (61) and (62) into (59) 


and (60) to find that the coefficient A(k) is expressed as 


4*0 


s 2 r t {rea)K l ( rea ) - sj t ( m)K' t (m ) ’ 


(63) 


where I t (z 0 ) - dl t (z)/d(z) \ and similarly for K t (z 0 ) . Now the potential G (1, (r|r , ) at an 

observation point exterior to the dielectric cylinder due to a point charge also outside of the cylinder is 
expressed as 
G (1, (r|r') = 

X { Jo Gq 1 ) (P’ P'i K ) cos [ a:(z - z')\die +2^ cos[/(^ - (/>')) J q G, (1) (p, p' , re) cos[/e(z - z')\dre 

The derived Green's function, dependent on the permittivity and radius of the dielectric rod under test, is 
used later to calculate the capacitance of the arc -electrode sensor. 
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Note on the choice of the Bessel function kernel 

Instead of using the identity in (55), one can also express 1 / r in terms of the of Bessel function of the 
first kind of order zero J 0 ( xp) [ 121 ]: 


1 f 00 

— = J 0 (xp)e KZ dK, (65) 

y JO 

and express the Green's function Gj 1 ( r | r') in the form of Bessel functions of the first and the second 
kind 


G ( J ] (r|r') = -^xj r K ( {) r> {p,p\ x)e~^ dx +2]T cos[t {</> - K , (1) (p,p', x)e~^ dx | 

4# t (=i I 


( 66 ) 


where 


Kf ] {p, p',K) = J t ( xp)J t {Kp') + Aj{x)Y t (i xp)J t {xp') / = 0,1,---, 


a A k )=- 


(s 2 - s l )J i {xa)j' t {xa) 
e 2 J' t {m)Y t {xa) - e x J t (xa)Y'(xa) 


(67) 

( 68 ) 


J t (xa) = dJfxp )/ d{xp) | a and Y t ( xa ) = dYfxp)/ d{xp) \ . However, the denominator in (68) 


contains an infinite number of zeros for x from 0 to oo , and increases the complexity in MoM numerical 
implementations. Therefore, the Green's function in the form of modified Bessel functions, (64), is a 
better choice here for calculating the sensor output capacitance. 


Numerical implementation 

The capacitance C between the two arc -electrodes is calculated numerically as follows. The Green's 
function derived above is used to set up the integral equation in MoM calculations, which leads to the 
solution for the surface charge density on each electrode. The two electrodes are oppositely charged in the 
numerical calculations. Because of the axisymmetry of the problem, it is only necessary to calculate the 
surface charge density on one of the electrodes. The output capacitance C is then calculated from 

c = 2, (69) 

v 

where the total charge Q on each electrode is obtained by integrating the surface charge density over the 
electrode surface and V is the potential difference between the electrodes. 

Calculation method 

Figure 69 shows the discretization of the arc -electrode surfaces into MxN elements of assumed 
constant surface charge density. Each electrode is discretized into M elements in the (j) direction and N 


97 


elements in the z direction. Denote the surface charge density on the left electrode as <7 s ((f>',z') and that 

on the right electrode as <J s ((f)' + . The potential at the observation point I" = (p {] ,(j), z) on the 

electrode surface due to the charged arc-electrodes can be expressed by integrating (64) over the electrode 


surfaces: 


^{(f), z ) = — JJ G (l) (r | r’)(j s ((fj , z’)p n d(p'dz' 

^ Left electrode 

- — || G (1) (r | r')<r s (^' + n, z')p 0 d(p'dz' 


(70) 


'0 Right electrode 

In the MoM calculations, the following expansion is used to approximate the continuous function 




MN 


°S$\ z ')=Y. (7 , b M'' z '^ 

j = 1 


(71) 


where <T is the unknown constant surface charge density on element j and b i(j)' , z') is the pulse basis 


function 


bU',z'): 


\ 1 on element j 

0 elsewhere. 


(72) 


To solve for the MN unknown coefficients o ■ , weighting (or testing) functions w.(0, z) are 

introduced to force that the boundary condition for the potential in (70) is satisfied for each element on 
the sensor surface. The point-matching method is used, in which the weighting functions are Dirac delta 
functions: 

w i {(/),z)-8{(j)-(j) i )8{z-z i ) on element/, (73) 

where i = 1,2 . Discretizing the integral equation using weighting functions in each of the MN 
elements, (70) is expressed as the following matrix equation: 
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Figure 69 Discretization of the arc -electrode surfaces into M xN elements of assumed constant surface charge density. 
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( 74 ) 


(75) 


All the elements in V share the same potential v that is the potential applied to one of the electrodes. The 
other electrode has potential -v. From (75) the surface charge density <7 s {(j)' ,z) on one of the electrodes 
is solved, and that for the other electrode is simply -<7 s (^ + 7T ,z') . The total charge Q on each electrode 

can be found by integrating <7 s (<f)',z') over the electrode surface. The sensor output capacitance C is 
ultimately calculated through (69). 


Example calculations 

When numerically calculating the matrix element given in (75), the zero to infinity summation and 
integral in G (1, (r | F ) (see (64)) need to be truncated. The convergence of the Green's function depends 
on values of 82 / 81 , alp and <j> Q . When these values are large, large truncation ranges for the summation 
and integral in (64) are needed. It is found that, for the case 82 / 81 = 5 , a / p = 1 , $, = 177 and 1 = 4 cm, if 

one truncates the summation in G (1) (r | G) with 40 terms and the integral with the range from 0 to 2000 
for the off-diagonal components in (75), and the summation with 300 terms and the integral with the 
range from 0 to 2000 for the diagonal components, accuracy to three significant figures can be achieved 
in the final calculated sensor output capacitance C. The cases calculated in Figure 70 to Figure 72 and 

Figure 74 have smaller e 2 /si, alp and <f > 0 values than those in the case calculated above. The truncation 

standard used here is adopted in all numerical calculations of sensor capacitance value in this section. It 
guarantees achieving convergence with accuracy to the third significant digit in all the cases discussed 
below. 

The dependence of sensor output capacitance on the electrode configuration is investigated as follows. 
In Figure 70, sensor output capacitance C is plotted as a function of the electrode length 1 and the arc- 
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angle . In this example calculation, the infinitely long dielectric rod is assumed to be in free space, with 

relative permittivity s r =2.5 and radius a = 9.525 mm (chosen to be similar to the radii of the rods used for 
experiments described later.) The arc-electrodes share the same radius as the cylindrical rod. It is seen 

from Figure 70 that for any fixed electrode arc-angle , there exists a linear relationship between the 
sensor capacitance C and the electrode length I. On the other hand, for any given electrode length /, the 
sensor output capacitance C increases as the electrode arc -angle </> 0 increases, and tends to infinity as </> 0 

tends to 180°. This is explained by the fact that the output capacitance C results from interaction between 
the sensor electrodes. The charge density on the electrodes is highest at the electrode edges, and increases 

as the electrode edges come closer together. As <f) 0 tends to 180°, the gaps between the edges of the two 

electrodes become infinitesimally small and therefore the resulting capacitance tends to infinity, in 
accordance with the singular behavior of the charge density at the electrode edges. Figure 70 shows that 
in order to achieve maximum sensor output signal, the ideal sensor electrodes would be as long as 

practically possible and with large arc -angle <f> 0 . 
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Figure 70 Calculated sensor output capacitance as a function of electrode length I and arc -angle (j> 0 . The dielectric rod is in free 

space, with a relative permittivity of 2.5 and a radius of 9.525 mm. 


Figure 71 shows an example of the sensor output capacitance C as a function of the ratio Cl / p 0 (see 
Figure 67). Rod parameters are as for Figure 70. The arc-electrodes each have fixed radius p 0 = 9.525 
mm, arc -angle = 174.44° and length 1=4 cm. In other words, Figure 71 shows the dependence of 
sensor capacitance on the cylindrical test-piece diameter, for a fixed arc -electrode sensor configuration. It 
is seen from Figure 71 that as the ratio O 7 p 0 increases, sensor output capacitance increases dramatically, 
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especially when this ratio tends to 1. This is because as a / p 0 increases, the average permittivity interior 
to the arc -electrodes increases and therefore C increases. On the other hand, the sensor's most sensitive 
area lies in the region close to the gaps between the two electrodes. As al p {] tends to 1, the arc- 
electrodes are more likely to detect increases in the average permittivity surrounding the sensor. This is 
why the sensor output capacitance C changes more rapidly as the ratio a / p 0 approaches 1 . The 

theoretical calculation in Figure 71 demonstrates that, during measurements, unidentified small air gaps 
existing between the arc electrodes and the dielectric rod under test can introduce relatively large 

uncertainty in the measured C, especially as a — > /7 0 . Therefore, in order to achieve the strongest sensor 

output signal and the smallest uncertainty due to possible air gaps between the electrodes and test-piece, it 
is desirable to have the arc -electrodes in tight surface contact with the test-piece. 



Figure 71 Calculated sensor output capacitance C as a function of the ratio of dielectric rod radius a to electrode radius /7 q . The 
electrode radius, arc-angle, and length are /7 0 = 9.525 nun, (j) = 174.44°, and 1 = 4 mm, respectively. 


The sensor output capacitance C as a function of dielectric rod relative permittivity e r2 is plotted in 
Figure 72, in which different sensor configurations are considered. A linear relationship between the 
sensor output capacitance and the test-piece permittivity is observed and has been verified numerically, 
by computation of a sufficient number of data points (seven in this case). It is seen that the slope of sensor 

output capacitance versus rod permittivity depends on both the sensor configuration and the ratio a / /?„ . 

For a given a / p 0 , the value of the slope increases as the electrode length / and arc-angle <j> {] increase. 
This is because the value of the slope represents changes in the absolute values of the capacitance for any 
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rod permittivity increment. These absolute value changes in capacitance are most obvious for sensors 
with large electrode length 1 and arc-angle (/> 0 values. This also explains why the value of the slope, for 

fixed 1 and <f> 0 values, increases as Cl / p {) increases. However, it is worth pointing out that for fixed 

electrode radius p , arc-angle <j > 0 , and a / p {] values, although increasing electrode length / increases the 

value of the slope, relative changes in capacitance as s r 2 changes stay the same, because of the linear 
relationship between the sensor output capacitance C and electrode length / (see Figure 70). 



Figure 72 Calculated sensor output capacitance as a function of dielectric rod relative permittivity. Cl / /7q — 1 , / = 4 cm and 
= 174.44° except where indicated. All the sensor electrodes have fixed radius /? 0 =9.525 mm. 

Experimental verification 

Capacitance experiments were performed to verify the validity of the developed theory. Two sets of 
rectangular planar electrodes (shown in Figure 73) were fabricated using photolithography by American 
Standard Circuits Inc.. The sensor shape was achieved by selectively etching a 18 p m thick copper 
cladding (14 mL standard) off a flat 25.4 // m thick Kapton ’ type 100 CR polyimide film. These flexible 
electrodes were fixed onto different cylindrical dielectric test-pieces later to form the arc -electrode 
capacitance sensors. The sensor dimensions are vv=29 mm and 1=20 mm for one set and w= 29 mm and 
/=40 mm for the other (see Figure 67). A Nikon EPIPHOT 200 microscope was used to independently 
measure the fabricated sensor dimensions, for the purpose of checking the difference between the 
fabricated dimensions and the specified ones, and therefore the accuracy of the fabrication process. The 
“traveling microscope” measurement method, with accuracy of 0.01 mm, was used to measure the 
relatively large sensor electrode dimensions. It was found that the measured dimensions of the fabricated 
electrodes are identical with the nominal values under such measurement accuracy. 
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Three 304.8-mm-long dielectric rods are used in the measurements to simulate the infinitely long 
cylindrical dielectric rod. The dielectric rods are long compared with the electrode lengths (factors of 
approximately 8 and 1 5 longer), and the edge effect due to finite rod length can be neglected if the sensor 
electrodes are placed at the center of the rods. The rod materials are Acetal Copolymer (Tecaform™), 
Cast Acrylic, and Virgin Electrical Grade Teflon® PTFE. A digital caliper, with accuracy of ± 0.01 mm, 
was used to independently measure the diameter of each rod. The permittivity of each rod was 
independently determined by cutting a slice from the end of each rod, and then measuring the permittivity 
of each slice using a Novocontrol Alpha Dielectric Spectrometer at 1 Mffz. In the Novocontrol 
measurements, both sides of each slice were brushed with silver paint to form the measuring electrodes. 



Figure 73 Agilent E4980A precision LCR meter and Agilent probe test fixture 16095 A used for sensor capacitance 
measurements. Subfigure: photography of the flexible rectangular planar electrodes fabricated using photolithography. 

The rectangular planar electrodes were attached to each dielectric rod by taping the thin Kapton ' sensor 
substrate tightly against the rod material, as shown in Figure 73. The electrodes were aligned carefully so 
that the upper and lower edges of the two electrodes were at the same height, the vertical edges of both 
electrodes were in parallel, and the two vertical gaps between the two electrodes were of the same size, as 
assumed in the theoretical model. Another layer of 25.4-//m-thick Kapton" fdm was wrapped tightly onto 
the outsides of the electrodes in order to minimize the air gap between the electrodes and the dielectric 
rod. Because the Kapton® films used were quite thin, influences from their permittivity on the 
measurement signal were negligible. 

For each dielectric test-piece used in the benchmark experiments, the test-piece material, test-piece 
diameter, independently measured test-piece relative permittivity, electrode radius p {] , and electrode arc- 

angle ()>q are listed in Table 18. The electrode radius p {) for each rod is obtained by summing the 
dielectric rod radius and the Kapton® substrate thickness. The electrode arc -angle (/> 0 in Table 18 is 
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different for each test-piece because of the fact that the diameter of each rod is different while the 
electrode width w is the same. The parameters shown in the table were used as the inputs in the numerical 
calculations. 

An Agilent E4980A precision LCR meter was used to measure the sensor output capacitance. The LCR 
meter operating frequency was set as 1 MHz, so that the measurement error from the LCR meter was less 
than 0.3% for a 1 pF capacitance while at the same time being a good approximation for the electrostatic 
assumption in the numerical model (results of the calculation do not depend on frequency). A static model 
can be applied for this configuration even at/= 1 MHz because, at this frequency, the corresponding 
wavelength A = 300 m. The diagonal dimensions of the electrodes in the measurement are smaller than 5 
cm, which means the maximum phase change over the electrode surfaces is less than 0.06° and the effect 
of scattering is therefore negligible. A 0.06° phase change is not detectable here, being below the 
measurement sensitivity of the Novocontrol dielectric spectrometer and the LCR meter. If a lower 
operating frequency is desired for practical capacitance measurements, an impedance measurement 
instrument with higher accuracy when measuring large impedance values should be used. (According to 
the relationship Z = 1 / j2n fC , the impedance Z resulting from measuring a given capacitance C under a 
lower frequency / will be larger). 


Table 18 Parameters of the dielectric test-pieces and the arc-electrode sensors used in benchmark experiments. The areas of the 
two sets of sensor electrodes are 29 X 20 mm 2 and 29 X 40 mm 2 , respectively. 


Test-piece 

Test-piece 

Measured test-piece 

Electrode radius 

Electrode arc- 

material 

diameter (mm) 

permittivity 

/>D (mm) 

angle (j ) o 

Tecaform™ 

19.08 ±0.01 

3.77 ±0.05 

9.565 ±0.005 

173.71° ± 0.10° 

Cast Acrylic 

19.03 ±0.01 

2.88 ±0.05 

9.540 ±0.005 

174.17° ±0.09° 

Teflon® 

19.10 ± 0.01 

2.23 ±0.05 

9575 ±0.005 

173.53° ±0.09° 


Sensor output capacitance C was measured by placing the probe of an Agilent probe test fixture 
16095A across the two sensor electrodes, as shown in Figure 73. The parts on the electrodes where the 
probe is in surface contact were not covered by Kapton ’ films. This probe test fixture was connected to 
the LCR meter and the measured capacitance was read from the LCR meter screen. Figure 74 shows the 
comparison between the calculated and measured sensor output capacitance for each test-piece material 
and the two different electrode configurations. Experimental data show excellent agreement with 
numerical results (to within 3%), and the maximum absolute difference in capacitance is less than 0.1 pF. 
It is worth pointing out that even if Kapton ” films are tightly wrapped around the electrodes to attempt to 
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eliminate the air gaps between electrodes and the test-piece, small gaps still exist. In particular, the 
vertical edges of the electrodes tend to bend up, giving rise to small air gaps, where the sensor is most 
sensitive. This points to the fact that the ideal way to achieve best agreement between theory and 
measurements is to deposit the arc -electrodes directly onto the cylindrical test-pieces. Thus errors coming 
from the misalignment of sensor electrodes and the existence of air gaps will be eliminated. On the other 
hand, deposition of electrodes directly onto the test-piece is costly, time-consuming and undesirable for 
most practical purposes. 

One purpose for these arc-electrode sensors is the inverse determination of the permittivity of materials 
under test from measured sensor capacitance. Comparisons are made in Table 19 to assess the capability 
and accuracy of the arc-electrode sensors in material dielectric property characterization. The sample 
permittivities inferred from measured capacitance values shown in Figure 74 are compared with those 



Electrode length I (cm) 

Figure 74 Measured and calculated C for various sensor configurations (see Table 18) in contact with different dielectric test- 
pieces. Measurement results and error bars are denoted by the black symbol. 


Table 19 Comparison of test-piece permittivity values between independently measured ones and inversely determined ones 

from measured capacitance using the arc-electrode sensors. 


Test-piece 

material 

Independently 

measured 

permittivity 

Electrode 
length (mm) 

Inversely 

determined 

permittivity 

Relative 

difference 

(%) 

Tecaform™ 

3.77 ±0.05 

20 

3.76 ±0.07 

0.3 

40 

3.82 ±0.08 

1.3 

Cast Acrylic 

2.88 ±0.05 

20 

2.88 ±0.01 

0.0 

40 

2.89 ±0.01 

0.3 

Teflon® 

2.23 ±0.05 

20 

2.16 ± 0.11 

3.1 

40 

2.18 ±0.09 

2.2 
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measured by a Novocontrol dielectric spectrometer, in the manner discussed in the second paragraph in 
Section 1.4. Again, excellent agreement between inferred and independently measured test-piece 
permittivities is obtained. Table 19 demonstrates the great potential of using the arc-electrode capacitive 
sensors for accurate and convenient permittivity measurements of cylindrical dielectric materials. 

Summary 

A numerical model, based on the electrostatic Green's function due to a point source exterior to an 
infinitely long cylindrical dielectric, has been developed to quantitatively evaluate the dielectric property 
of cylindrical dielectric test-pieces. The quantitative dependence of the sensor output capacitance on the 
test-piece permittivity and radius has been demonstrated numerically and verified experimentally. The 
permittivity of various cylindrical test-pieces has been inferred from measured capacitance to within 1% 
accuracy, on average. 

2. Analysis of a capacitive sensor for the evaluation of circular cylinders with a conductive 
core 

A capacitive sensor has been developed for the purpose of measuring the permittivity of a cylindrical 
dielectric that coats a conductive core cylinder. The capacitive sensor consists of two identical curved 
patch electrodes that are exterior to and coaxial with the cylindrical test-piece. The permittivity of the 
cylinder is determined from measurements of capacitance by means of a physics-based model. In the 
model, an electroquasistatic Green's function due to a point source exterior to a dielectric-coated 
conductor is derived, in which the permittivity of the dielectric material may take complex values. The 
Green's function is then used to set up integral equations that relate the unknown sensor surface charge 
density to the imposed potentials on the electrode surfaces. The method of moments is utilized to 
discretize the integral equation into a matrix equation that is solved for the sensor surface charge density 
and eventually the sensor output capacitance. This model enables the complex permittivity of the 
dielectric coating material, or the geometry of the cylindrical test-piece, to be inferred from the measured 
sensor capacitance and dissipation factor. Experimental validation of the numerical model has been 
performed on three different cylindrical test-pieces for two different electrode configurations. Each of the 
test-pieces has the structure of a dielectric coated brass rod. Good agreement between measured and 
calculated sensor capacitance (to an average of 7.4%) and dissipation factor (to within 0.002) was 
observed. Main sources of uncertainty in the measurement include variations in the test-piece geometry, 
misalignment of sensor electrodes, strain-induced variation in the test-piece permittivity, and the 
existence of unintended air gaps between electrodes and the test-piece. To demonstrate the effectiveness 
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of the sensor, measurements of capacitance have been made on aircraft wires and the permittivity of the 
insulation inferred. A significant change in permittivity was observed for thermally degraded wires. 

Introduction 

This part of the report describes development and benchmark testing of a model-based capacitive 
method for complex permittivity measurement of a cylindrical dielectric that coats a conductive core 
cylinder. The work is motivated by the effective evaluation of degradation status of air- and space-craft 
wiring insulation. Degradation in electrical wiring insulation has the potential to cause aviation 
catastrophe due to consequent short-circuiting or loss of control function [88], One effective approach of 
evaluating insulation degradation state is through permittivity measurements, which can be achieved 
using capacitive methods. 

Capacitive methods offer a favorable solution to the accurate characterization of material dielectric 
properties at low costs. For example, model-based interdigital sensors allow the inverse determination of 
test-piece material properties from measured sensor output capacitance [109]. Applications of interdigital 
sensors include humidity and moisture sensing, electrical insulation properties sensing, chemical sensing, 
biosensing, and others. Rectangular coplanar capacitive sensors have been developed to detect water 
intrusion in composite materials, on the physical basis that changes in material dielectric properties lead 
to variations in the sensor capacitance [117]. Rectangular capacitance sensors also find application in 
damage detection in laminated composite plates [118], evaluation of moisture content in reinforced 
concrete covers [122], and rain sensing [123], Rectangular capacitive sensor arrays have been reported in 
[119], and shown to be capable of detecting surface and subsurface features of dielectrics. Circular shaped 
capacitive sensors have also been developed for the quantitative characterization of material properties. 
Multichannel fringing electric field sensors [114], cylindrical geometry electroquasistatic dielectrometry 
sensors [115] and concentric coplanar capacitive sensors [116] are some examples of these. In addition, 
parallel plate capacitors formed by two circular discs with coplanar guard electrodes have been used to 
detect strength-limiting defects (large voids) in cellular glasses [124], Capacitive techniques that have the 
potential to be integrated into aircraft wiring test systems are presented in [125] to [99]. Arc -electrode 
capacitive sensors have been developed to characterize material and structural properties of cylindrical 
dielectric rods [125]. A so-called ‘meander' coil and a quarter-circular interdigital sensor have been used 
in detection of cable insulation damage [113]. In [99], linear relationships between the capacitance of 
open-circuited aircraft wires and their length have been demonstrated and enable the determination of 
wiring conductor length from measured capacitance values. Further references on capacitive methods can 
be found in [116] and [125], 
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In addition to capacitive techniques, other electrical testing methods have been developed to evaluate 
the wiring conductor condition. Reflectometry is one of the most commonly used techniques for aircraft 
wiring testing, in which a high frequency electrical signal is transmitted along the wire and any 
impedance discontinuities result in reflected signals whose interpretation may give an estimate of the flaw 
position. An excellent review paper that compares different reflectometry methods is [91]. Low-voltage 
resistance tests and dielectric-withstand-voltage tests are two qualitative methods that can detect faults, 
but are not suitable for inspection of aging aircraft wiring because of the difficulty of miniaturization and 
pinpointing the fault [90], 

The fundamental relationships between the capacitance of a curved patch capacitive sensor and the 
dielectric property of its test-piece are presented. The test-piece discussed here has the structure of a 
dielectric-coated conductor. The sensor consists of two identical and symmetric curved patch electrodes, 
which are located exterior to and coaxially with the cylindrical test-piece. For the modeling, an 
electroquasistatic Green's Junction due to a point source exterior to a cylindrical dielectric-coated 
conductor is derived in cylindrical coordinates, in which the permittivity of the dielectric material may 
take complex values, Section 2.2. The Green's function is then used to set up integral equations that relate 
the unknown sensor surface charge density to the predetermined electric potentials on the sensor 
electrodes. The method of moments (MoM) is utilized to discretize the integral equation into a matrix 
equation that is solved for the sensor surface charge density, Section 2.3. The total charge Q on each 
electrode is then calculated and the capacitance between the two curved patch electrodes is obtained 
through C = Q/V, in which V is the potential difference between the two electrodes. Large-scale 
experimental verification of the theory has been carried out on copper rods coated with different dielectric 
materials, for two different electrode configurations, Section 2.4. Measured and numerically calculated 
sensor capacitance values agree to an average of 7.4%, whereas the dissipation factors agree to within 
0.002. Major factors contributing to the measurement uncertainty are variations in the test-piece geometry, 
misalignment of sensor electrodes, strain induced variation in the test-piece permittivity, and the existence 
of unintended air gaps between the sensor electrodes and cylindrical test-pieces. These are discussed in 
detail in Section 2.4. 

This part of Chapter V is focused on the development and verification of a physics-based measurement 
method. In related research that will be published in a later article, a prototype capacitive probe has been 
built based on the method presented in this section [126]. This probe has been applied for quantitative 
characterization of insulation degradation on actual aircraft wires. Changes in the insulation complex 
permittivity, induced by thermal and hydrolytic exposures, have been clearly detected. One of these 
results is shown here, Section 2.5, in order to demonstrate the feasibility of quantitative evaluation of 
wiring insulation permittivity using the model described in this section. The capacitance technique 
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developed in this section has the potential to be built into smart embedded wiring test systems of the 
future. It is complementary to large-scale inspection techniques. For instance, two ultrasonic transducers 
can be used in a pitch/catch configuration to generate and receive an ultrasonic guided wave in a wire, 
and obtain an overall indication of the wire insulation condition [101]. When faults are indicated by such 
system-level inspection techniques, high accuracy local inspections using the curved patch capacitive 
sensors may follow. 

Modeling 

The method of moments (MoM) is utilized in the modeling process instead of the finite element method. 
The adopted numerical method has the advantage of only needing to discretize surfaces, rather than the 
volume, to obtain the sensor capacitance. The number of unknowns to solve is proportional to N 2 , where 
N is the number of unknowns in one dimension. Most commercially available finite element tools, 
however, have to discretize a truncated space, in which the number of u nk nowns to solve is proportional 
to N\ This means using the MoM in the modeling significantly reduces the number of unknowns to solve 
and therefore the required computation time and computer memory. This feature is important in 
quantitative evaluation of wiring insulation permittivity using the capacitive method. Because the 
insulation permittivity is determined inversely based on measured capacitance through the numerical 
model, faster numerical models can significantly reduce the amount of time needed during this process. 

Sensor configuration 

The configuration of the curved patch capacitive sensor is shown in Figure 75. The two identical and 
symmetric curved sensor electrodes are assumed in the theory to be infinitesimally thin. The cylindrical 
test-piece is modeled as an infinitely long dielectric-coated perfect conductor. The conductor is kept at 
ground potential in the modeling and the patch electrodes are held at equal and opposite potentials, ± V/2. 
The sensor output capacitance is calculated in the electroquasistatic regime, in which the permittivity of 
the dielectric coating may take complex values while the resulting potential still satisfies the Laplace 
equation. 


109 



Sensor electrodes 


Figure 75 Curved patch capacitive sensor. The radii of the sensor electrodes, the conductor, and the cylindrical test-piece are 
denoted /? 0 , a, and b, respectively. The arc -angle of each sensor electrode is (j>Q (rad). The length of each electrode in the 

vertical direction is l and the width in the horizontal direction is W = (j)^ X p . 


Derivation of Green's function in cylindrical coordinates 

Figure 76 shows a point source placed at ( p',(j)',z ') exterior to a cylindrical dielectric-coated 
conductor. This configuration is used in the following Green's function derivation. The outer radius of the 
cylindrical test-piece is b and the radius of the inner conductor is a. The complex permittivities of the 
dielectric coating and the background medium are s 2 *= e 2 ’ -ys 2 ” and Si*= sf -je i”, respectively, where 

j - \[—\ . Considering the canonical shapes of the sensor electrodes and the test-piece, the 
electroquasistatic Green's function is derived in cylindrical coordinates. The electric potential at an 
observation point (p,(/),z) due to the point source in Figure 76 satisfies the Laplace equation in each 
homogeneous region exterior to the conductor: 


1 d 


r 


\ 


P- 


+ - 


1 d 2 


p dp ^ dp ) p l d(/) 2 dz 2 


T (i) (p, </,, z) = — —S(p — p’)S(</> - f)S{z - z’) i = 1, 2, 

P 


(76) 


and is subject to a Dirichlet boundary condition at the surface defined by p = a: 


x ¥ (2) {p = a,(j),z) = 0, (77) 

and the interface conditions at the surface defined by p = b : 

x V m {p = b,<f>,z) = '¥ (2 \p = b,<f>,z ) 

£ f^ (1) (A </>, z ) I p=b=s -^-'F (2) (/7, </>, z ) | b 
dp dp 


( 78 ) 
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where superscripts (1) and ( 2 ) correspond to the regions defined by p > b and a < p <b , respectively. 
Solution to the potential Green's function can be found by following similar procedures in Section 2.2 
Chapter V, in which the Green's function due to a point source exterior to a homogeneous dielectric 
cylinder is derived. Steps taken to derive the electroquasistatic Green's function are summarized as 
follows. 

Begin with the free space Green's function in cylindrical coordinates due to a point charge at the 
origin. The potential at an observation point ( p,(f),z ) is [120] 

1 If 00 

G 0 (p,z) = - — = - 2-1 K 0 (Kp)cos(tcz)dK, ( 79 ) 

4 nr in Jo 

where r = yjp 2 + z 2 and K t) (kp) is the modified Bessel function of the second kind of order zero. The 

potential due to a source point at arbitrary location ( p' ,(f)' ,z ') can be obtained by applying the addition 
theorem given in ( 80 ) [ 120 ] 

OO 

(< ) = h (* 7 2 < )K 0 (Kp y ) + 2 ^ cos[m(</> - f)]I m (/cp < )K m (k’p > ), ( 80 ) 

m = 1 

to K 0 (tcp ) in (79) and replacing z by z-z’, where C, = p 2 + p' 1 -2 pp' cos(0 - (/>') , p < is the lesser of 
p and p' , and p y is the greater. Now the Green's function in free space is rewritten as 


] ( poo ~ 

G 0 (p, <j), z | p\ (/)', z) = — T x G 0 (p, p\ k) cos[x-(z - z')]dK 

2 n r° 

+2 2 cos[/??(^ - ^')]| (| G m (p, p\ k) cos[x-(z - z')]dK L 
m = 1 J 


( 81 ) 
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where G p (p,p',K) = I p (x-p < )K p (>cp > ) p = 0,1,2, •••, and I p {Kp K ) is the modified Bessel function 

of the first kind. As pointed out in [125], the difference between the Green's function due to a point 
charge in free space and that due to a point charge exterior to an infinitely long dielectric-coated 
conductor Figure 76 is due to interface conditions that are applied at surfaces of constant p , 
corresponding to the physical interfaces of the test-piece. To find the Green's function for this case, one 
needs only to modify the integral kernel G (p, p' , k) in (81) so that the interface conditions in (77) and 
(78) are satisfied: 


G ( p \p = a,p',/c) = 0 


G^(p,p',K) = G ( ;\p,p',K) 

£ i^-G ( p ) {p,p\k) = £ 2 ^-G ( p ) (p,p\k:) 

dp dp 


(82) 

(83) 

(84) 


Following the same steps as in [125], the Green's function G{p,(f),z | p',<f)',z') at an observation point 


( p,(f),z ) due to a point charge at ( p', (j)\ z) that is also exterior to the dielectric-coated conductor is 
derived as: 


where 


1 f°° ~ 

<j( r I r ') = ^G 0 {p,p',K)cos[K{z-z'^fiK 


0 l 0 


(85) 
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( 86 ) 


(87) 


a p (k) - I p (kq) / K p (jca), I p (jcb) = dl P (icp) / d(icp)\ _ and similarly for K p (tcb) . The Green's 


function in the form of modified Bessel functions, (85), is used in the following calculations of the sensor 
output capacitance. The sensor capacitance is computed later using the derived test-piece geometry and 
permittivity dependent Green's function. 

Note that (85) can be simplified to the case of a homogeneous dielectric rod, described in [125], by 
assigning a = 0 . The Green's function (85) can also be expressed in the form of Bessel functions of the 
first and second kind. Flowever, the denominator of the integrand in the Green's function contains an 
infinite number of zeros and increases the complexity in the numerical implementations. 
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Numerical implementation 


Calculation method 

The sensor output capacitance C is calculated numerically using the method of moments (MoM). 
Calculation procedures used here are similar to those in [125]. In summary, the following steps are taken 
to compute the sensor capacitance. 

First, the Green's function (85) is used to set up the integral equation that relates the u nk nown sensor 
surface charge density cr s (0',z') to the imposed potential 'F(^,z) on the sensor electrodes 


v F(^,z) = [[ G(p, tj), z | p' , tj)' , z')cr (</>', z')p 0 dtj)'dz’ 

J J Left electrode /oo\ 

G{p,<j),z\p' ,tj)' ,z')g^' ^n^p^dtp'dz’. 

Right electrode 

In order to solve for the sensor surface charge density numerically, i.e., to use discrete functions 
approximating the continuous function cr s ((f)',z') , each electrode in Figure 77 is discretized into MxN 

rectangular elements. The charge density on each element is assumed to be constant and can be different 
from others. Mathematically, this approximation is expressed as 

MN 

(89) 

H 


where b .((/>' ,z') is the selected pulse basis function and <J . is the unknown coefficient to be determined. 

Note the axisymmetry of the problem, it is only necessary to calculate the surface charge density on one 
of the electrodes. 

To solve for the MN unknown coefficients cr j , weighting functions H’ ( (/), z) are introduced to force 

the integral equation (88) to be satisfied for each element on the sensor surface. The point-matching 
method is used in this process, in which the weighting functions are Dirac delta functions. Expressions for 
b j {(j)' , z') and w^tj), z) can be found in [125], Discretizing the integral equation using weighting 


functions in each of the MV elements, (88) is expressed as the following matrix equation: 
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where 


Gj = { j element . G(p i , (j), , z, | p'j , (ft) , z' )bj (f, z')p 0 d<j>'dz', 


(91) 
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L = MN and all the elements in V share the same potential that is applied to the electrode. The 
unknown coefficients cr are obtained by solving the matrix equation, and the total charge Q on each 

electrode is calculated by integrating the surface charge density over the electrode surface. The 
capacitance C between the two electrodes is obtained using the relationship C=Q/V. 



Figure 77 Curved patch capacitive sensor is divided into MxN elements on each electrode, each with assumed constant surface 

charge density. 


Example calculations 

The dielectric-coated conductor is assumed to be in free space in the following calculations, i.e., s r i*=l. 

When numerically evaluating elements in the G matrix, the zero to infinity summation and integral in 
G{p,(j),z | p',(/)',z') (see (85)) need to be truncated. The convergence of the Green's function depends on 
values of 82, tanS = 82/si, alb, b / p and (f) () . When these values are large, large truncation ranges for the 
summation and integral in (85) are needed. It is found that, for the case s r 2 ,= 5, tanS = 0.02 , a/b = 0.9, 
b / p = 1 , </> 0 - 178 and 1 = 4 cm, if one truncates the summation in G(p, (j), z \ p' , (j)' , z ) to 80 terms and 

the integral with the range from 0 to 6000 for the off-diagonal components in G matrix, and the 
summation to 400 terms and the integral with the range from 0 to 6000 for the diagonal components, 
accuracy to three significant figures can be achieved in the final calculated complex sensor output 
capacitance C, for both real and imaginary parts. The cases calculated in Figure 78 to Figure 80 and 

Section 2.4 have smaller e r2 ’, tan 8 , a/b, b / p and <j> 0 values than those in the case calculated above. The 

truncation standard used here is adopted in all of the following numerical calculations, which guarantees 
convergence to three significant figure accuracy is achieved for all the numerically calculated capacitance 
values in this section. 

The dependence of sensor capacitance on the electrode configuration is shown in Figure 78. The sensor 
output capacitance is plotted as a function of the electrode length 1 and the arc-angle <f > 0 . The dielectric 
coating has a relative permittivity s r2 ’=2.5. The radius of the core conductor a = 8 mm and the outer 
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radius of the dielectric-coated conductor b = 9 mm. The sensor electrodes are assumed to be right on the 
cylindrical test-piece: p 0 — b . In Figure 78, a linear relationship between the sensor capacitance C and the 

electrode length 1 is observed for any fixed electrode arc -angle <f) 0 . On the other hand, the sensor output 
capacitance C increases as the electrode arc-angle (/> 0 increases for any given electrode length /, and tends 

to infinity as </> 0 tends to 180°. This is explained by the fact that as <j> 0 tends to 180°, the gaps between the 
two electrodes become infinitesimally small and the resulting capacitance tends to infinity. Figure 78 
shows that the sensor output capacitance changes dramatically when <j> 0 and / have large values. When 
performing dielectric measurements in practice, it is usually helpful to have large sensor output signal and 
therefore to have large </> 0 and / values. Flowever, when 1 and <f> 0 are large, C changes rapidly, and it is 

important to have accurate sensor configuration information to infer accurately material dielectric 
properties from measured C. 



Figure 78 Calculated sensor output capacitance as a function of electrode length I and arc -angle (f) 0 . The rod is in free space 
with conductor radius Q — 8 mm, dielectric radius b — p Q = 9 mm and dielectric permittivity s r2 ’=2.5. 

The dependence of sensor capacitance on the test-piece geometry is shown in Figure 79. The sensor 
output capacitance C is plotted as a function of the ratio b / p 0 and the ratio a/b (see Figure 75). The 

dielectric coating permittivity s r2 ’ is as for Figure 78. The sensor electrodes have fixed radius p {) = 9 mm, 
arc -angle <f> 0 = 174° and length 1 = 4 cm. It is observed that for fixed b / p () values, the sensor capacitance 

increases as the ratio a/b increases. Such a trend is more obvious when the ratio b / p 0 tends to 1. This is 

because the calculated capacitance C is actually the series capacitance of the capacitance between one 
electrode and the core conductor and the capacitance between the core conductor and the other electrode. 
When the ratio a/b increases, the distance between the sensor electrodes and the core conductor decreases 
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and the resulting total capacitance increases. In particular, when b / p 0 =\, the output capacitance tends 
to infinity as a/b tends to 1 , in which case the gaps between the sensor electrodes and the core conductor 
approaches zero. This also explains why the sensor capacitance increases as the ratio b! p {) increases for 
given a/b values, and why such changes in capacitance are more rapid for large a/b values. The fact that 
the overall permittivity of the region between the electrodes and the conductive core increases as b / p 0 

increases also contribute to increases in the sensor output signal. In summary, the existence of the 
conductive core in the test-piece increases the output capacitance for given sensor configurations, and as 
the conductive core radius a approaches zero, the test-piece reduces to a homogeneous dielectric rod. 



Figure 79 Calculated sensor output capacitance C as a function of the ratio of cylindrical test-piece outer radius b to electrode 
radius /7 q and the ratio of conductive core a to cylindrical test-piece outer radius b. The electrode radius, arc -angle and length 

are 9 mm, 174° and 4 cm, respectively. 

Figure 80 shows the sensor capacitance C and dissipation factor D as functions of the dielectric coating 
real permittivity s r2 ’ and imaginary permittivity s r 2 ”, respectively. Different sensor configurations are 
considered. In Figure 80a), a linear relationship between C and s r2 ’ is observed for all sensor 
configurations. It is seen that the slope of sensor capacitance versus dielectric coating real permittivity, 
i.e., the sensor sensitivity, depends on both the sensor configuration and the geometry of the cylindrical 
test-piece. The largest slope in Figure 80 occurs when a/b, b / p {) , electrode length / and electrode arc- 
angle (j) Q are the largest of the values considered. Flowever, it is worth pointing out that although 

increasing electrode length / increases the value of the slope, relative changes in capacitance stay the 
same, because of the linear relationship between the sensor output capacitance C and electrode length / 
(see Figure 78. For practical inspection of wires, a/b is fixed, and it is therefore important to keep b / p (] 

close to 1 to achieve the highest sensitivity. In the selection of (f) 0 , a trade-off exists. Larger (j) {) gives rise 
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to larger sensitivity as well as capacitance. On the other hand, larger (f) {] means that the inter-electrode 

gap decreases and the penetration of the field into the insulation decreases as a consequence. 

Similar relationships between D and s r2 ” are observed in Figure 80b). The major difference between the 
response of C and D to the investigated factors is that D is less sensitive than C to changes in /, (j) {] and 
a/b. 



Figure 80 Calculated sensor capacitance and dissipation factor as a function of the dielectric coating real permittivity s r 2 ’and 
imaginary permittivity s r 2 ”. Sensor configuration: /? 0 = 9 mm, alb — 0.8 , b / /? 0 = 1 , / = 4 cm and (j)^ — 170 except 
where indicated in other lines. s r2 ” and the material dissipation factor are assumed to be zero in a). s r 2 ,= 2 in b). 

Dependence of capacitance on test-piece permittivity and sensor configuration 
If / — > oo , a — > 0 and 6/ p 0 - 1, the case of Figure 75 becomes a two-dimensional problem. An 

analytical expression for the capacitance per unit length between the two curved patches has been derived 
in [127] and takes the following general form 

C = F(f+£ 2 ), (92) 

where F s is a shape factor that depends solely on the capacitor geometry. Considering the linear plots of 

Figure 80a), the following relationship is found to hold in general for the problem discussed in this 
section (Figure 77): 

C = F s (s l + a£ 2 ), (93) 

where a is a constant showing the dependence of C on Si and e 2 . a >\ means C is more dependent on e 2 
than Ei, and vice versa, a = 1 means C depends equally on 8i and s 2 . The factors a and F s in (93) may 
be determined by selecting two data points on any of the lines in Figure 80a). It is found that a and F s 

are constant for any given sensor configuration, independent of the particular data points selected for the 
calculation. 
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As an example shown in Figure 80a), C obtained based on (93) fits nicely on the dashed line for 
a = 2.61 and F = 1.01 m. Similar results are observed for all the other sensor configurations as well. It 

is found that a increases as b / p 0 , alb and / increase, and as (/> () decreases. The latter relationship can 
be explained by the fact that as (j> {] decreases, more electric field penetrates the dielectric coating, and 82 
will therefore have larger impact on C . The product CxF s is the slope in the C versus 82 plot and shows 
the sensor sensitivity, whose dependency on the sensor configuration has been discussed earlier. 

Experiment 

Experimental setup and measurement procedures 

Benchmark experiments comparing the measured sensor capacitance with numerically-predicted values 
were performed to verify the validity of the developed theory. Experiments were conducted at frequency 
1 MHz. Note that although the numerical model is developed in the electroquasistatic regime, i.e., the 
wave length X (approximately 300 m at 1 MHz) is much greater than the dimension of the test-pieces in 
the experiment, the dielectric coatings still have complex permittivities. This is due to losses arising in the 
materials due to the polarization response of the polymers lagging behind the switching of the applied 
electric field at 1 MHz. For this reason, complex permittivities are considered in the following benchmark 
experiments. 

Two sets of rectangular planar electrodes (shown in Figure 81) were fabricated using photolithography 
by selectively etching a 1 8-/./m-thick copper cladding (14 mL standard) off a flat 25.4-/un-thick Kapton 8 
type 100 CR polyimide film. These flexible electrodes were attached to different cylindrical test-pieces 
later to form the capacitance sensors. The sensor dimensions are w = 29 mm and / = 20 mm for one set 
and w = 29 mm and / = 40 mm for the other (see Figure 75). A Nikon EPIPHOT 200 microscope was 
used to measure the fabricated sensor dimensions, for the purpose of checking the difference between the 
fabricated dimensions and the nominal ones, and therefore the accuracy of the fabrication process. The 
“traveling microscope” measurement method, with accuracy of ± 0.01 mm, was used to measure the 
dimensions of the relatively large sensor electrodes. It was found that the measured dimensions of the 
fabricated electrodes are the same as the nominal ones under such measurement accuracy. 

Three cylindrical test-pieces, each being a dielectric tube coated copper rod, were used in the 
measurements to simulated the infinitely long dielectric -coated conductors in theory (see Figure 75). The 
dielectric tube materials were Acetal Copolymer (Tecaform™), Acrylic, and Virgin Electrical Grade 
Teflon" PTFE, respectively. The dielectric tubes were hollowed from homogeneous rods so that the inner 
diameter of the tube matched the diameter of the copper rod as closely as possible. All the three dielectric 
tubes were in tight surface contact with the central copper rods. The cylindrical test-pieces used are 152.4 
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mm in length. These test-pieces are long compared with the electrode lengths (factors of approximately 4 
and 7 longer). The edge effect due to finite rod length can be neglected if the sensor electrodes are placed 
sufficiently far from the ends of the test-piece. Prior to being hollowed out, the permittivity of each 
dielectric tube was measured, by cutting a disc from the end of each rod and measuring its permittivity 
using a Novocontrol Alpha Dielectric Spectrometer at/= 1 MHz and room temperature. The measured 
test-piece permittivity values, together with uncertainty, are provided in Table 20. In the Novocontrol 
measurements, both sides of each disc (around 19 mm in diameter) were brushed with silver paint to form 
the measuring electrodes. The capacitance values resulting from these discs were between 2 and 5 pF. The 
impedance measurement accuracy of the Novocontrol at 1 MHz and room temperature for such 
capacitance values is 0.1% for the magnitude and 0.06° for the absolute phase accuracy. Note since 
Teflon" are low loss materials (loss tangent below 1 0 4 in the frequency range 20 kHz to 1 MHz at room 
temperature [128]), accurate measurements of their imaginary permittivities were not achieved using 
either the Novocontrol or the Agilent E4980A LCR meter with the Agilent 1645 IB dielectric test fixture. 

The diameter of the copper rods, the outer diameter of the dielectric-coated conductors and the 
thickness of the dielectric coatings were measured using a digital caliper with accuracy of ± 0.01 mm. 
Each of these values was measured at ten different locations on the test-piece. An average value and the 
corresponding deviation were obtained, in which the deviation is defined as the maximum absolute 
difference between the average value and the ten measured values. 



Figure 81 Agilent E4980A precision LCR meter and Agilent probe test fixture 16095A used for sensor capacitance 
measurements. Subfigure: photography of the flexible rectangular planar electrodes fabricated using photolithography. 

The experimental arrangement for the capacitance measurements is shown in Figure 81. The 
rectangular patch electrodes were conformed to each cylindrical test-piece by taping the thin Kapton ' 
sensor substrate tightly against the dielectric material. The thickness of the film (25.4 //m) is accounted 
for in the numerical modeling while effects of its permittivity are neglected. This is because the sensor 
capacitance is much more sensitive to small variations in b I p 0 , when this ratio is very close to 1, than 
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those in the test-piece permittivity. The permittivity of the substrate was hence assumed to be that of the 
test-piece, an assumption that introduces negligible uncertainty. The electrodes were aligned carefully in 
order to achieve the sensor configuration in the theoretical model as closely as possible. The goal in the 
alignment was to keep the upper and lower edges of the two electrodes at the same height, the vertical 
edges of both electrodes in parallel, and the two vertical gaps between the two electrodes of the same size. 
Another layer of 25.4-//m-thick Kapton ’ film was wrapped tightly onto the outsides of the electrodes to 
further minimize the air gap between the electrodes and the cylindrical test-piece, leaving part of each 
electrode exposed to make electrical contact with the probe test fixture later. 

In the experimental verifications, two groups of capacitance measurements were performed for each 
cylindrical test-piece: one group using the 20-mm-long electrodes and the other using the 40-mm-long 
electrodes. For each cylindrical test-piece, the tube material, brass rod diameter, outer diameter of the 

dielectric-coated brass, variation in the dielectric tube thickness, electrode radius p 0 and electrode arc- 

angle p {) are provided in Table 21, with uncertainties. Because the two types of electrodes were attached 
at different locations on each test-piece, the outer diameters of the test-pieces in Table 21 were measured 
at those individual locations and vary slightly. The electrode radius p 0 for each cylindrical test-piece is 
obtained by summing the outer radius of the dielectric-coated brass rod and the Kapton" substrate 
thickness. The electrode arc -angle p 0 and its uncertainty in Table 21 are calculated from the electrode 

width w, electrode radius p () and its uncertainty. Each test-piece has its own electrode arc -angle </> 0 

because the fabricated electrodes share a fixed width w while the cylindrical test-pieces have different 
radii. The parameters shown in Table 20 and Table 21 were used as the inputs in the numerical 
calculations. 

An Agilent E4980A precision LCR meter was used to measure the sensor output capacitance at room 
temperature. The ECR meter operating frequency was set at 1 MHz to approximate the electroquasistatic 
assumption in the numerical model. Under these conditions, the measurement accuracy of the LCR meter 
for a 4 pF capacitance is 0.15% and the absolute accuracy for the dissipation factor is 0.0015, whereas 
those for a 13 pF capacitance are 0.13% and 0.0013, respectively. The measured capacitance values in 
this section are all within 4 and 13 pF. If a lower operating frequency is desired in capacitance 
measurements, an impedance measurement instrument with higher accuracy when measuring large 
impedance values should be used. (According to the relationship Z = 1/ j2n fC , the impedance Z 
resulting from measuring a given capacitance C under a lower frequency / will be larger). The sensor 
capacitance C was measured by placing an Agilent probe test fixture 16095 A across the two sensor 
electrodes, as shown in Figure 81. This probe test fixture was connected to the LCR meter and the 
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measured capacitance was read from the LCR meter screen. Note that in the modeling the two electrodes 
are assumed oppositely charged and the conductive core of the test-piece is kept at ground potential. The 
calculated capacitance is the series capacitance of the capacitance between one electrode and the 
conductive core and the capacitance between the conductive core and the other electrode. When 
performing capacitance measurements, however, one needs only to place the probe test fixture across the 
two sensor electrodes, as shown in Figure 81. The potential on the conductive core is the average of the 
potentials on the two sensor electrodes, due to the symmetry of the problem, and the capacitance picked 
up by the probe is the series capacitance calculated in the numerical model. 


Table 20 Measured complex permittivity values of the dielectric coating materials. 


Dielectric tube material 

Measured dielectric 
tube real permittivity £ r2 ’ 

Measured dielectric 

tube imaginary permittivity £ r2 ” 

Tecaform™ 

3.77 ±0.05 

0.014 ±0.002 

Acrylic 

3.02 ±0.05 

0.055 ±0.009 

Teflon® 

2.21 ±0.04 

Below instrument capability 


Table 21 Parameters of the test-pieces and the capacitive sensors used in benchmark experiments. The three copper rods used as 
the conductive cores in the cylindrical test-pieces had a unifonn diameter of 15.90 ± 0.01 mm. 


Dielectric 

material 

/ (mm) 

2b (mm) 

Variation in 
(b-a) (mm) 

Po (mm) 


Tecaform™ 

20 

19.28 ±0.03 

0.06 

9.64 ±0.02 

172.4 ±0.3 

40 

19.22 ±0.05 

0.06 

9.61 ±0.03 

172.9 ±0.5 

Acrylic 

20 

19.18 ± 0.04 

0.03 

9.59 ±0.02 

173.3 ±0.4 

40 

19.17 ±0.07 

0.03 

9.59 ±0.04 

173.4 ±0.6 

Teflon® 

20 

19.31 ±0.02 

0.07 

9.65 ±0.01 

172.1 ±0.2 

40 

19.3 ±0.1 

0.07 

9.67 ±0.07 

171.9 ± 1.2 


Results and discussion 

Comparison between calculated and measured sensor capacitance C and dissipation factor D for each 
cylindrical test-piece and the two different electrode configurations is made in Table 22. Because of 
material complex permittivity, complex sensor output capacitance is expected. The sensor dissipation 
factor is defined as the ratio of the imaginary part of the capacitance to its real part. Since accurately 
measured imaginary permittivity was not achieved and Teflon®, the sensor dissipation factor cannot be 
calculated and comparison between its calculated and measured D is not made in Table 22. 
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Table 22 Measured and calculated capacitance for various sensor configurations in contact with different cylindrical test-pieces. 


Dielectric 

tube 

material 

Electrode 

length 

(mm) 

Calculated 
C( P F) 

Measured C 
(pF) 

Relative 
difference 
in C(%) 

Calculated 

D 

Measured 

D 

Relative 
difference 
in D 

Tecaform™ 

20 

6.58 

6.09 ±0.12 

-7.5 

0.003 

0.002 ± 0.002 

0.002 


40 

12.34 

11.47 ± 0.19 

-7.1 

0.003 

0.002 ± 0.002 

0.001 

Acrylic 

20 

5.64 

5.30 ± 0.1 1 

-6.0 

0.015 

0.014 ±0.006 

0.001 


40 

10.36 

9.47 ±0.17 

-8.6 

0.015 

0.015 ±0.007 

0.000 

Teflon® 

20 

4.20 

3.92 ±0.13 

-6.7 

- 

- 

- 


40 

7.63 

6.98 ±0.16 

-8.5 

- 

- 

- 


Measured and calculated capacitance C agree to within an average of 7.4% in Table 22. All the 
measurement results are smaller than the numerically-predicted ones, by between 0.3 to 0.9 pF. Two 
factors contributing to lower measured C are identified. First, during sample preparation, the brass rods 
were inserted into hollowed dielectric rods to achieve tight surface contact. After insertion, the dielectric 
coatings were subjected to normal stress exerted by the brass rods and circumferential strain was 
introduced in the dielectrics. In an independent study, it is found that the real permittivity of Teflon" 
(PTFE) decreases as its strain increases [129]. This means the actual permittivity of the Teflon" coating 
the brass rod is likely to be lower than the value used in the numerical calculations, which was measured 
on a sample slice cut from the rods prior to their being hollowed out, i.e., before strain was introduced. 
Consequently, measured C will be smaller than that predicted numerically for dielectrics with strain- 
induced reduction in permittivity. Second, although Kapton" films are tightly wrapped around the 
electrodes to attempt to eliminate air gaps between the electrodes and the test-piece, small gaps still exist 
due to the nonuniform diameters (surface roughness) of the cylindrical test-pieces. For example, the two 
largest percentage differences between calculated and measured capacitance values in Table 20 are 
observed for the cases of Acrylic and Teflon " coated test-pieces with 40-mm-long electrodes, which also 
show the largest variations in test-piece outer diameter (Table 20). It is seen from Figure 79 that the 

sensor output capacitance changes rapidly when the ratios b / p {) and a/b are lose to one. This indicates 

that for alb w 0.8 as in these measurements, the existence of air gaps can affect measurement results 
significantly. 

Apart from the above two factors, two sources contributing to measurement uncertainty (not necessarily 
lower C) exist. One is the misalignment of sensor electrodes. For example, the vertical edges of the 
electrodes tend to bend up, giving rise to small air gaps in the vicinity of inter-electrode gaps, where the 
sensor is most sensitive. This points to the fact that the ideal way to achieve best agreement between 
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theory and measurements is to deposit the sensor electrodes directly onto the cylindrical test-pieces, and 
errors coming from the misalignment of sensor electrodes and the existence of air gaps will be eliminated. 
However, deposition of electrodes directly onto the test-piece is costly, time-consuming and undesirable 
in some cases although may be useful for real-time monitoring of structures by in-situ sensors. The other 
source of uncertainty is variation in the dielectric tube thickness, Table 20. These variations can be traced 
to surface roughness of the dielectric tubes and/or non-concentricity between the conductive core and the 
dielectric tube. The effects of rough test-piece surfaces have been discussed above. The effects of 
eccentricity are discussed here. For the ideal case that the conductor and the dielectric are concentric, as 
shown in Figure 75, the total sensor capacitance my be approximated by a relationship of the following 
form: 


c = 


c 

c +c 


= C 0 /2, 


(94) 


where C 0 is the capacitance between the conductive core and either electrode. When the conductor and 

the dielectric become non-concentric, the capacitance between one electrode and the conductor changes to 

C 0 + mC 0 while the capacitance between the conductor and the other electrode changes to C 0 — fiC 0 , 

where w,n > 0 . The total capacitance C changes to 

c = (C 0 +/wC 0 )(C 0 -wC 0 ) = c l + m-n + mn 
C 0 + mC 0 + C 0 - nC 0 0 2 + m-n 

It is found that when m<n , C' is always smaller than C. However, when m > n , C' is not necessarily 
smaller than C. In other words, non-concentricity between the copper rods and the dielectric coatings does 
not necessarily results in smaller sensor output capacitance values, but does introduce uncertainty to the 
capacitance measurement. 

As shown in Table 22, the absolute difference between the numerically predicted and the measured 
sensor dissipation factor D is found to be 0.002 or less. This is close to the LCR meter measurement 
accuracy, e.g., the measurement accuracy of the dissipation factor of a 5 pF capacitance is ± 0.001 at 1 
MHz and room temperature. Large variation in measured D is observed, however. One important reason 

for this is that large variation (± 16%) in the measured test-piece 0. 2 is observed. Table 20, due to the 

measurement accuracy of the instrument. This fact also introduces uncertainty into the calculated D and 
contributes to the difference between numerically-predicted and measured results, given in Table 22. 
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Summary 

The electroquasistatic Green's function due to a point source exterior to a dielectric-coated conductor 
has been derived in cylindrical coordinates. The capacitance of the curved patch capacitive sensor has 
been calculated numerically using the method of moments based on the Green's function. The quantitative 
dependence of the sensor capacitance on test-piece geometry and the dielectric coating permittivity has 
been demonstrated numerically and verified experimentally. A discussion of measurement uncertainty is 
provided. The curved patch capacitive sensor developed in this section has the potential to evaluate 
effectively the condition of wiring insulation and is complementary to other inspection methods that are 
focused on evaluating the conductor's condition. 
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3 . A Capacitive Probe for Quantitative Nondestructive Evaluation of Wiring Insulation 

In this section, a capacitive probe is presented that has been developed for quantitative evaluation of 
wiring insulation permittivity. The probe consists of two patch electrodes that conform to the curvature of 
the wire under test. A previously-developed numerical model is utilized for inverse determination of 
insulation complex permittivity from measured probe response. Experimental studies on thermally and 
hydrolytically exposed wire samples show that the resulting insulation degradation is successfully 
detected using the described capacitive probe, for the wire type MIL-W-81381/12, which is 
predominantly Kapton ' coated. Permittivity changes in the wiring insulation detected by the capacitive 
probe are shown to be in accordance with results of research conducted previously on thermally and 
hydrolytically exposed Kapton” film samples. Thus the feasibility of quantitative evaluation of wiring 
insulation permittivity using capacitive techniques is demonstrated, which is of particular interest to the 
aerospace industry. 

Introduction 

Efforts have been made for years to guarantee the functionality of key electrical systems on aircraft. 
However, the perfonnance of the wiring that connects these key systems was not a strong focus of 
attention until the crashes of TWA 800 and Swissair 111 [88], attributed to aging wiring. In [89], causes 
and modes of failure in legacy aircraft wiring have been categorized. These causes include chemical 
degradation such as corrosion of current carriers and hydrolytic scission of polymer chains in the 
insulation; electrical degradation of the insulation that may be due to concentrated electric fields at sites 
of electrical stress; and mechanical degradation due to vibration, over-bending and other kinds of 
mechanical stress. 

Visual inspection is the most widely used method for aircraft wiring inspection. It is highly laborious 
while giving little quantitative information about the condition of the wires. Different physics-based wire 
inspection techniques have been developed over the past decade to replace this traditional inspection 
method, of which a summary is given here. 

One of the most commonly used physics-based techniques for the inspection of the conductor in aircraft 
wiring is reflectometry. A high frequency electrical signal is transmitted along the wire and any 
impedance discontinuities in the conductor result in reflected signals. The location of the fault can be 
determined from the time and/or phase delay between the incident and reflected signals. An excellent 
review paper that compares different reflectometry methods is [91]. Reflectometry methods are 
distinguished by the types of incident voltage applied. Time domain reflectometry (TDR) uses a short rise 
time voltage step as the incident voltage. This method is susceptible to noise and is not optimal for live 
wire testing [92] [93] [94], Frequency domain reflectometry (FDR) uses a set of stepped-frequency sine 
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waves as the incident voltage. A conceptual design of a smart wiring system based on FDR methods that 
can be used for on-board testing of aging aircraft wiring has been described in [90], Phase-detection 
frequency-domain reflectometry (PD-FDR) has also been applied for locating open and short circuits in a 
Navy F-18 flight control harness [95]. Sequence time domain reflectometry (STDR) and spread spectrum 
time domain reflectometry (SSTDR) use pseudo noise sequence and sine wave modulated pseudo noise 
code as the incident voltage, respectively. Testing systems based on these two techniques are capable of 
testing live wires and therefore have the potential to be used on energized aircraft to locate intermittent 
faults. Parameters that control the accuracy, latency, and signal to noise ratio for SSTDR in locating 
defects on live cables have been examined in [96], and the feasibility of spread-spectrum sensors for 
locating arcs on realistic aircraft cables and live wire networks has been demonstrated in [97] and [98], 
Aside from reflectometry methods, capacitive and inductive methods have also been applied for 
assessment of wiring conductor condition. In [99], linear relationships between the 
capacitance/inductance of open-/short- circuited wires (parallel insulated round wires, twisted-pair wires, 
and coaxial cables) and their conductor length have been demonstrated. This relationship enables the 
inverse determination of wire length from measured capacitance/inductance values. These techniques all 
inspect for so-called ‘hard' faults in the metal wire conductor itself and are not capable of inspecting the 
insulation conditions. 

Techniques have also been developed for the evaluation of wiring insulation condition. Infrared 
thermography and pulsed X-ray systems have been developed for nondestructive testing of wiring 
insulation [100], Infrared thermography offers the advantage of rapidly examining large areas of wiring 
and can serve as a global testing method, whereas a portable pulsed X-ray system can be used to obtain a 
radiographic image of a portion of the wire or cable. Ultrasonic methods have also been developed for 
quantitative assessment of degradation in wiring insulation condition caused by heat damage, by 
modeling insulated wires as cylindrical waveguides [101]. Moreover, acoustic and impedance testing 
methods aimed at locating intermittent faults in aircraft wires have been reported in [102], Partial 
discharge (PD) analysis methods for diagnosing aircraft wiring faults are explored in [104], in which a 
simulation of PD signal based on a high-voltage insulation testing standard [105] has been detailed, 
followed by wavelet-based analysis to de-noise the PD signals. Deficiencies of the above methods include 
the need of complex instruments in the measurement and not being able to provide quantitative 
information about the insulation condition at specific locations. A favorable solution to these deficiencies 
is the capacitive method, from which quantitative measurements of the permittivity of wiring insulation at 
specific locations can be made, from which its condition can be inferred, using relatively simple 
equipment. 
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This section describes a prototype capacitive probe that has been fabricated and applied for the 
inspection of wiring insulation condition of type M1L-W-81381/12 but in principle can be extended to 
any single-conductor wire. The probe is designed based on a previously developed physical model, in 
which a curved patch capacitor is exterior to and coaxial with a cylindrical dielectric that coats a 
conductive core cylinder (Chapter V Section 2). In the model, a quantitative relationship between the 
complex permittivity of the dielectric coating material and the capacitance and dissipation factor of the 
capacitor is described. This relationship is utilized here for quantitative assessment of wiring insulation 
permittivity, based on measured probe capacitance. To demonstrate the feasibility of this proposed 
technique, groups of wire samples have been thermally and hydrolytically exposed, under different 
conditions, to induce dielectric property changes in the insulation. Capacitance measurements were 
performed on the samples, and complex permittivity values of their insulation determined inversely by 
means of the numerical model. Comparisons made between the complex permittivity of the damaged 
wires and the control wires show that both the real and imaginary parts of the insulation permittivity of 
the damaged wires increase as the thermal exposure temperature/time and hydrolytic exposure time 
increase, and are higher than those of the control wires. Especially, changes in the imaginary permittivity 
are more significant than those in the real part. For example, following thermal exposure, the imaginary 
permittivity was observed to increase by up to 39% while the real part by up to 17%, for exposures at 
temperatures between 400 and 475 °C for various times up to 5 hours. In the hydrolytic exposure 
experiment, the imaginary permittivity was observed to increase by up to 60% and the real part by up to 
12%, for exposure in water for various times up to 4 days. Permittivity changes in the wire insulation 
detected by the capacitive probe are in accordance with results of independent measurements conducted 
previously on planar thermally and hydrolytically exposed Kapton® film (Chapter II). These proof-of- 
concept experiments have demonstrated the excellent capability of the capacitive probe for quantitative 
evaluation of insulation condition for wires of type MIL-W-81381/12, and in principle can be applied to 
any single-conductor wire. 

Apart from the capacitive probe discussed in this research, many other capacitive techniques have been 
developed and applied for NDE of dielectric materials [130]. For example, capacitive arrays have been 
developed for robotic sensing using ‘scanning' and ‘staring’ modes, [131]. Detailed literature surveys of 
capacitive NDE methods can be found in [1 16]and [125]. 

Summary of the physical model 

The capacitive probe is designed based on the sensor model described in (Chapter V Section 2). The 
curved patch capacitive sensor consists of two identical and symmetric electrodes, Figure 75. The wire 
under test is modeled as an infinitely long dielectric-coated perfect conductor. A relationship between the 
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complex permittivity s* of the dielectric coating, the sensor capacitance C and dissipation factor D is 
established in the electroquasistatic regime, in which s* = s' -j s' \ s' being the real permittivity, s" being the 

imaginary permittivity, and j = V-T . It is worth pointing out that, by adopting this model, the multilayer 
wire insulation Figure 84 is modeled as a one-layer structure. This means the insulation status assessed by 
the capacitive probe will be the overall condition of all the insulation layers. 

In the numerical model described in Section 2.2, the insulation real permittivity s' is in linear 
relationship with the sensor capacitance C, and the insulation loss tangent tan 6 = s'7 s' is linearly related 
to the sensor dissipation factor D. In this section, complex permittivity values of wire samples are inferred 
from measured probe capacitance and dissipation factor values, based on this numerical model. 

In the model, the conductor is kept at ground potential and potentials on the two electrodes are kept at 1 
V and -IV respectively. The calculated capacitance is the series capacitance of the capacitance between 
one electrode and the conductive core, with the capacitance between the conductive core and the other 
electrode. In the measurements, however, the conductive core does not have to be grounded. This is 
because the potential on the conductive core will be the average of that on the two electrodes, and the 
measured capacitance under this circumstance is the series capacitance calculated through the model. 

The validity of this physics-based model has been verified by good agreement between numerical 
calculations and results of benchmark experiments (Chapter V Section 2). A detailed description of 
numerical modeling and experimental verification can be found in (Chapter V Section 2). 

Probe and measurement system 

Probe fabrication 

Figure 82 shows the capacitive probe. The probe consists of two 2x4 cm 2 acrylic plates and an acrylic 
rod on which the two plates are mounted so that their surfaces remain parallel. The lower plate is attached 
to the acrylic rod using a plastic screw, whereas the upper plate can glide up and down by adjusting 
another plastic screw perpendicular to the two plates. The curved sensor electrodes are formed by 
brushing a layer of silver paint onto the symmetric grooves in the two plates. The two electrodes are 
connected to two pins, which are then connected to an LCR meter for capacitance measurements. In order 
to ensure that the two plates are in parallel, two acrylic dowels are attached to the upper plate using epoxy 
glue, and inserted into the lower plate. These two acrylic dowels, together with the plastic screw, assure 
that the two acrylic plates remain parallel during measurements. The subfigure in Figure 82 shows the 
probe holding a wire under test. The probe and the wire are in tight surface contact with each other. 
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Figure 82 Photograph of fabricated capacitive probe with curved patch electrodes. Subfigure: capacitive probe holding a wire 

sample under test. 

Parameters 

Parameters of the probe are (see Figure 75): electrode radius p (] - 2.50 mm and electrode length / = 20 

mm. The electrode radius is taken to be the same as the specified wire sample outer radius. The electrode 
length is measured directly from the probe. 

The region exterior to the curved electrodes in Figure 75 is assumed in the model to be homogeneous, 
e.g., air. This assumption is not satisfied, however, for the probe depicted in Figure 82 due to the 
existence of the acrylic plates exterior to the electrodes. To account for this effect, along with possible air 
gaps between the electrodes and the testing wire due to surface roughness, an effective electrode arc -angle 

<f> 0 is introduced for the capacitive probe. Steps taken to determine (j> 0 are shown in Figure 83. First, an 

assumed electrode arc -angle </> 0 , the wire sample geometry information, and the assumed wire insulation 

real permittivity ( s / = 2.7) are input into the numerical model, from which a computed probe capacitance 
C mo dei is obtained. Second, capacitance measurements are performed on the virgin wire samples. Third, 
C m odei is compared to the measured probe capacitance C meas . If C modd and C meas agree to three significant 

figures, </> 0 is considered as the effective electrode arc -angle. Otherwise, a different value is assumed for 
the effective electrode arc-angle and the above steps are repeated until the stop criterion is satisfied. For 
testing wires of type MIL-W-81381/12, the effective arc-angle is determined as (f> 0 = 80.1 ± 0.5°. This 

inferred (f is quite close to the physical electrode arc -angle, which is in the range 80 to 90°. Uncertainty 
in (f>Q is due to variations in the measured capacitance (2.1 1 ± 0.01 pF). Probe parameters are tabulated in 


Table 23. 
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Table 23 Probe parameters. Figure 82. 


/ (mm) 

20.0 ±0.06 

Po (mm) 

2.50 ±0.13 

(/) 0 (degrees) 

80.1 ±0.5 


Select 
effective 
electrode 
arc-angle <p0 



Do not 
agree? 


Calculated sensor 
capacitance 

Compare 


Measured probe capacitance 
on a virgin wire 


Agree 


(pO is the 
effective 
arc-angle 


Figure 83 Algorithm for determination of the effective electrode arc-angle . 


Measurement system and uncertainty analysis 

For capacitance measurements, the probe was connected to an Agilent LCR meter 4980A by an Agilent 
probe test fixture 16095A, whose probe pins were connected to the curved patch electrodes of the 
prototype probe. Capacitance measurements in this section were performed at 1 MHz and room 
temperature. Detailed discussion on the selection of operating frequency is provided in (Chapter V 
Section 2). 

Uncertainties in the measurement procedure are attributed to achievable precision in the measurement 
equipment, uncertainty in (f > 0 , and variation in the geometry of the wires under test. These uncertainties 

contribute to uncertainty in the insulation complex permittivity determined inversely utilizing the 
numerical model. 

The achievable precision of the LCR meter depends on the magnitude of the capacitance measured. In 
the LCR meter operation manual, relative uncertainty is provided for the measured capacitance whereas 
absolute uncertainty is provided for the measured dissipation factor. Measured capacitance values in this 
paper are within the range 2 to 3 pF and for these values the corresponding uncertainty is 0.18% for the 
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capacitance and 0.0018 for the dissipation factor. This means that, in order to detect insolation 
degradation by capacitive measurements, the degradation-induced changes in the measured probe 
capacitance and dissipation factor must be that greater than these values. 

As seen in Table 23, the uncertainty in p 0 and / are approximately 5%. When the arc -angle of the 
curved patch electrodes is far from 180°, i.e., the distance between the edges of the oppositely-charged 
electrodes is relatively large, the capacitance and dissipation factor change slowly as (j> () varies (Chapter V 

Section 2). This is because the interaction between the electrodes is not intense under this circumstance. 
For the degradation cases studied in this paper, the uncertainty in the inferred insulation real permittivity 
s r ' resulting from the uncertainty in (j) {] is within 0.01, i.e., less than 0.5%, whereas that in s r " is 0.0002, 
i.e., less than 0.2%. 

Real wires exhibit variations in their dimensions, surface roughness, roundness and curvature even 
when they may appear macroscopically similar. For this reason, in the experiments that follow, 
capacitance measurements were performed on three samples, for each degradation condition. Error bars 
are included in all the following measurement results. It is worth pointing out that the outer diameters of 
the thermally and hydrolytically exposed wire samples still lie in the range between 2.41 and 2.63 mm. As 
will be seen in Section 3.6, the standard deviation in the measurements deriving from these sources is the 
dominant source of uncertainty in the measurements, larger than that due to the other sources discussed 
above. Nonetheless, changes in the insulation complex permittivity, due to thermal and hydrolytic 
exposure, are clearly observed even when uncertainty is taken into account. 

Parameters of the wire under test 

Capacitance measurements throughout this paper are performed on aircraft wire samples of type M1L- 
W-8 138 1/12. The wire is composed of a nickel-coated copper conductor, wrapped with two layers of 
polyimide 150FW-N019 film and one layer of aromatic polyimide coating. Each layer of the polyimide 
150FWN019 film is constructed by attaching a 13- // m-thick fluorinated ethylene propylene (FEP) 
fluoropolymer film to a 25- p m-thick polyimide FN film, Figure 84. Nominal conductor and outer 
diameters of the wire are 2.09 and 2.50 mm, respectively. Actual measured wire outer diameters for wires 
examined in this study range between 2.41 and 2.63 mm. 
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Figure 84 Schematic diagram of fluorinated ethylene propylene (FEP)-coated polyimide 150FWN019 film. Nominal thicknesses 
of each layer of polyimide 150FWN019 film, FEP fluoropolymer film, and the liquid H lacquer film are 25 pm, 13 pm, and 129 
pm, respectively. The FEP provides adhesion between the layers of polyimide. Z>c=2.09 mm, Z>w=2.50 mm (nominal). 


Real permittivity 

In order to compare the inferred permittivity of the thermally/hydrolytically exposed wires to that of the 
control wires and therefore provide quantitative assessment of the condition of the insulation, an initial 
value for the permittivity of the control wires must be assigned. 

The real part of the undamaged wire insulation permittivity is assumed to be 2.7 at 1 MHz and room 
temperature, the condition under which capacitance measurements were performed. Reasons for this 
assumption are as follows. Considering the schematic shown in Figure 84, we note that there are two 
main components: polyimide 150FWN019 and an aromatic polyimide outer layer. The manufacturer, 
DuPont, does not provide a specific permittivity value for polyimide 150FWN019 films, mentioning only 
that it is less than 3 at 1 kHz and room temperature. According to DuPont, the processing conditions of 
polyimide 150FWN019 film is, however, very similar to those of the standard polyimide FN films, whose 
real permittivity is 2.7 at 1 kHz and room temperature. Therefore, we assume the real permittivity of 
polyimide 150FWN019 to be 2.7 at 1 MHz and room temperature, taking into account the fact that the 
real permittivity of polyimide does not change significantly over the frequency range from 1 kHz to 1 
MHz. Regarding the aromatic polyimide portion of the insulation, the dielectric behavior of aromatic 
polyimide, types a through/ is presented in [132], As will be seen in Section 3.6, degradation of the wire 
insulation initiates between 400 °C and 450 °C. Only type / of the aromatic polyimides studied in [132] 
presents a thermal degradation initiation temperature (T; = 430 °C) within this range, whereas Tj of the 
other types are greater than 450 °C. The aromatic polyimide coating used in the wire insulation is 
therefore inferred to be type/ whose measured real permittivity is 2.7 at 1 MHz and room temperature 
[132], the same as that assumed for the polyimide 150FWN019 films. Finally, since the FEP adhesive 
layer is a relatively minor constituent of the insulation, accounting for approximately 13% of the 
thickness, and exhibits no unusual dielectric properties (s' =2.01 at 1 MHz and room temperature 
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according to DuPont), it is reasonable to assume that overall real permittivity of the undamaged insulation 
is 2.7 at 1 MHz and room temperature. 

Imaginary permittivity 

In order to compare the complex permittivity of thermally and hydrolytically exposed wires to the 
control/virgin wires, imaginary permittivity s' for the control wires has to be determined. Figure 85 shows 
the steps taken to determine the imaginary permittivity s' for the virgin wires. These steps are similar to 
those described in Figure 85, except that the quantity compared here is the probe dissipation factor D, 
instead of the capacitance C. The inversely determined s r " for the control wires is 0.016 ± 0.002. 
Uncertainty in s r " is due to variations in the measured dissipation factor of the control wires (0.0055 ± 
0.0006) and </> 0 . The inferred s r " is in accordance with the results for type / aromatic polyimide given in 

[132] ( s r " ~ 0.01 at 1 kHz and room temperature), with the value for polyimide 150FWN019 given by 
DuPont {s," ~ 0.006 at 1 kHz and room temperature), and with the value for FEP adhesive given by 
DuPont (s r " ~ 0.004 at 1 MHz and room temperature). 



Figure 85 Algorithm for determination of the imaginary permittivity s' for the virgin wires. 

Case study: evaluation of polyimide-coated wires after thermal and hydrolytic exposure 

The influence of thermal degradation and saline exposure on the complex permittivity of polyimide HN 
films has been studied in (Chapter II). The samples studied were 125-jt/m-thick. Explanations of how 
thermal degradation and saline exposure affect the complex permittivity of polyimide are provided. The 
work of (Chapter II) guided the choice of experimental parameters for this study. 

For the thermal exposure experiment, the exposure temperatures were selected as 400, 425, 450 and 
475 °C. For each exposure temperature, five groups of wires, each with three 4-cm-long samples, were 
isothermally heated in a muffle furnace, Figure 86, for 1, 2, 3, 4 and 5 hours. 
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Note, the temperature distribution in the furnace is not uniform. To ensure that all the samples are 
heated at the selected temperature in each experiment, a small ceramic bowl was used to accommodate 
the wires, and the temperature in the ceramic bowl was measured independently using a thermometer. 
The upper right figure in Figure 86 shows some of the heat-damaged wire samples in comparison with a 
control wire. Thermally exposed wires were sealed in plastic bags immediately after being taken out of 
the furnace to avoid moisture absorption. These samples were then cooled down to room temperature. 

The lower right figure in Figure 86 shows a wire sample used in the hydrolytic experiment. Both ends 
of the wire were sealed with wax to prevent water from migrating into the samples from its ends and via 
the conductor insulator interface. Five groups of wires, each having three 4-cm-long samples, were 
immersed in water at room temperature for 0.5, 1, 2, 3, and 4 days, respectively. Capacitance 
measurements were performed immediately after the samples were removed from the water. Sample 
surfaces were wiped dry with a soft cloth before capacitance measurements. 



Figure 86 Left: muffle furnace used for thermal exposure. Upper right: wire samples after heat exposure (brown) and a control 
wire (yellow). The samples are 4 cm long. Lower right: For hydrolytic exposure, both ends of the sample are sealed with wax. 

Results and discussion 

Thermal exposure 

Figure 87 shows probe capacitance and dissipation factor measured on the thermally exposed wires, as 
a function of exposure temperature and time. Measurements were performed on the three normally 
identical wire samples in each group and the mean value and standard deviation obtained. 

As can be seen in Figure 87, measured capacitance increases as exposure time and temperature increase. 
A similar behavior was observed in the measured dissipation factor. In Chapter II, it was found that the 
dissipation factor of polyimide HN films did not change significantly unless exposed at 475 °C for more 
than 3 hours. This suggests that increases in the dissipation factor of the thermally degraded wires 
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observed here results largely from degradation of the aromatic polyimide coating, rather than the 
polyimide 150FWN019 layers. 



Thermal exposure time (hrs) 



Thermal exposure time (hrs) 


Figure 87 a) measured capacitance and b) dissipation factor for heat exposed wires. Uncertainties derive from the standard 
deviation of measurements on three separate samples. Physical degradation of the sample heated beyond 2 hours at 475 °C 
prevented accurate capacitance measurement for those conditions. 

Complex permittivity of the wire samples was inferred from the measured probe capacitance and 
dissipation factor in the following way. First, an initial guess of the sample complex permittivity s* is 
input into the numerical model, from which particular values of probe capacitance C ca i c and dissipation 
factor D ca!c are obtained. These values are compared with the measured values C meas and D meas and s* 
adjusted until C cafc ( D cak: ) and C meas (D meas ) agree to within three significant figures. Then, e is considered 
as the inferred sample complex permittivity. 
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Figure 88 shows the inferred insulation real permittivity values in comparison with results presented in 
Chapter II. The way in which s' increases with time and temperature of heat exposure, Figure 88, is in 
very good agreement with results presented in Chapter II. The inferred imaginary permittivity values of 
the wires are presented in Figure 89, in which the existence of aromatic polyimide coating contributes 
predominately to the observed increases in the insulation imaginary permittivity. As can be seen from 
these two figures, both the real and imaginary parts of the insulation increase as heat exposure 
temperature and time increase. Especially, the relative change in the insulation imaginary permittivity is 
greater than that in the real permittivity; the imaginary part increases by up to 39% and the real part by up 
to 17% compared to the control wires. 



Figure 88 Inferred real permittivity s' of the thermally exposed wires in comparison with that of polyimide HN film (Chapter II). 

a) 400 and 425 °C; b) 450 and 475 °C. 
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Thermal exposure time (hrs) 

Figure 89 Inferred imaginary permittivity e" of the thermally exposed wires. 

Hydrolytic exposure 

The capacitance and dissipation factor for hydrolytically exposed wires are shown in Figure 90. In 
accordance with results of previous studies [53] [133], it is observed that both the measured capacitance 
and dissipation factor increase as water immersion time increase in the first three days, and do not change 
significantly afterwards. Inferred wire insulation complex permittivity is shown in Figure 91 along with 
measurement results presented in Chapter II for 125 -/un-thick polyimide HN film. The real part of the 
wire insulation permittivity increased by up to 12% and the imaginary part by up to 60% compared to the 
control wires. As can be seen from Figure 91 permittivity change in wire insulation as a function of 
hydrolytic exposure time is in accordance with the results of independent measurements conducted 
previously on hydrolytically exposed planar polymer insulation samples. 


137 




Figure 90 As for Figure 87 but for hydrolytically exposed wires. 
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Water immersion time (days) 



Figure 91 Inferred complex permittivity, a) real part and b) imaginary part, of the hydrolytically exposed wires in comparison 

with that of polyimide HN film (Chapter II). 


4. Summary 

Complex permittivity is an effective indicator of wiring insulation condition. Changes in wiring 
insulation permittivity, induced by thermal degradation and hydrolytic exposure, have been successfully 
detected using the capacitive probe presented here, and observations are in accordance with results of 
previous research (Chapter II). Furthermore, the one-to-one correspondence that exists between the 
insulation permittivity and measured capacitance justifies the monitoring of changes in insulation status 
by detecting changes in insulation permittivity, even though the absolute value of the control wire 
permittivity is assumed from the literature rather than directly measured. 

This section focuses on quantitative evaluation of the insulation permittivity after degradation. In 
practice, it may be more convenient under some circumstances to make comparative measurements of the 
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capacitance measured on a reference wire and that on the wire under test. The probe presented here is 
capable of both tasks. Additionally, prior knowledge of the locations at which wiring insulation 
degradation is most common in any particular aircraft is usually available. The capacitive approach 
described in this section can therefore be applied for accurate characterization of localized insulation 
degradation. 

For practical measurements, a capacitive device similar to this one can be mounted on a mechanical clip 
which can be simply applied to an exposed wire. This is the subject of future work. 
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Chapter VI. Conclusion 

The research presented in this report is motivated by current concerns for aging of electrical wiring 
insulation materials used in space/air-crafit. The influence of various environmental aging processes on 
insulation capabilities of three insulation polymers has been investigated. Also, a capacitive 
nondestructive sensor is developed to detect flaws in wiring insulation. 

Chapters II, III and IV of the report present the investigation of the deterioration of insulating function 
of three insulation polymers in association with changes in their material characteristics during various 
environmental aging processes. The findings of these investigations have successfully answered questions 
related to how electrical signatures of insulation material change over time with respect to environmental 
aging, and can be used to develop a library containing relationships between material characteristics and 
insulation performances. Thermal degradation process of the three polymers has been verified by kinetic 
reaction models. 

Thermal exposure in air and immersion in water/saline of PI has been found to lead to significant 
deterioration of insulating performance of PI, Chapter II. After either the heating process or immersion in 
water/saline, the real permittivity and loss factor of PI are increased substantially, there is a significant 
decrease in dielectric breakdown strength and a higher density of electrically weak points is observed. 
These observations are explained in terms of chemical degradation due to pyrolysis of imide groups 
during thermal exposure and formation of ionic side groups during the hydrolysis process of PI while 
immersed in water. However, dissolved sodium chloride shows no significant influence on dielectric 
properties of PI immersed in saline solutions. 

Results of a study on the influence of thermal exposure and tensile strain on permittivity of PTFE are 
presented in Chapter III. An increase of approximately 2% in the real permittivity of PTFE is observed 
following isothermal heating at 340 °C for 96 hr, due to associated increase in crystallinity that enhances 
the interfacial polarization between the amorphous and crystal phases of the polymer. On the other hand, 
the real permittivity of PTFE was observed to decrease by approximately 19% as a consequence of 150% 
mechanical strain. The observation is attributed to the more ordered structure which arises due to 
alignment of the polymer chains along the direction of loading, which limits motion of polar groups and 
weakens the interfacial polarization. An increase in the degree of crystallinity was also observed in this 
case, but it should be noted that value of crystallinity is not the primary factor in determining interfacial 
polarization relaxation intensity (and therefore permittivity). Rather, the nature and degree of perfection 
in the crystallized regions play a more dominant role in determining the dielectric properties of the 
polymer. The degradation mechanisms studied here (thermal exposure and mechanical stress) give rise to 
crystallizations that are significantly different from each other in nature and degree of perfection. Thus, 
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the trends in the measured values of crystallinity and in the permittivity for each of the degradation 
mechanisms described above are not in conflict but rather explained by significant differences in the 
crystal nature of those samples. 

Results of a study on the effect of thermal exposure at 160 °C for 96 hr in air on permittivity of ETFE 
are presented in Chapter IV. The 2% increase in the real permittivity and 11% decrease in the loss factor 
of ETFE are associated with oxidation and dehydrofluorination of ETFE during heating, which give rise 
to formation of polar groups and enhance polarization in the polymer in the presence of an electric field. 

Chapter V described the development of a prototype capacitive probe for quantitative NDE of wiring 
insulation. A numerical model, based on the electrostatic Green's function due to a point source exterior 
to an infinitely long cylindrical dielectric, and exterior to a dielectric-coated cylindrical conductor, was 
developed to evaluate quantitatively the permittivity of the dielectric in these cases. The capacitance of a 
capacitive sensor with two opposing curved patch electrodes was calculated numerically using the method 
of moments based on the Green's function. The quantitative dependence of the sensor output capacitance 
on the test-piece permittivity and radius was demonstrated numerically and verified experimentally. The 
permittivity of various cylindrical test-pieces was inferred from measured capacitance to within 1% 
accuracy for pure dielectric cylinders. For dielectric-coated conductors, numerical calculations of 
capacitance and measured agreed to within 7%, on average. 

A prototype capacitive probe for quantitative NDE of wiring insulation was designed and fabricated, 
based on the numerical model, and the insulation complex permittivity was determined inversely from 
measured probe response. Thermal and hydrolytic exposures were used to induce changes in insulation 
status for groups of wire samples. Experimental studies on both the damaged and undamaged wire 
samples demonstrate that insulation status changes for wire type M1L-W-8 138 1/12 can be successfully 
detected and quantified using the capacitive probe described in this research. In principle, the same 
technique can be applied to evaluate other wires composed of dielectric -coated conductors. 
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